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SUMMARY 


This  Final  Report  under  Contract  No.  F19628-82-C-0060  summarizes  the  results  of  several  research  studies  directed 
toward  the  development  of  models  of  the  optical  properties  of  the  atmosphere  which  are  suitable  for  the  operational 
support  of  tactical  weapons  and  sensors  operating  in  the  visible  and  infrared  portions  of  the  electromagnetic  spectrum. 

The  overall  contractual  effort  has  been  an  integral  part  of  the  overall  Air  Force  Geophysics  Laboratory 
(AFGL/AFSC)  program  to  develop  a  thorough  understanding  of  atmospheric  effects  on  visible  and  infrared  radiation 
through  the  atmosphere,  and  ultimately,  the  concomitant  limitations  on  the  performance  of  electro-optical  systems  such  as 
target  acquisition  devices  and  precision  guided  munitions. 

This  Final  Report  contains  a  brief  review  of  four  previously  issued  interim  reports  which  were  submitted  during  this 
contractual  interval  in  support  of  the  primary  model  development  task.  It  also  contains  the  complete  documentation  for  “ 
the  FASCAT  code  and  a  discussion  of  the  FASCAT  model  performance  under  fractional  cloud  meteorological  conditions. 
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THE  FASCAT  MODEL  PERFORMANCE  UNDER  FRACTIONAL  CLOUD  CONDITIONS 

AND  RELATED  STUDIES 

by 

Wayne  S.  Hering  and  Richard  W,  Johnson 


1.0  INTRODUCTION 

This  Final  Report  has  been  prepared  under 
Air  Force  Geophysics  Laboratory  Contract  No.  F 19628 - 
82-C-0060  and  addresses  several  studies  related  to  the 
development  of  models  Tor  atmospheric  optical  properties 
within  the  visible  and  infrared  portions  of  the  electromag¬ 
netic  spectrum. 

Initial  research  studies  and  their  subsequent 
reported  Andings  are  reviewed  In  section  two  of  this  Final 
Report.  Section  3,  the  major  element  within  this  report,  is 
reserved  for  a  discussion  of  the  FASCAT  (Fast  Atmos¬ 
pheric  SCATtering)  model  in  its  current  and  most  sophis¬ 
ticated  conAguration.  This  model,  Arst  presented  in  Her¬ 
ing  (1981)  has  been  extended  to  address  singly  and  multi¬ 
ply  scattered  radiances  not  only  under  clear  and  fully  over¬ 
cast  meteorological  conditions,  but  under  the  more 
difficult  and  substantially  more  general  conditions  typiAed 
by  fractional  cloud  cover. 

2.0  A  REVIEW  OF  RELATED  STUDIES 

There  were  several  research  studies  related  to  visible 
and  infrared  optical  properties  within  the  lower  tropo¬ 
sphere  that  were  conducted  as  part  of  this  overall  contrac¬ 
tual  effort.  The  subject  areas  ranged  from  the  generalized 
development  for  equation  of  transfer  formalisms,  Gordon 
(1982),  to  the  speciAc  analyses  of  experimentally  meas¬ 
ured  aerosol  size  distributions,  Fitch  (1983). 

In  order  to  briefly  summarize  these  analysis  and 
modelling  efforts,  four  of  the  resultant  scientiflc  reports 
which  were  Issued  during  this  contract  interval  are 
reviewed  in  the  following  paragraphs.  In  general,  the  for¬ 
mat  will  contain  a  very  brief  intoductory  comment,  fol¬ 
lowed  by  an  abstract  of  the  original  report's  primary  result 
and/or  recommendation. 

2.1  AFGL-TR-82-0223:  Implications  of  the  Equation  of 
Transfer  Within  the  Visible  and  Infrared  Spectrum 

by  Jacqueline  I.  Gordon 


Note:  Subsequent  pre-publication  review  has  modifled  the 
derivation,  but  not  the  results  discussed  in  Section  2.1  as 
related  to  Eq.  4.17  through  4,22.  See  App.Opt.  Vol.  24, 
No.  16. 


Introductory  Comment ;  As  documented  in  a  series  of 
preceding  AFGL  reports  which  are  tabulated  in 
Appendix  C,  the  Visibility  Laboratory  has  conducted  a 


substantial  program  of  visible  spectrum  measurement, 
analysis,  and  modeling.  The  radiative  transfer  formalisms 
used  in  these  programs  have  been  visible  spoctrum 
oriented,  and  in  general  did  not  contain  emission  and 
absorption  terms.  This  report,  -0223,  extends  the  work  by 
including  these  terms  in  the  analytic  formalisms  and  then 
uses  these  expressions  to  re-evaluate  previously  obtained 
experimental  measurements  and  assesses  their  use  In 
speciAc  modelling  applications. 

EQUATION  OF  TRANSFER 
(Section  2,  AFOL-TR-82-0223) 

Radiance 

The  basic  equations  developed  herein  are  mono¬ 
chromatic.  Applicability  to  broadband  sensors  will  be  dis¬ 
cussed  where  appropriate. 

The  most  basic  equation  in  radiation  theory  is  the 
equation  of  transfer  (Eq.  10  Duntley  etal.  (1937), 
Eq.  (3)  Sec.  3.15  Prelsendorfer  (1976),  and  Eq.  (46) 
Chandrasekhar  (I960)) 

dL(z,0,4)/</r  -  -  a(z)L(z,0,*)  +  L«(z,  0,6).  (2.1) 


This  equation  relates  the  incremental  change  in  radiance 
rfL(z,0,$)  at  altitude  z  in  direction  zenith  angle  0  and 
azimuth  6  over  the  incremental  path  length  dr  to  the 
attenuation  coefficient  a(z),  the  radiance  L (z,  0,  <t>)  and 
the  path  function  L,(z,0,$).  The  Arst  term  is  the  loss 

term  and  the  second  is  the  gain  term. 

The  attenuation  or  extinction  coefficient  a(z)  is 
equivalent  to  the  Chandrasekhar  (1960)  mass  attenuation 
coefficient  k  times  the  density  p.  The  attenuation 
coefficient  is  the  sum  of  the  total  scattering  coefficient 
s  (z)  plus  the  absorption  coefficient  a  (z) 

a(z) -  r(z)  +  o(z),  (2.2) 

The  path  function  L*(z,0, 0)  is  equivalent  to  the 

Chandrasekhar  (1960)  mass  emission  coefficient  J  times 
the  density  p.  The  path  function  is  the  sum  of  a  scattered 
component  and  an  emitted  component  which  is  related  to 
the  absorption 

f.,(z,0,+)- L.  <r,0,0>  +  L,  (z).  (2.3) 

i  o 
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The  emitted  component  due  to  absorption  is  isotropic  and 
hence  is  shown  without  direction  modifiers. 


Scattering 

The  scattering  coefficient  has  two  components.  Ray¬ 
leigh  or  molecular  scattering  Rs  is  highly  wavelength 
dependent  being  proportional  to  the  inverse  fourth  power 
of  the  wavelength  \~4.  Mie  scattering  (due  to  atmos¬ 
pheric  particulates)  is  less  wavelength  dependent.  Both 
tend  to  be  smooth  continuous  functions  with  wavelength. 

The  path  function  due  to  scattering  is  the  Integral  of 
the  incoming  radiance  l(z,0\  ^')  in  all  4»r  directions 
(including  the  sun  where  appropriate)  times  the  volume 
scattering  function  cr(z,0‘) 


Ltf  CM,  *)  -  Or  ,#',  *')  <r  (: &)dn ,  (2,4) 

The  scattering  angle  0'  is  equivalent  to  the  Chandrasekhar 
(1960)  angle  0.  It  is  a  function  of  the  incoming  radiance 
direction  angles  9\  <j>'  and  the  sensor  direction  angles  9,  4> 
as  follows 


cos0'  -  sin#  sin#'  cos  +  cos#  cos#'  (2.5) 

The  Integral  of  the  volume  scattering  function  over  4  n  is 
the  total  scattering  coefficient 


The  volume  scattering  function  Is  equivalent  to  the  Chan¬ 
drasekhar  (1960)  phase  function  p(cos0)  times  the 
attenuation  coefficient  «  divided  by  4  rr. 

Absorption  and  Emission 

Atoms  and  molecules  in  the  gas  phase  absorb  and 
emit  in  line  and  band  spectra.  The  atomic  spectra  are  line 
spectra  and  tend  to  be  at  the  shorter  wavelengths.  The 
band  spectra  are  molecular  and  tend  to  be  at  the  longer 
wavelengths.  The  continue  are  essentially  part  of  the  band 
spectra,  weak  but  broad  spectrally. 

Molecules  In  the  liquid  or  solid  phase  cannot  emit 
line  or  band  spectra  but  can  only  emit  in  a  continuous 
spectrum,  the  distribution  of  which  is  determined  by  the 
ambient  temperature,  in  other  words  like  a  black  or  grey 
body  (incomplete  radiator).  This  implies  that  the  absorp¬ 
tion  spectrum  for  liquids  and  solids  is  also  continuous 
spectrally. 

An  aerosol  is  defined  as  a  mixture  of  gas  and  small 
solid  or  liquid  particles.  Although  the  solid  particulates 
and  water  droplets  with  dissolve^  particulates  absorb  and 
emit  like  a  grey  body  with  a  continuous  absorption  spec¬ 
trum,  the  gas  or  air  absorbs  and  emits  in  line  and  band 
spectra.  Thus  the  aerosol  absorbs  and  emits  like  a  spectral 
or  colored  body, 


Emission  Mechanisms.  The  principal  emission 
mechanisms  above  the  mesosphere  are  electrolumines¬ 
cence  and  chemiluminescence  resulting  in  line  and  band 
spectra.  These  emissions  are  called  airglow.  These  are 
Important  In  the  visible  part  of  the  spectrum  at  twilight 
and  night  but  will  not  be  dealt  with  herein. 

The  principle  emission  mechanism  in  the  tropo¬ 
sphere  and  in  the  atmosphere  at  or  below  the  mesosphere 
is  temperature  radiation  which  is  photon  emission  caused 
by  atomic  or  molecular  collision.  The  atmosphere  is 
assumed  to  be  in  local  thermodynamic  equilibrium,  hence 
KlrkhofTs  law  applies.  The  path  function  due  to  emission 
is  thus  [Eq.  (38)  Chandrasekhar  (I960)] 

L,  (z)-a(r)4(X,D  (2.7) 

0 

where  L(X,D  is  the  blackbody  radiance  at  wavelength  X 
and  T  is  the  temperature  in  degrees  Kelvin.  The  path 
function  due  to  emittance  at  and  below  the  mesosphere  is 
highly  wc  :  .igth  dependent  and  should  be  dealt  with 
monocht  natically, 

A  black  body  In  thermodynamic  equilibrium  absorbs 
and  emits  as  a  continuous  function  of  wavelength  and 
temperature  according  to  the  classical  equation  [Wolfe 
(1978)  Table  1.7] 


£(x,n-c,/  [irX5(**-D]  (2.8) 

where 

c,  -  2ir  he2  -  3.741382E-I6  Wm2  (2.9) 

and 

x-ci/an.  (2.io) 

The  form  3.74E-16  Is  an  alternate  format  for  3.74x  10~“. 
This  computer  form  Is  used  throughout  this  report.  The  c 
is  the  speed  of  light,  h  is  the  Planck  constant  and 

ca  -  ch/k  -  1.438786E-2  mK  (2.11) 

where  k  is  the  Boltzmann  constant.  The  constants  are 
from  Driscoll  and  Vaughn  (1978)  Table  A.l. 

Visible  Spectrum.  Blackbody  radiance  and  hence  exi- 
tance  (emittance)  is  negligible  in  the  visible  spectrum. 
Even  at  lMw  at  300°K  (27°C  or  80*F), 

Z,  Om*. 300*)-  I.76E-7  W/flm2,  At  shorter 
wavelengths  (all  the  visible  wavelengths)  and/or  lower 
temperatures  (the  normal  range  of  temperature  below  the 
mesosphere)  the  emittance  is  still  less.  Hence  for  the  visi¬ 
ble  spectrum  Eq.  (2.3)  becomes 

L,(r,#,«t)  -  m*,#.*)  .  (2.12) 
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Equilibrium  Radiance 


U  {*.»,,  4) -f  L,<*,,0s,4)Trl(z,9,)dr  .  (3.1) 


The  equilibrium  radiance  is  defined  as  the  radiance 
when  the  incremental  change  in  radiance  over  the  incre¬ 
mental  path  length  dL/dr  equals  zero,  therefore  from 
Eq.  (2.1)  (Eq.  (11)  Duntleyefa/.  (1957)) 

0- L,(z,0,<p)- a(z)  Lq(z,0,<t>)  (2.13) 

or  (Eq.  (40)  Chandrasekhar(1960)) 


L?(z,0,4)-L,(z,fl,<i)/a(z).  (2.14) 


Thus  the  equilibrium  radiance  is  equal  to  the  Chan¬ 
drasekhar  (1960)  source  function  7. 

No  Absorption.  In  a  medium  or  at  wavelengths  with 
no  absorption,  the  attenuation  coefficient  would  equal  the 
scattering  coefficient  and  the  equilibrium  radiance  would 
be  equal  to 


L,(z,fl,*)- (.^(z.fl^lMz).  (2.15) 


Scattering  Negligible.  When  scattering  is  negligible, 
such  as  at  longer  wavelengths  in  the  infrared,  the  attenua¬ 
tion  is  equivalent  to  the  absorption.  Substituting  Eq.  (2.7) 
Into  Eq.  (2.14)  we  have 


Z.,(i,0,*)-£(x,n.  (2.16) 

Equation  of  Transfer,  An  alternate  form  of  the  equa¬ 
tion  of  transfer  is  to  divide  both  sides  by  the  attenuation 
coefficient  and  let  dr  equal  sec  6  dz  (Eq.  63  Chandrasekhar 
(I960)], 


cosfl  dL  (2,9.  <t>)  /  [a(r )  di\ 

--Z.<z,fl,*)  +  L,(z,S,*).  (2.17) 

The  cosfl  Is  equal  to  the  Chandrasekhar  (1960)  function  p 
and  «(z)  dz  is  equal  to  dt(z)  the  differential  of  the  optical 
thickness  r(z). 


DETERMINATIONS  FROM 
THE  SOLAR  ALMUCANTAR 
(Section  3.  AFGL-TR-82-0223) 

Sky  Radiance 

The  almucantar  is  the  part  of  the  sky  at  a  constant 
zenith  angle,  all  azimuths.  The  solar  almucantar  is  the  sky 
at  the  zenith  angle  of  the  sun  fl, .  The  sky  radiance  in  the 
solar  almucantar  is  the  path  radiance  out-of-the- 
atmosphere  to  the  sensor  at  altitude  z .  The  sky  radiance 
can  be  expressed  as 


The  path  function  L«(z,0,,4)  can  be  expressed  In 
terms  of  a  sun  component,  a  sky-terrain  component  and 
an  emitted  component 

L»(z,fl,,<£)  -  ,«(z)a(z,/3) 

+  £  Z.(*,A',*’)<r(i,/9')4ft  +  L,fl(z)  ,  (3.2) 


L'lx,8s,  4)  -  ,*  Cz ) (cr  (z,/3) 

+  jf  Z.(z.0',*')<r(z.j8’)rfn/,«(z) 

+  £.»  (*)/,«(*))  ,  (3.3) 

where 

,c  (z)  is  the  solar  scalar  irradlance  at  altitude  z. 

Let  the  last  two  terms  be  designated  by  the  function  C(z) 
such  that 

C(z)-Jf  Mz,fl\*')(r(z,0’)dn/J«<z)  +  L*fl(z )/,«(!).  (3.4) 
Substituting  Eqs.  (3.3)  and  (3.4)  into  Eq.  (3.1)  wo  get 
Lt.  (*,0,,d>)  -  J  Ja(z()7,f<(*,9J)[a>(z/,/J)  +  CIZ/)]dr .  (3,5) 


The  sun  scalar  irradiance  and  the  transmittance  can 
be  expressed  as  a  function  of  the  sun  irradiance  at  the  top 
of  the  atmosphere ,« (<»)  as  follows 


,«(z()rr(  (z,fl, )  -  ,« (co)exp(-  J  « (z  )dr-  J"  o(zWr)  .  (3.6) 


The  two  integrals  can  be  combined  and  expressed  in  terms 
of  the  limits  «  to  z  which  is  equivalent  to  the  total 
transmittance,  thus 

,«(z, )r,((z,«,)-  ,€(»)r«,(z,0J)  -  ,«(:) ,  (3.7) 

Now  substituting  Into  Eq.  (3.5),  the  sun  irradiance  can  be 
taken  out  of  the  integral  and  we  get 

Lm  (z  .A, .  ^l/j*  (:)  “  J"  l(ir(z  ,fi)  +  C'iz  )ldr  .  (3.8) 


iy>>vyvy.  w  yy  vvwv.vw.yv 


The  quantity  in  the  square  brackets  is  a  function  of 
altitude  and  not  of  zenith  angle,  therefore 

J*  [tr Cr J3)  +  CU)\dr 

am 

“[mm{z,9s)/rnm(j, 0)]  +  C(r)]<*  .  (3.9) 


where  mw(z,0,)/m..(z,O)  is  the  relative  optical  airmass  at 
9,.  For  altitudes  up  to  6  kilometers  we  can  use  the  sea 
level  value  for  sun  zenith  angles  0  to  86  degrees  so  for 
convenience  we  will  shorten  this  to  m(6,).  Substituting 
back  into  Eq.  (3.8)  we  get 

£(*,«,,<*)/[,«(*  )«•(<>,))  -  j*  ((y(*,0)+  Ch))dz  .  (3,10) 

m 

The  right  hand  term  can  be  separated  into  two  parts 
f  [or(z,0)  +  C(x)]dz  -  J*  <r(z,0)A  +  J*  C(j)dz  .  (3.11) 

as  mm 


Now  let  us  define  an  optical  thickness  function  rM(z,0) 
such  that 


r„(»,0)  -  J*  <r(g,fi)dx  .  (3.12) 

m 


The  integral  of  the  optical  thickness  function  over  4  tr  is 
the  optical  thickness  due  to  scattering ,  tm  (r ) 

,/«.( i)- f  T„(*,0)rfn  (3.13) 


since  the  optical  thickness  due  to  scattering  is  the  Integral 
of  the  total  scattering  coefficient  with  altitude 


,'«(*)"  J*  i Car)rfr  .  (3.14) 

•• 

Substituting  Eq.  (3.12)  into  Eq.  (3.10)  we  get 


LL(:,9t,4>)/\st(j)m(9t))  -  3)  +  J"  C(z)A  .  (3.15) 


This  Is  the  basic  equation  for  the  radiance  in  the 
solar  aimucantar  and  how  It  relates  to  the  sun  scalar  irradi* 
ance,  the  optical  airmass,  the  optical  thickness  function 
and  one  additive  component  defined  by  Eq  (3.4). 


Optical  Stability 

Transmittancea  from  Solar  or  Solar  Aimucantar  Radiances 

The  atmosphere  is  considered  to  be  optically  stable 
when  the  atmospheric  transmittance  does  not  change  with 
time.  When  the  day  is  optically  stable,  measurements  of 
the  apparent  sun  radiance  ,I„(z,0, ,0)  made  over  a  large 
range  of  sun  zenith  angle  can  yield  a  good  measure  of 
total  transmittance  f.(*,  0)  and  the  inherent  sun  radiance 

,  Lm (i ,  0)  - , L0  (-iC^  (r,0) .  (3.16) 

Taking  the  log  of  both  sides  of  Eq.  (3.16)  we  get 

log,£„.Cr,0,,O)-  log ,Lt(oo)  +  m(6t)  log  r„(i,0) .  (3.17) 

Graphs  of  apparent  sun  radiance  on  semi-log  paper  of  log 
apparent  sun  radiance  versus  airmass  yield  the  typical 
Langley  (or  Bouguer)  graph  where  log  inherent  sun  radi¬ 
ance  is  the  Intercept  and  log  transmittance  is  the  slope. 

Variations  of  Eq.  (3.15)  give  both  another  method 
for  obtaining  the  total  transmittance  and  several  ways  of 
checking  the  optical  stability,  Equation  (3.15)  can  also  be 
written  as 


C  (r,0],/})/m(9J) 

-  r.^.O)"®'’ (1,0)  +  /  C(*)<6] .  (3.18) 

m 

Taking  the  log  of  both  sides  we  get 

log(L»  (z,0,,0)/«(9,) I  -  mtoj)  logr.U.O)  +  log/l  (3.19) 
where 

A  -  ,«(oo)[r„(j,0)  +  J*  C(x)dt] .  (3.20) 

m 

For  an  optically  stable  day,  A  should  be  constant  for 
a  given  angle  from  sun.  Thus  the  sky  radiance  at  a  con¬ 
stant  0  in  the  solar  aimucantar  can  be  used  to  obtain  the 
beam  transmittance  when  measured  over  a  large  range  of 
airmass  values  on  an  optically  stable  day.  A  semi-log 
graph  of  !og(C  Cr,0,,0)/«(0,)1  versus  airmass  would  be 
linear  with  a  slope  of  log  transmittance  and  an  intercept  at 
log/4 .  The  transmittance  obtained  from  the  sun  radiance 
using  Eq,  (3.17)  and  the  transmittance  obtained  from  the 
sky  radiance  in  the  solar  aimucantar  are  then  averaged  to 
obtain  the  most  accurate  value. 

Tashenov  (1970)  used  this  method  with  measure¬ 
ments  of  the  solar  aureole  to  obtain  spectral  transmittance 
in  the  region  410  to  735  nanometers. 


Tests  to  Determine  Optical  Stability 

Pyaskovskaya-Fesenkova  (1970)  outlines  two  tests 
for  optical  stability  as  well  as  a  unique  equation  for 
transmittance  which  is  valid  only  for  optically  stable  days. 

The  first  test  is  to  graph  for  a  constant  scattering 
angle  the  ratio  L'm  Or,0,  ,/3)/(  ,«(z)m(0,)]  versus  w(0,), 
As  can  be  seen  from  Eq.  (3.15),  if  the  day  is  optically 
stable  the  right  hand  term  is  constant.  Therefore  this 
graph  should  result  in  a  horizontal  straight  line.  That  is, 
for  an  optically  stable  day 

(3.21) 

where 

B  -  r„(t ,0)  +  j  C(i)dz  .  (3.22) 


The  second  test  for  optical  stability  is  to  graph  the 
ratio  Ll  (z,0,,/3)/,e (z)  versus  relative  airmass  for  a  con¬ 
stant  13,  The  graph  should  result  In  a  straight  line  going 
through  the  origin  on  an  optically  stable  day.  To  under¬ 
stand  this,  Eq.  (3.15)  is  rewritten  in  the  form 

Li(i>9s,0)/s*U)-m(6a)B .  (3.23) 

The  B  is  now  the  slope  and  the  Intercept  is  zero. 

If  the  day  is  optically  stable  a  graph  of  solar  almu- 
cantar  radiance  at  a  constant  scattering  angle  versus  air- 
mass  will  yield  a  value  of  transmittance  as  follows.  If  we 
write  Eq.  (3.21)  in  the  form 

C  (z  . 0, ,  p )  - ,« (oo ) Tm  (i ,  0)" ' *» * '  m  (0, )  B  (3.24) 

and  differentiate  (only  the  sky  radiance  and  the  air  mass 
are  variables,  the  inherent  sun  scalar  Irradlance,  B  and  the 
transmittance  are  constant),  we  get 

dil  (2,9,  ,0) 

-  B  s e (oo  1 7*0. (z , 0)m W J *  d  m(0j)(l  +  m(9s)lnTm(i  ,0)]  ,  (3.25) 

Equation  (3.25)  Is  only  valid  when  the  angle  from  sun  is 
constant,  the  vertical  transmittance  from  space  to  sensor  Is 
constant  and  the  day  is  optically  stable,  Setting 
dLl  ( 2 ,0,  ,0)  equal  to  zero  and  rearranging  we  get 

In  T„  (2  ,0)  -  -  1/m  (flj )  (3  26) 

where  6,  is  the  sun  angle  where  L'm  (:,9,fi)  is  a  max¬ 
imum.  Therefore  a  graph  of  solar  almucantar  radiance  at 
a  constant  angle  from  the  sun  versus  airmass  will  increase 


and  then  decrease  with  the  maximum  at  m(0,)  when  the 
day  is  optically  stable.  If  the  curve  is  reasonably  well 
defined  with  sufficient  airmass  values,  a  reasonable 
transmittance  value  can  be  obtained.  It  can  also  be  used 
as  a  third  test  for  optical  stability.  The  values  of  transmit¬ 
tance  obtained  by  means  of  Eqs.  (3.17),  (3.19)  and  (3.26) 
should  be  in  agreement  on  an  optically  stable  day. 

Pyaskovskaya-Fesenkova  (1970)  describes  an  instru¬ 
ment  for  measuring  the  quantities  in  the  preceding  equa¬ 
tions.  It  Is  essentially  a  radiance  photometer  with  a  three 
degree  field  of  view  for  the  solar  measurement.  For  the 
aureole  measurement  a  two  degree  to  three  degree  circular 
zone  is  used.  Her  article  also  includes  a  theoretical  test  of 
the  error  in  obtaining  transmittance  from  the  sun  radiance 
when  the  aureole  in  the  three  degree  field  is  necessarily  a 
part  of  the  sun  radiance  measurement.  She  concluded  the 
error  to  be  negligible. 


from  Horizon  Sky 

Roessler  and  Faxvog  (1981)  developed  equations 
for  the  horizon  sky  luminance  as  a  function  of  the  single 
scattering  albedo  after  assuming  both  an  isotropic  scatter¬ 
ing  function  and  an  isotropic  luminance  distribution.  The 
technique  can  be  made  more  generally  applicable  by 
developing  the  equations  for  the  horizon  at  a  scattering 
angle  of  55  derrees  .rom  the  sun. 


The  equation  for  the  path  function  at  55  degrees 
scattering  angle  can  be  expressed  us  a  function  of  the  sun 
component  separate  frum  the  sky-terrain  component  in  a 
manner  similar  to  Eq.  (3.2) 

/.*(z,0,  55/9)  -  4u(r )tr (r , 55) 

+  f  L(z%9‘,d>')<T(z,l}')d{\  +  L,a(i) .  (4.7) 


Substituting  Eq,  (4.6)  Into  Eq.  (4.7)  and  assuming  the 
sky-terrain  or  diffuse  scalar  irradiance  j*  (z )  Is  isotropic  we 
get 

L,(2,9,  55 1})  -,«(.-)*  <2)/(4tt) 

+  rf«(z)j(r)/(4tr)  +  Lt  (z) .  (4.8) 

a 

Since  the  total  scalar  irradiance  is  the  sun  plus  the  sky- 
terrain  scalar  irradiance,  Eq.  (4.8)  can  be  written  as 

L, (2,9,  55/3)  -«(z)j (r)/(4rr)  4  L,a(:)  .  (4.9) 


ADDITIONAL  DETERMINATIONS 
(Section  4,  AFOL-TR-82-0223) 

Method  of  Obtaining  Single  Scattering  Albedo 


Horizon  Sky  at  55°  from  Sun 
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Dividing  both  sides  by  the  attenuation  coefficient  and 
expressing  the  emitted  path  function  according  to 
Eq.  (2.7)  we  get  an  equation  for  the  equilibrium  radiance 

Lq  (r,9,  55  0)  —  «(r)r(r)/ta(r)4ir] 

+  fl(z)t(t, D/o(z) .  (4.10) 

This  can  be  expressed  in  terms  of  the  single  scattering 
ulbedo  to  (from  Eq.  (2.39))  as 

Lq(j,8,  550)  -  oi  (r  )«(*)/ (4  ir) 

+  (l-«(z)U.0t,D.  (4.11) 

The  general  equation  for  the  horizon  sky  radiance  at 
55  degrees  scattering  angle  from  Eq.  (2.S7)  is 

/I (a 90 550)-  L,U9<t550)[l-  r„ 0(90)1 .  (4.12) 

Therefore  from  Eq.  (4.11)  the  horizon  sky  radiance  at  55 
degrees  is 

Ll  (r, 90. 550) 

-U(r)«(z)/<4ir)+  ll-w<z)]i,(X,r))tl-  r„(r,90)] .  (4.13) 

The  above  equation  assumes  that  the  attenuation 
coefficient  is  horizontally  isotropic  and  therefore  the  hor¬ 
izon  sky  must  be  cloud  free. 


Visible  Spectrum 

In  the  visible  spectrum  L(X,T)  is  negligible.  Also 
for  the  photopic  at  sea  level  T*(. 0,90)  is  negligible  and 
Eq  (4.13)  becomes 


t„  (0,90,550)  -  w(0)« (0)/4ir  .  (4.14) 

Thus  the  single  scattering  albedo  u  can  be  determined  by 
measuring  the  horizon  sky  at  55  degrees  and  the  total 
scalar  irradiance  at  sea  level 

<u(0)  -  Ll  (0,90, 550)4  ir/«(0) .  (4.15) 

Since  this  is  an  approximative  method,  the  precision  level 
of  the  resultant  <u  (0)  is  not  expected  to  be  high.  How¬ 
ever,  it  can  still  be  quite  useful  as  one  measure  of  the  sin¬ 
gle  scattering  albedo. 

When  there  is  no  absorption  w  -  1  and  the  horizon 
sky  should  equal  the  average  scalar  irradiance  e  (0)/  (4ir) 

L*  (0,90, 550) -«(0)/(4n)  .  (416) 


Method  of  Obtainina  Scattering  Transmittance 
from  Sky  Radiance  Ratios  at  55  Pearces 

When  absorption  is  negligible,  the  general  equation 
for  the  equilibrium  radiance  at  55  degrees  scattering  angle 
Eq.  (4.11)  reduces  to 


£,(*. 8  $50)  -  « (r)/4ir . 


Also  when  absorptio.  .  negligible  the  total  scalar  irradi¬ 
ance  c  (z )  is  constant  with  altitude,  therefore  the  equili¬ 
brium  radiance  at  55  degrees  is  also  constant  with  altitude. 
Rearranging  the  equilibrium  radiance  form  of  the  equation 
of  transfer  Eq.  (2.17)  we  get 

• 

[L(r,fl,  55 0)-L?(*,fl,  550)1-'  *(*,«, 55  0) 

- -s(z)dr ,  (4.18) 

where  the  total  scattering  coefficient  s(z)  has  been  substi¬ 

tuted  for  the  attenuation  coefficient  since  absorption  is 
negligible.  Since  the  equilibrium  radiance  is  a  constant 
Eq.  (4.18)  can  be  integrated  with  respect  to  altitude 

* 

t  [L(z.S.  550)  -£,(*,«,  550)]-' 0L(r, 9,550)  . 

--/i(r)<fr.  (4.19) 

0 

The  result  of  the  integration  is 

L,  fa  ,fl,  55  0)  -  L0  Or, ,  8,  55  0)  f  Tr  (r  .•> 

+  L,(M,  55 0)tl-,rf(r, 0)1  (4.20) 

where ,  Tr  (z  ,9)  is  the  transmittance  due  to  scattering. 

Sky  Radiance  Ratio  at  55s 

For  sky  radiance  at  55  degrees  scattering  angle,  the 
inherent  radiance  £,  is  equal  to  zero  and  Eq.  (4.20) 
becomes 

Ll  (z,0,  550)  —  Lq  (*,#,  550)11-  .  (4.21) 

Now  when  55  degree  sky  radiances  at  two  zenith 
angles  6  and  9'  are  ratioed,  since  the  equilibrium  radiance 
is  equivalent,  the  sky  radiance  ratio  is  solely  a  function  of 
the  scattering  transmittance 

Ll(z,9,  550)/ Lm  U,8\  550) 

-  [i- Jr„u,o)m»>]/  li- ,r.. u,o)'"»'>] .  (4.22) 


**«  **v**Or-  /o'  •  "pv 


Although  Eq.  (4,22)  cannot  be  solved  for  directly,  it  can 
be  solved  by  iterative  means.  Error  analysis  indicates  that 
the  zenith  angle  difference  9-9'  should  be  large  to 
minimize  the  error  in  the  resultant  vertical  transmittance 


Visible  Spectrum 

In  the  visible  spectrum,  absorption  is  negligible 
except  for  ozone.  The  total  transmittance  in  the  visible 
spectrum  would  thus  be  approximated  by  the  product  of 
the  scattering  transmittance  times  the  ozone  transmittance 

rMOt.O)-  eTm(t,0) .  (4.23) 

To  illustrate  the  use  of  Eqs.  (4.22)  and  (4.23),  we 
again  refer  to  the  sky  radiance  and  solar  transmissometer 
data  taken  on  the  Visibility  Laboratory  rooftop  in  1964. 
The  atmosphere  on  2  September  1964  was  optically  stable 
during  the  afternoon  as  illustrated  in  Section  3.3.  Thus 
evaluation  of  the  scattering  transmittances  throughout  the 
afternoon  should  give  us  an  estimate  of  the  precision  of 
the  estimation  method  for  broad  band  sensors  in  the  visi¬ 
ble.  Equation  (4.22)  was  evaluated  for  all  sky  radiances  at 
55  degrees  plus  or  minus  2.5  degrees  scattering  angle  for  9 
from  81.6  to  64.7  degrees  and  9'  from  64.7  to  2.8  degrees. 
Later  error  analyses  indicated  that  some  of  these  zenith 
angle  combinations  are  less  error  prone  than  others. 
Hence,  the  averages  in  Table  4.1  are  less  accurate  than  can 
be  obtained  with  a  smaller  yet  better  selection  of  zenith 
angle  combinations.  However,  they  are  presented  herein 
u  a  first  approximation. 

The  transmittance  ratios  in  Table  2.1  are  the 
transmittances  based  on  Eq.  (4.22)  times  the  ozone 
transinittance  divides!  by  the  transmittance  from  the  solar 
transmissometer.  The  ozone  transmittance  is  noted  at  the 
bottom  of  each  column.  The  transmittance  comparison  is 
best  for  Filter  1  which  is  a  relatively  narrow  band  filter 
with  the  least  absorption.  The  near  noon  data  compare 
least  well.  The  accuracy  is  probably  sufficient  to  warrant 
further  development  of  the  method  to  improve  the  preci¬ 
sion  level. 

Use  of  sky  radiances  at  55  degrees  scattering  angle 
should  be  an  improvement  on  the  sky  ratio  method 
described  in  Duntley  eta!.  (1972)  Section  2.1  and  in 
Duntley  etaL  (1978)  Section  2.2.  That  method,  which 
stemmed  from  the  nomographic  method  of  Kushpil'  and 
Petrova  (1971),  used  sky  radiances  at  all  scattering  angles. 
Kushpil*  and  Petrova  (1971)  suggested  use  of  ratios  at 
57.2  degrees  scattering  angle  for  the  visible  spectrum  and 
53.9  degrees  for  the  near  infrared  portion  of  the  spectrum 
but  did  not  give  a  theoretical  basis  for  the  nomograph  or 
for  the  selection  of  these  angles. 

SUMMARY 

(Section  5,  AFGL-TR-82-0223) 

Some  implications  of  the  equation  of  transfer  as  it 


Table  2.1.  Comparison  of  transmittance  determination  pro¬ 
cedures.  (sky  radiance  ratio  method  w  solar  transmissometer). 
(Table  4.1  from  AFGL-TR-82-0223) 
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relates  to  a  scattering  and  absorbing  medium  have  been 
explored.  The  mgjor  implication  is  that  a  measurement  of 
the  4ir  radiance  distribution  can  yield  a  great  deal  of  infor¬ 
mation  about  the  atmosphere  both  in  the  visible  and  the 
near  infrared.  If  the  solar  almucantar  is  cloud-free  the 
aerosol  optical  thickness  can  be  derived.  If  the  4w  radi¬ 
ance  distribution  is  measured  at  several  altitudes,  the  radi¬ 
ance  arrays  can  be  tested  for  consistency  and,  if  con¬ 
sistent,  a  measure  of  absorption  obtained.  The  Visibility 
Laboratory  data  catalog  of  radiances,  both  airborne  and 
ground-based  are  a  mine  of  information  waiting  to  be 
tapped.  This  catalog  is  useful  for  testing  methods  of 
retrieval  of  basic  scattering  and  absorption  information 
about  the  atmosphere  as  well  as  for  development  of 
atmospheric  models. 

2.2  AFGL-TR-83-0236:  Analytic  Techniques  for 
Estimating  Visible  Image  Transmission  Properties 
of  the  Atmosphere 

by  Wayne  S.  Hering 

Introductory  Comments:  A  fast  operationally  oriented 
technique  for  estimating  visible  spectrum  contrast 
transmittance  was  introduced  in  Hering  (1981).  This  1981 
report,  recognized  that  the  complexity  of  multiple  scatter¬ 
ing  and  absorption  processes,  in  conjunction  with  the 
Inherent  variability  of  the  aerosol  particle  distribution  and 
composition  would  require  effective  simplifications  in 
order  that  estimates  of  contrast  transmittance  through  the 
atmosphere  be  made  rapidly  and  consistently. 


Thus,  a  series  of  modelling  approximations  that 
relate  optical  scattering  properties  to  meteorological  vari¬ 
ables  were  derived  from  experimental  data,  and  combined 
with  available  analytic  approximations  for  radiative 
transfer  calculations  to  develop  the  operational  technique 
discussed  in  that  report. 

A  subsequent  report,  AFGL-TR-83-0236,  Hering 
(1983),  was  prepared  during  this  current  contract  interval 
which  extends  the  methodology  of  Hering  (1981)  to 
include  not  only  clear  and  hazy  day  meteorological  condi¬ 
tions,  but  full  overcast  cloud  conditions  as  well.  Since  the 
work  reported  in  Hering  (1983)  lays  the  fundamental 
ground  work  for  the  FASCAT  model  discussed  in  Sec¬ 
tion  3  of  this  Final  Report,  the  text  of  the  1983  document 
has  been  included  herein  as  Appendix  A,  and  need  not  be 
discussed  further  here. 


2.3  AFGL-TR-83-0270:  Characteristics  or 
Aerosol  Volume  Distributions  Measured  at 
Meppen,  W.  Germany 
by  Bruce  W.  Fitch 

Introductory  Comment:  One  of  the  most  vexing  prob¬ 
lems  facing  those  whose  goal  is  the  full  characterization  of 
the  atmosphere's  influence  upon  optical  propagation,  is 
the  reliable  specification  of  the  aerosol  size  distribution. 
This  report,  AFGL-TR-83-0270  addresses  this  problem  via 
the  analysis  of  a  composite  data  base  collected  by  two 
different  measuring  systems,  both  of  whose  pertinent 
characteristics  are  described.  The  analysis  uses  the 


volume  distribution  methodology  discussed  in  several  ear¬ 
lier  publications,  Whitby  et  at.  (1972),  Shettle  (1975).  In 
the  discussion,  three  size  modes  are  defined  as  the  accu¬ 
mulation  mode  (sub-micron  region),  the  course  mode 
(greater  than  ip  region)  and  the  middle  mode  (occurring 
at  particle  radii  of  approximately  one  micron).  In  general, 
log  normal  fits  to  each  of  the  size  distributions  are  used 
throughout. 


A  CASE  STUDY 

(Section  5,  AFGL-Tk-83-0270) 

From  4  March  through  13  March,  1978,  Royco 
measurements  were  made  several  times  each  day  around 
local  noon.  Figure  2-1  shows  the  values  Af|,  Mi,  and 
along  with  wind  speed,  wind  direction  and  relative  humi¬ 
dity  for  these  ten  days.  Each  column  in  the  figure 
represents  a  full  24  hr  day.  At  the  start  of  the  series  on 
4  March,  a  large  middle  mode  is  present  but  there  is  no 
accumulation  mode.  Note  that  for  this  day,  the  values  of 
are  plotted  on  a  X10  scale.  A  large  coarse  mode 
appears  in  two  samples  measured  near  1300  hours.  From 
0$00  hours  up  through  1430  hours,  the  time  of  the  last 
sample,  the  relative  humidity  is  >85%  and  the  wind  is 
from  the  north  at  a  speed  of  <  2  m/sec.  The  visible  (pho- 
topic)  extinction  coefficient  is  2.9  km-1  at  1100  hours  and 
1.9  km-1  at  1400  hours  indicating  the  presence  of  a  heavy 
mist  or  light  fog.  A  mist  is  considered  present  when  visi¬ 
ble  extinction  is  >  1  km'1  and  relative  humidity  is  >80%; 
a  fog  is  present  when  visible  extinction  is  — 4  km'1.  Fig. 
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FI*.  2-1.  Maximum  volume  concentrations  M  of  the  three  modes  with  relative  humidity  (RH),  wind  direction,  wind 
speed  and  periods  of  rain  (R)  for  ten  days  in  March.  (Fig.  5-1  from  AFGL-TR-83-0270). 


2-2 (a)  shows  the  volume  distribution  measured  four 
different  times  on  4  March.  The  figure  shows  a  decrease 
in  the  volume  Mi  over  the  time  interval  1005  through 
1428  hrs  corresponding  to  a  decrease  in  the  extinction 
coefficient.  The  mist  persists  throughout  4  March  and  is 
gone  by  OSOO  hours  on  5  March. 

The  middle  mode  is  gone  by  1400  hours  on 
5  March.  High  relative  humidity  and  visible  extinction 
through  the  night  of  4  March  and  early  morning  of 

5  March  suggest  the  presence  of  a  middle  mode  till  about 
1000  hours.  At  this  time  there  is  a  sharp  increase  in  wind 
speed  followed  by  a  decrease  in  relative  humidity.  This 
suggests  that  the  air  parcel  was  dried  through  either  verti¬ 
cal  mixing  or  advection.  It  is  likely  that  the  drop  in  rela¬ 
tive  humidity  to  values  <65%  dried  the  middle  mode. 
Prior  to  the  aerosol  measurements  on  5  March  there  were 
three  occurrences  of  rain.  It  is  doubtful  that  the  rain 
washed  out  the  middle  mode  because  measurements  on 
several  other  days  show  a  middle  mode  present  immedi¬ 
ately  following  an  episode  of  rain.  The  accumulation  and 
coarse  particle  modes  remain  fairly  constant  from  4  to 

6  March. 

There  is  a  marked  increase  in  the  values  of  M\ ,  M: 
and  M)  on  7  March  as  compared  to  their  values  on  5  and 

6  March.  Corresponding  to  the  change  in  the  magnitude 
of  the  three  maximum  volumes,  there  is  a  shift  in  wind 
direction  from  300-350°  to  200-250°.  The  major  source  of 
aerosol  in  the  300-350°  direction  is  the  North  Sea; 
whereas  n  the  200-250°  direction  there  are  several  aerosol 
sources  including  the  industrial  center  of  Dusseldorf  at 
about  200°  and  Amsterdam  and  the  North  Sea  near  250°. 
This  suggests  that  the  accumulation  and  coarse  particle 
modes  are  associated  with  an  air  mass  from  an  industrial 
center.  The  middle  mode,  whose  presence  is  dependent 
on  high  relative  humidity,  is  probably  a  marine  aerosol  or, 
possibly,  a  very  hydroscopic  component  of  the  accumula¬ 
tion  mode.  The  three  volume  modes  measured  three 
times  on  7  March  are  shown  in  Fig.  2-2(b).  The  relative 
humidity  is  >75%  for  the  aerosol  measurements  on 

7  March  and  the  visible  extinction  is  about  1  knr1 ,  Dur¬ 
ing  the  time  of  aerosol  measurement,  the  wind  speed  is 
about  5m/sec,  suggesting  an  advection  mist.  Conditions 
on  8  March  are  similar  to  those  on  7  March  except  7j 
increases  to  ~85 ^m}cm~3.  The  composition  of  the 
coarse  mode  is  not  known  but  the  modest  relative  humi¬ 
dity  present  at  the  time  of  measurement  suggests  the  par¬ 
ticles  are  fairly  dry.  Late  on  8  March  the  wind  direction 
shifts  from  around  200°  to  330°. 

The  first  measurements  of  9  March  show  the  pres¬ 
ence  of  the  middle  mode  but  not  the  accumulation  and 
coarse  particle  modes.  The  rapid  decrease  in  Mi  after 
1100  hours  is  in  concert  with  a  decrease  in  relative  humi¬ 
dity.  Fig.  2-2(c)  shows  the  decrease  in  the  size  of  the 
middle  mode  from  1015  to  1655  hours.  The  sharp 
decrease  in  relative  humidity  ir.  associated  with  an  increase 
in  wind  speed  similar  to  what  happened  on  5  March  at 
about  the  same  time  of  day.  In  the  early  evening  on 
9  March  the  wind  direction  shifts  to  200-250°  and  the  rela¬ 
tive  humidity  increses  sharply.  With  this  change,  Fig.  2- 


2(c)  shows  a  significant  increase  in  the  large  particles  at 
1655  hours.  The  following  day  a|l  three  volume  modes 
are  again  present^ 

The  appearance  of  the  three  modes  on  10  Msrch 
with  the  wind  from  200-250°  and  relative  humidity  >  80% 
is  similar  to  the  situation  on  7  March.  Fig.  2-2(d)  shows 
the  three  volume  modes  measured  on  10  March.  A  meas¬ 
ure  of  the  visible  extinction  on  10  March  is  not  available 
but  the  value  of  the  scattering  coefficient,  2.5  km-1,  indi¬ 
cates  a  mist.  Since  absorption  is  small  at  visible 
wavelengths,  the  scattering  and  extinction  coefficient  will 
be  used  interchangeably  to  indicate  a  mist  or  fog.  An 
advection  mist  is  indicated  since  wind  speed  is  >  3  m/sec. 
During  10  March,  the  wind  direction  slowly  becomes  more 
northerly  indicating  that  the  next  day  may  show  a  middle 
mode  if  relative  humidity  is  high. 

On  1 1  March  the  relative  humidity  is  >  85%  and  the 
wind  direction  is  fluctuating  from  about  250  to  330*.  The 
visible  scattering  coefficient  decreases  rapidly  from  a  value 
of  —10  kmr1  at  0500  hours  to  about  1  km-1  at  1100 
hours.  By  1300  hours  it  is  about  0.6  km"1.  The  values  of 
M  indicate  that  all  three  modes  are  present  but  in  small 
number  concentration.  A  small  number  concentration 
would  explain  why  the  middle  mode  is  not  large  with  rela¬ 
tive  humidity  persisting  above  85%.  Thus,  it  appears  that 
an  advection  mist  was  present  in  the  morning  but  by  mid¬ 
day  was  replaced  with  an  air  mass  high  in  humidity  and 
low  in  particle  concentration. 

On  12  March,  the  wind  is  from  100-160°,  relative 
humidity  is  >75%  and  all  three  volume  modes  are 
present.  The  wind  direction  and  large  accumulation  mode 
indicate  the  influence  of  an  industrial/urban  aerosol 
source.  Late  in  the  day  the  wind  direction  changes  to 
200-250°  suggesting  all  three  modes  will  again  be  present 
on  13  March. 

Fig.  2-2(e)  shows  that  all  three  modes  are  present 
on  13  March.  On  this  day,  the  values  of  M,  for  the  accu¬ 
mulation  mode  are  smaller  compared  to  previous  cases 
where  this  mode  was  present  for  wind  from  200-250°. 
Clean  air  generally  associated  with  the  passage  of  a  cold 
front  In  conjunction  with  a  small  production  of  industrial 
particulates  on  Sunday,  12  March,  may  be  responsible  for 
M\<nm1cm~i  on  13  March.  A  frontal  passage  is  indi¬ 
cated  by  cooler  daytime  temperatures  on  13  March  com¬ 
pared  to  12  March  and  a  long  period  of  rain  during  the 
night  of  12  March.  There  is  a  sharp  increase  i  t  wind 
speed  at  about  1000  hrs  on  13  March  followed  by  a 
decrease  in  relative  humidity  and  M}.  This  same  meteoro¬ 
logical  phenomenon  occurred  on  5  and  7  March  and  simi¬ 
larly  resulted  in  a  decrease  in  Mi. 

The  ten  day  case  study  presented  in  Fig.  2-1  shows 
that  a  knowledge  of  local  wind  direction,  wind  speed  and 
relative  humidity  is  not  always  sufficient  to  predict  the 
characteristics  of  the  local  aerosol  size  distribution.  High 
relative  humidity  is  important  to  aerosol  growth  but  the 
measured  volume  modes  are  not  equally  affected.  Parti¬ 
cles  of  the  middle  mode  are  much  more  sensitive  to 
changes  in  relative  humidity  than  those  of  the  other  two 
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}  volume  modes.  The  effect  of  an  increase  in  relative  humi¬ 
dity  on  extinction,  however,  not  only  depends  on  which 
modes  are  present  but,  also,  on  the  initial  number  concen¬ 
tration.  At  Meppen,  particles  comprising  the  three  modes 
are  advected  in  from  remote  sources.  Thus,  an  estimate 
of  the  values  M  cannot  be  based  solely  on  a  measure  of 
the  local  wind  direction  but  must,  also,  include  an  esti¬ 
mate  oT  the  type  and  strength  of  the  downwind  aerosol 
sources.  Local  sources  of  aerosol  associated  with  farm 
operations,  such  as  plowing  and  burning,  may  become 
important  additional  contributors  during  the  summer 
months.  The  case  study  shows  that  a  measure  of  local 
meteorology  is  useful  for  estimating  the  presence  or 
absence  of  the  three  modes  and  for  predicting  changes  in 
the  values  of  M  but  is  not  useful  for  predicting  the  magni¬ 
tudes  of  M. 
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SUMMARY  A  CONCLUSIONS 
(Section  6.  AFOL-TR-83-0270) 

This  report  presents  the  results  of  analysis  of  a  data 
base  consisting  of  aerosol  size  distributions  measured  con¬ 
currently  with  the  local  meteorology  at  Meppen.  The  pur¬ 
pose  of  the  analysis  was  to  corroborate  an  extensive  set  of 
airborne  aerosol  measurements  described  in  Fitch  and 
Cress  (1981  and  1983)  and  to  determine  the  utility  of  local 
meteorological  measurements  in  uniquely  defining  the 
local  aerosol  size  and  number  distribution. 

The  characteristics  and  behavior  of  the  ground 
measured  aerosol  volume  distributions  agree  well  with  the 
airborne  measured  distributions.  Volume  distribution 
plots  for  both  data  sets  show  three  volume  modes,  each 
well  fit  by  a  log-normal  distribution.  The  behavior  and 
values  of  the  log-normal  fit  parameters  for  the  accumula¬ 
tion  and  middle  modes  are  in  good  agreement  as  shown  in 
Table  2.2.  The  table  presents  the  average  value  of  stan¬ 
dard  deviation  <r  and  the  range  of  mode  radii  7  for  each  of 


Table  2.2.  Summary  of  the  volume  mode  parameters, 
(Table  6.1  from  AFGL-TR-83-0270). 
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the  three  modes  that  appear  in  both  the  ground  and  air¬ 
borne  data  bases.  The  greater  frequency  of  occurrence 
and,  generally,  bigger  size  of  the  ground  coarse  modes 
compared  to  the  airborne  coarse  modes  indicates  a  larger 
average  particle  number  concentration  for  the  ground 
modes.  The  lack  of  agreement  between  or3  for  the  two 
data  sets  is  probably  a  result  of  low  particle  counts  in  the 
airborne  data  and  an  incomplete  sampling  of  the  entire 
coarse  mode.  A  fit  of  three  log-normal  curves  to  the 
measured  distributions  of  both  data  sets  explains  more 
than  90%  of  the  variance.  The  recurring  behavior  and, 
often,  constant  value  of  the  fit  parameters  7  and  <r  for  the 
three  modes  show  that  the  aerosol  distribution  over  the 
measured  particle  radius  range  could  be  estimated  using 
only  a  measure  of  the  maximum  volumes  M\,  A/2,  and 
A/j. 

Concurrent  measurements  of  the  aerosol  distribu¬ 
tion  and  local  meteorology  show  that  the  formation  of 
large  particle  volumes  associated  with  mist  and  fog 
requires  high  values  of  relative  humidity.  High  relative 
humidity,  however,  does  not  always  produce  a  large  parti¬ 
cle  volume.  A  case  study  of  the  aerosol  distribution  and 
meteorology  indicates  conditions  which  promote  high 
humidity  do  not  always  insure  the  presence  of  an  aerosol 
necessary  for  the  formation  of  a  mist  or  fog.  The  PMS 
data  show  that  the  local  meteorology  is  not  a  definitive 
predictor  of  aerosol  presence.  During  a  mist,  the  Royco 
data  show  that  the  aerosol  distribution  can  be  a  middle 
mode  solo  or  all  three  modes  in  combination.  It  is  antici¬ 
pated  that  as  more  aerosol  and  meteorological  data 
become  available,  the  task  of  determining  the  presence 
and  characteristics  of  the  aerosol  distribution  based  on  a 
measure  of  the  meteorology  will  be  continued.  It  is 
recommended  that  future  studies  concentrate  on  the 
behavior  of  the  three  volume  modes. 


The  data  analysis  indicates  that  accurate  short  term 
prediction  of  aerosol  size  distribution  morphisis  is  possible 
by  combining  methodologies  of  predicting  local  episodes  of 
mixing  and  changes  in  relative  humidity  with  measure¬ 
ments  of  A/|,  A/j  and  A/3.  This  technique  would  avoid 
the  complexities  and  expense  involved  in  developing  an 
elaborate  scheme  which  uses  air  mass  trajectory  analyses 
in  conjunction  with  a  location  and  time  dependent  aerosol 
source  inventory.  With  a  local  measure  of  the  values  A/, 
there  would  be  no  need  to  develop  or  maintain  an  aerosol 
source  inventory  or  to  predict  changes  in  the  aerosol  dis¬ 
tribution  during  long  passages  from  its  source  to  the  site 
of  interest.  The  data  suggest  that  extinction  in  the  visible 
and  infrared  could  be  predicted  using  LOWTRAN,  a 
measure  of  the  present  values  of  M  and  a  six  to  twelve 
hour  estimate  of  future  meteorological  conditions.  Such  a 
methodology  would  not  require  a  complex  aerosol  source 
inventory  and  could  provide  a  mesoscale  map  of  the  opti¬ 
cal  environment  when  implemented  over  a  network  of 
meteorological  stations. 
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2.4  AFGL-TR-84-0096:  A  Traniportablc,  Machine 
Oriented  Library  of  European  Sky  and  Terrain 
Radiance  Distributions  with  Contemporary 
Radiometric  and  Meteorological  PnlUes 

by  Richard  W.  Johnson  k  Miriam  K.  Oleinik 

Introductory  Comment:  The  radiometric  and  meteoro¬ 
logical  data  base  upon  which  most  of  the  optical  properties 
modelling  discussed  in  previous  sections  of  this  Fmal 
Report  has  been  based,  was  obtained  by  the  C-130 
mounted  system  described  in  AFGL-TR-81-0275,  Johnson 
(1981).  These  data,  similar  in  nature  to  much  of  that 
obtained  since  1968  by  the  Visibility  Laboratory,  have  not 
until  now  been  readily  available  to  the  general  public.  A 
large  of  sample  of  these  radiometric  and  meteorological 
measurements  representing  nearly  500  4rr  radiance  distri¬ 
butions  at  5  degree  angular  resolution  are  now  in  a  tran¬ 
sportable  tape  format. 

For  those  readers'  not  familiar  with  the  nature  of 
these  radiometric  measurements,  two  illustrative  examples 
are  included  in  this  Final  Report,  AFGL-TR-84-0168,  for 
quick-look  orientation.  Data  flight  number  C-379  was 
conducted  south  of  Lolland  Is.,  Denmark  on  17  May  1976. 
It  was  a  mid-day  flight  spanning  local  apparent  noon  under 
clear  day  meteorological  conditions.  The  data  for  this 
flight  are  included  in  the  herein  described  transportable 
library  and  are  illustrated  photographically  in  Fig.  2-3,  and 
radiometrically  in  Fig.  2-4.  Both  of  these  figures  originally 
appeared  in  AFGL-TR-81-0317  an  early  report  analyzing 
the  overall  data  set, 

This  report,  AFGL-TR-84-0096,  describes  both  the 
data  base,  and  the  tape  format,  with  sample  extraction 
routines  included. 


DESCRIPTION  OF  DATA  BASE  CONTENT 
(Section  2,  AFGL-TR-84-0096) 

The  tape-oriented  data  base  that  is  discussed  in  this 
report  contains  two  independent  but  related  sets  of  meas¬ 
urements  which  are  provided  on  several  separate  storage 
tapes  whose  format  will  be  discussed  in  the  following  sec¬ 
tion.  In  this  section,  we  will  confine  ourselves  to  a  brief 
discussion  of  the  pertinent  characteristics  of  the  measured 
data  itself  rather  than  that  of  the  storage  medium. 

The  two  data  sets  comprising  this  data  base  are 
described  herein  as;  a)  the  Sky  and  Terrain  Radiance 
Data,  and  b)  the  Scattering  Coefficient  and  Related 
Meteorological  Data.  For  convenience,  the  Sky  and  Ter¬ 
rain  Radiance  Data  will  often  be  referred  to  as  the 
"Scanned  data,  and  the  Scattering  Coefficient  and  Related 
Meteorological  Data  as  the  "  Profile"  data.  As  a  general 
rule,  the  data  are  sorted  into  sub-sets  identified  by  flight 
number  as  indicated  in  Tables  2.3  and  2.4. 

As  with  all  of  the  radiometric  data  associated  with 
the  reports  listed  in  Table  2.5,  the  radiometric  measure¬ 
ments  contained  in  this  data  library  may  be  characterized 
by  the  illustrative  information  in  Fig.  2-5  and  Table  2.6. 
More  detailed  information  regarding  the  radiometric  cali¬ 
bration  procedures  applied  to  these  data  is  available  in 
several  preceaing  reports,  AFCRL-70-0137,  AFCRL-72- 
0461,  and  AFGL-TR-80-0207  among  others  and  thus  will 
not  be  included  here. 

Users  should  also  be  aware  that  all  flight  altitudes 
associated  with  these  data  sets  are  listed  in  meters  above 
ground  level,  m(AGL),  and  have  therefore  been  corrected 
for  local  meteorological  conditions  extant  at  the  time  of 
the  measurements. 


FLIGHT  C-379 
Rodby  Track 


Upper  and  Lower  Hemisphere 
2R2m  AGL  Filler  4 


Upper  and  Lower  Hemisphere 
6265m  AGL  Filter  4 


FI*.  2-3.  Right  C-379  photographs  From  AFGL-TR-81-0317. 
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Table  2.3.  Summary  of  iky  and  terrain  radiance  measurements 
(Table  1.2  from  AFQL-TR-84-OOW) 
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Ttbla  2.4.  Summary  of  atmospheric  vertical  profile  measurements. 
(Table  1.3  from  AFGL-TR-84-O096). 
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Table  2.S.  OPAQUE  related  aircraft  data  reports, 
(Table  1.1  from  AFGL-TR-84-0096). 
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Fig.  2-4.  Flight  C-379  (Rodby,  Denmark),  clear  sky  over  ocean  with  glitter,  300m  AGL. 
(Fig,  3-7  from  AFGL-TR-8I-03I7). 


Table  2.4.  Spectral  characteristics  summary  for  AVIZ  filters. 
(Table  2,1  from  AFGL-TR-84-0096), 
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Fl|.  2-S.  Standard  spectral  responses,  Peak  wavelengths  are: 
2-47Jnm  (Blue),  3«660nm  (Red),  4-SS0nm  (Photoplc), 
5-7SOnm  (NIR).  (Fig.  2-1  from  AFGL-TR-  84-0096), 


Sky  and  Terrain  Radiance  Data 

As  discussed  in  Johnson  (1981)  and  Johnson  and 
Hering  (1981),  both  of  which  should  be  considered  com¬ 
panion  reports  to  this  current  effort,  a  data  set  containing 
nearly  500  arrays  representing  measurements  of  sky  and 
terrain  radiance  values  has  been  developed  for  general 
application  to  image  propagation  studies  within  the  lower 
troposphere.  The  arrays  have  been  organized  for  storage 
by  flight  number  primarily,  although  date  and  site  are  also 
listed  in  the  header  information  as  will  be  noted  in  Sec¬ 
tion  3,  As  a  general  rule,  the  scanner  data  arrays  are 
organized  as  illustrated  in  Table  2.7. 

Each  of  the  individual  arrays  indicated  in  Table  2.7 
contains  1080  data  points  representing  either  an  upper  or 
lower  hemisphere  radiance  distribution.  Each  data  point 
represents  the  apparent  radiance  of  the  sky  or  terrain  as 
"seen"  by  the  5®FOV  of  one  of  the  airborne  scanning 
radiometers  described  in  Johnson  (1981).  There  is  a  data 
point  for  every  6°  in  azimuth  and  every  5®  In  zenith  angle 
within  each  hemisphere.  All  azimuth  notations  are  taken 
with  respect  to  the  sun  as  illustrated  in  Fig.  2-6.  It  should 
be  noted  that  for  each  entry  in  the  data  summary  listed  in 
Table  2.3,  there  are  two  data  strays  in  the  Horary,  one  for 
the  upper  hemisphere  radiances,  and  one  for  the  simul¬ 
taneously  measured  lower  hemisphere  radiances, 

The  radiance  values  included  in  this  data  library  are 
the  same  as  those  used  to  create  the  graphical  displays 
included  in  Johnson  (1981),  and  Johnson  and  Hering 
(1981),  and  should  be  used  with  the  same  caveats, 


Table  2.7.  Sky  &  Terrain  Radiance  Data  Organization  Per  Flight. 
(Table  2.2  from  AFGL-TR-84-0096. 

I.  Flight  Identification:  Flight  No.,  Date,  Site,  etc. 

A.  Altitude  Number  One  (See  Table  1.2) 

1.  Spectral  Band  No.  2  (see  Fig.  2-1) 

a.  Upper  Hemisphere  Radiance  Array 

b.  Lower  Hemisphere  Radiance  Array 

2.  Spectral  Band  No.  3 

a.  Upper  Hemisphere  Radiance  Array 

b.  Lower  Hemisphere  Radiance  Array 

3.  Spectral  Band  No.  4 

a.  Upper  Hemisphere  Radiance  Array 

b.  Lower  Hemisphere  Radiance  Array 

4.  Spectral  Band  No.  5 

a.  Upper  Hemisphere  Radiance  Array 

b.  Lower  Hemisphere  Radiance  Array 

B.  Altitude  Number  Two 

1.  ) 

ij'  k  same  as  In*  A*  above 

A.  i 

C.  Altitude  Number  Three 
I.  I 

j  t  same  as  In  "A"  above 

D.  Altitude  Number  Four 
2  1 

,  k  same  as  In  "A"  above 


The  software  required  to  create  these  displays  Is  not 
included  as  part  of  this  report. 


A 


THE  DATA  ARRAYS  CONTAINED  IN  THIS  LIBRARY 
HAVE  NOT  BEEN  CORRECTED  FOR  THE  NEAR  SUN 
STRAY  LIGHT  EFFECTS  AS  RECOMMENDED  IN  SECTION 
4  OF  JOHNSON  (1981),  AFGL-TR-81-0275. 

An  example  of  the  clear  day  sky  and  terrain  radi¬ 
ance  distributions  that  may  be  extracted  from  these  data 
arrays  is  shown  In  Appendix  A  for  flight  C-379.  These 
graphical  displays  have  been  abstracted  from  Johnson  and 
Hering  (1981)  as  a  convenient  reference  for  the  reader. 


Fig.  2-6.  Sky  and  terrain  coordinate  system.  (Fig,  2-2  from 
AFGL-TR-84-0096):  w-',-t  > 


Vertical  Profile  Data 

During  each  of  the  data  flights  listed  in  Table  2.3, 
there  were  radiometric  and  meteorological  measurements 
conducted  at  fixed  altitudes  as  well  as  during  ascent  and 
descent  modes.  The  data  which  were  associated  with  the 
ascents  and  descents  of  the  aircraft  have  been  processed 
and  reported  in  a  series  of  preceding  reports  which  are 
listed  in  Table  2.5.  However,  most  of  these  ascent  and 
descent  data,  l.e.  profile  data,  are  also  included  as  part  of 
the  general  purpose  library  described  herein.  These  profile 
data,  summarized  in  Table  2.4,  are  contained  on  a 
separate  data  tape  from  the  radiance  arrays  and  thus  can 
be  conveniently  used  either  with  or  without  the  larger 
radiance  data  set. 

As  with  the  scanner  radiance  data,  these  profile  data 
also  have  been  organized  for  storage  by  flight  number. 
The  data  available  for  each  flight  is  summarized  in 
Table  2.8.  Due  to  a  variety  of  data  processing  problems, 
there  is  not  an  exact  one  to  one  correspondence  between 
the  scanner  and  profile  data  sets  listed  in  Tables  2.3  and 
2.4.  The  raw  data  for  exact  correspondence  exists  on  ori¬ 
ginal  storage  tapes,  but  has  not  been  fully  retrieved  for 
inclusion  In  this  library. 

It  should  be  noted,  as  discussed  in  the  earlier  report 
series  listed  in  Table  2.5,  that  the  profile  measurements 
were  often  made  over  a  substantial  period  of  time,  and 
thus  represent  the  spatial  and  temporal  variabilities  occur¬ 
ring  along  a  flight  track  roughly  20  to  30  miles  long,  over 
a  period  of  30  to  120  minutes  in  time.  A  brief  description 
of  the  flight  profile  used  to  collect  these  data  has  been 
abstracted  from  AFQL-TR-80-0207,  Johnson  and  Oordon 
(1980)  and  is  reproduced  as  Appendix  C. 

The  radiometric  portion  of  the  profile  data  contains 


Tabic  2.8.  Scattering  coefficient  &  related  meteorological  data 
organization  per  flight.  (Table  2.3  from  AFGL-TR-84-0096) . 

I.  Flight  Identification:  Flight  No.,  Date,  Site,  etc. 

A.  Profile  No.  I  (Ascent)  (See  Appendix  C,  Fig,  C-l) 

1.  Scatt.  CoefT.  In  Spectral  Band  No.  2  (See  Fig,  2-1) 

2.  Ambient  Temperature 

3.  Dewpolnt/Froitpoint  Temperature 

4.  Relative  Humidity  (computed) 

5.  Absolute  Pressure 

6.  Atmospheric  Density  (computed) 

7.  Real  Time 

B.  Profile  No.  2  (Descent) 

I.  Scatt.  CoefT.  in  Spectral  Band  No.  3 


C. 


?} 
Profit* 
I.  Set 

?} 

Profllt 

1.  Set 

?} 


same  as  In" A"  above 

Profile  No.  3  (Ascent) 

I.  Scan.  CoefT.  in  Spectral  Band  4 

same  as  In  *  A"  above 


D.  Profile  No,  4  (Descent) 

1.  Scan.  CoefT.  In  Spectral  Band  3 

same  as  In  "A"  above 


measurements  of  total  volume  scattering  coefficient  as  a 
function  of  altitude  in  each  of  the  spectral  bands  indicated 
in  Table  2.4. 

The  meteorological  portion  of  the  profile  data  con¬ 
tains  measurements  of  ambient  temperature, 
dewpoint/frostpoint  temperature,  and  absolute  pressure. 
The  meteorological  measurements  were  made  simultane¬ 
ously  with  the  radiometric  set,  and  are  thus  listed  in  asso¬ 
ciation  with  a  specific  spectral  band  for  temporal  synchron¬ 
ization,  although  real  time  data  at  one  second  intervals 
were  also  included  in  the  original  raw  data  sets. 

In  both  sets,  the  specific  altitude  above  ground  level 
has  been  calculated  from  measurements  of  absolute 
atmospheric  pressure,  and  the  mean  ground  elevation  as 
determined  from  local  navigation  charts. 

Although  the  profile  data  points  are  all  listed  at  even 
30m  altitude  intervals,  users  should  be  aware  that  the  data 
values  represent  the  average  of  several  measurements 
taken  throughout  that  interval.  Scattering .  coefficient 
measurements  were  collected  at  a  rate  of  four  samples  por 
second,  and  meteorological  data  were  collected  at  a  rate  of 
two  samples  per  second  during  all  ascents  and  descents. 
The  rircrafts  average  rate  of  ascent/descent  throughout 
the  lower  3  km  AGL  was  approximately  1200  feet  per 
minute  (—6  meters/sec). 

Also  included  in  the  profile  data  listings  are  values 
for  several  derived  quantities.  Relative  humidity  has  been 
calculated  from  the  ambient  and  dewpoint/frostpoint  tem¬ 
perature  measurements,  and  atmospheric  density  lias  been 
calculated  from  the  temperature  and  pressure  measure¬ 
ments. 

*.n  example  of  the  measured  profile  data  for  flight 
C-379  is  illustrated  in  Appendix  B.  These  plots  of  scatter¬ 
ing  coefficient,  ambient  temperature  and  computed  relative 
humidity  have  been  abstracted  from  AFGL-TR-77-0078, 
Duntley  el  al.  (1977).  The  software  for  producing  these 
plots  Is  not  included  as  part  of  this  report. 

Users  should  also  note  that  the  profile  listings  in 
this  library  do  not  Include  extrapolations  above  or  below 
the  last  measured  data  point.  Thus  those  data  points 
reported  in  previous  reports  (Table  2.5)  which  were 
parenthetical  l.e.  representing  extrapolations,  are  not 
present  In  this  data  library.  Deletions  from  the  original 
data  base  for  instrumentation  reasons  are  likewise  zeroed 
out,  and  have  not  been  re-interpolated. 

THE  DATA  LISTINGS  CONTAINED  IN  THIS  LIBRARY 
HAVE  NOT  BEEN  CORRECTED  TOR  THE  HIGH  ALTITUDE 
BIAS  INDUCED  BY  INCOMPLETE  AERODYNAMIC  PURG¬ 
ING  OF  THE  AIRBORNE  INTEGRATING  NEPHELOMETER 
AS  DISCUSSED  IN  SECTION  2.3  OF  JOHNSON  AND  FITCH 
(1911),  AFGI -TR-S2-0049. 
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3.0  A  REVIEW  OF  FASCAT  DEVELOPMENT 
AND  STATUS 

The  development  of  techniques  for  the  practical 
determination  of  the  atmospheric  effects  on  visible  image 
transmission  is  a  central  aspect  of  the  optical  measurement 
and  analysis  program.  The  objective  is  to  establish  com¬ 
putational  techniques  that  are  general  enough  and  fast 
enough  to  cover  a  broad  span  of  real-time  field  applica¬ 
tions,  yet  accurate  enough  to  utilize  effectively  all  available 
data  that  are  relevant  to  the  calculations.  Studies  over  the 
past  few  yeais  have  produced  a  fast  atmospheric  scattering 
model  (FASCAT)  for  direct  determination  of  the  apparent 
spectral  radiance  of  distant  objects  and  backgrounds.  It 
represents  a  natural  extension  to  the  overall  Visibility 
Laboratory  program  of  theoretical  and  experimental 
research  into  the  structure  and  behavior  of  atmospheric 
optical  properties. 

The  development  and  testing  of  the  basic  multiple 
scattering  model  for  calculating  background  radiance  fields 
from  conventional  meteorological  observations  and 
environmental  data  are  described  in  two  previous  reports 
(Hering,  1981,  and  Hering,  1983).  For  ready  reference, 
the  second  of  these  reports  is  included  as  Appendix  A  of 
this  report.  It  provides  background  information  on  the 


modeling  procedures  established  for  the  contrast  transmit¬ 
tance  calculations  and  describes  performance  tests  involv¬ 
ing  comparisons  with  calculations  derived  from  other  radi¬ 
ative  transfer  models  and  comparisons  with  measured 
radiance  data  gathered  in  a  wide  range  of  environmental 
conditions. 

Recent  refinements  and  additions  to  FASCAT  have 
increased  the  scope  of  applications  that  can  be  addressed 
with  the  modeling  procedures.  Building  upon  techniques 
for  radiance  field  calculations  in  clear  and  overcast  sky 
conditions,  approximate  solutions  are  introduced  for  deal¬ 
ing  with  the  determination  of  image  transmission  charac¬ 
teristics  in  partly  cloudy  atmospheres.  As  described  in 
Section  3.1,  solutions  are  obtained  as  a  function  of  cloud 
amount  for  the  conventional  types  of  high,  middle,  and 
low  clouds. 

In  a  model  extension  described  in  Section  3.2,  ana¬ 
lytic  techniques  are  Introduced  for  the  direct  calculation  of 
Inherent  target  radiance  as  determined  from  the  intrinsic 
target  properties  and  the  calculated  irradiance  distribution. 
When  combined  with  calculations  of  directional  path  radi¬ 
ance  and  beam  transmittance,  this  model  extension  pro¬ 
vides  companion  calculations  of  the  apparent  spectral  radi¬ 
ances  of  both  target  and  background  along  any  predeter¬ 
mined  slant  path,  and  in  turn  the  apparent  target  contrast 
at  the  point  of  observation. 

Other  modifications  to  FASCAT  Involve  substantial 
revisions  to  data  input  procedures  in  an  attempt  to  sim¬ 
plify  the  entry  of  meteorological,  optical  and  target  param¬ 
eters  into  the  computer  program.  Detailed  step  by  step 
Instructions  for  the  .Input/output  data  files  are  included  In 
the  set  of  user  instructions  for  FASCAT  given  In 
Appendix  B.  A  brief  discussion  of  some  special  input 
options  is  given  in  Section'3.3  including  the  use  of 
LOWTRAN  aerosol  models  (Shettle  and  Fnnn,  1979)  for 
the  specification  of  optical  property  profiles. 

Finally,  in  Section  3.4,  calculations  of  visual  detec¬ 
tion  range  for  simplified  targets  and  backgrounds  are  used 
as  a  yardstick  to  explore  the  relative  impact  of  natural 
variations  in  atmospheric  and  target  factors  on  visible 
Image  transmission.  The  examples  illustrate  the  variation 
In  the  target  detection  range  as  a  function  of  changes  in 
boundary  layer  extinction  coefficient,  cloud  cover,  solar 
zenith  angle,  type  of  aerosol,  and  target  Illumination. 

3.1  Extension  of  the  Radiance  Model  to 
Partly  Cloudy  Atmospheres 

Explicit  determination  of  the  instantaneous  upper 
and  lower  hemisphere  radiance  fields  In  a  partly  cloudy 
atmosphere  requires  a  detailed  knowledge  of  the  actual 
cloud  structure,  orientation  and  illumination  that  is  nor¬ 
mally  not  available.  Even  if  precise  determinations  were 
possible,  the  resultant  Hculations  may  be  representative 
foi  nly  short  time  periods.  However,  we  can  consider  the 
development  of  computational  schemes  that  provide  esti¬ 
mates  of  the  average  (most  probable)  effects  of  changes  In 
cloud  cover  on  the  apparent  target  and  local  background 
radiance  fields.  For  partly  cloudy  atmospheres,  the  calcu¬ 
lation  of  both  the  direct  sun  contribution  to  the  irradiance 
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incident  on  the  target  and  the  singly  scattered  sunlight 
component  of  the  path  radiance  can  be  handled  as  before 
with  the  basic  FASCAT  model.  Cloud-free  path  algo¬ 
rithms  (Allen  and  Malick,  1983)  define  the  approximate 
probability  that  the  target,  background  or  observation  path 
will  be  in  sunlight  or  cloud  shadow.  The  important  com¬ 
plicating  factor  is  of  course  the  determination  of  the 
diffuse  components  of  the  irradiance  and  radiance  fields 
and  their  dependence  on  cloud  type  and  amount.  After 
consideration  and  rejection  of  many  approaches,  we  settled 
on  a  highly  simplified  method  which  solves  for  the  diffuse 
radiance  field  as  the  weighted  average  of  dear-sky  and 
overcast-sky  diffuse  field  calculations.  The  analytical  tech¬ 
niques  as  derived  through  reference  to  SOLMET  (1977) 
data  base  help  ensure  that  the  results  are  consistent  with 
the  average  observed  surface  irradiance  as  a  function  of 
cloud  type  and  amount  and  solar  zenith  angle.  The  SOL- 
MET  dat«i  base  and  its  application  for  the  determination  of 
the  average  optical  thickness  of  overcast  cloud  layers  for 
individual  cloud  types  are  described  in  Appendix  A,  Sec¬ 
tion  4.  Reference  is  made  also  to  the  discussion  of  the 
development  and  testing  of  techniques  for  the  determina¬ 
tion  of  the  radiance  fields  in  overcast  sky  conditions  given 
in  Appendix  A.  Extension  of  the  techniques  for  use  in 
partly  cloudy  atmospheres  is  discussed  in  the  following 
paragraphs. 

The  SOLMET  data  were  utilized  by  Shapiro  (1982) 
in  a  different  way  with  other  modeling  procedures  to 
determine  the  average  flux  transmittance  and  average 
albedo  for  individual  layers  of  partial  cloud  cover  for  many 
cloud  types  as  a  function  of  cloud  amount  and  solar  zenith 
angle.  With  additional  simplifying  assumptions,  the 
results  of  these  studies  can  be  applied  in  part  for  the  cal¬ 
culation  of  the  approximate  average  diffuse  irradiance  and 
radiance  fields  in  partly  cloudy  atmospheres. 


exhibit  a  systematic  behavior,  but  depart  significantly  from 
a  straight  linear  dependence  on  n.  Shapiro  derived 
bi-quadratic  polynomial  expressions  which  provide  close 
approximations  of  the  derived  weighting  factors  as  deter¬ 
mined  for  the  individual  cloud  types.  With  due  considera¬ 
tion  to  the  expected  time  and  space  variability  of  the 
albedo  of  partly  cloudy  layers  and  the  lack  of  strong  sys¬ 
tematic  variability  in  f(n,/ut)  as  a  function  of  cloud  type, 
we  combined  solutions  to  obtain  a  general  polynomial 
representation  given  by 

F(n,fi,)  -  n(  1.43- 1.21m,  -2.0Q« 

+  1.21  a**,  +  l.57/iJ ).  (3) 


It  should  be  emphasized  that  while  Eq.  (3)  as  depicted  in 
Fig.  3.1  applies  to  all  cloud  types,  the  albedo  of  the  partly 
cloudy  layer,  as  calculated  from  Eq.  (2),  is  dependent 
upon  the  actual  optical  properties  of  the  cloud  used  for  the 
calculation  of  A0y.  (See  Appendix  A,  Section  4).  Note 
also  that  F(n,n,)  is  zero  for  clear  sky  and  one  for  over¬ 
cast  cloud  cover. 

Since  the  total  downward  (solar  plus  diffuse)  irradi¬ 
ance  E(zt,d)  incident  on  the  top  of  a  single  cloud  layer 
does  not  vary  significantly  with  n,  we  can  write  approxi¬ 
mately  through  combination  of  Eqs.  (1),  (2),  and  (3), 

EKiza.u)  -  F(n,n,)  Epy  izg  ,u) 

+  (1-F(«,M,»£c*  <*.«>,  (4) 

where  E0y  (z„  ,  u  )  and  ECr  (r, ,  u )  are  the  upward  diffuse 
irradiances  as  determined  for  overcast  and  clear  sky 
respectively. 


DETERMINATION  OF  UPWARD  IRRADIANCE 

The  albedo  at  altitude  z„  at  the  top  of  the  cloud 
layer  is  given  by 

A  (z, )-  £(v  u)/E(za, d)  ,  (1) 

where  E(z<,,u)  and  Eiza,d)  are  the  upward  and  down¬ 
ward  irradiance,  respectively.  Shapiro  (1982)  assumed 
that  the  albedo  of  the  partly  cloudy  layer,  AKt  is  given  by 
the  weighted  average  of  the  albedo  for  a  corresponding 
clear  atmosphere,  ACr  ,  and  the  albedo  for  the  overcast 
sky,  A0v ,  having  the  same  cloud  form,  thus 


Arc  (*» )  “  Fin  )  Aqv  iz„ ) 

+  [I-FOi.m,  >]/<<:*<*<,) 


(2) 


Selecting  SOLMET  data  subsets  with  a  single  cloud  type 
and  amount  present,  Shapiro  determined  weighting  factors 
which  define  average  values  of  F(n,n ,)  as  a  function  of 
cloud  type,  cloud  fraction,  n ,  and  cosine  of  the  solar  zen¬ 
ith  angle  m , ■  The  derived  empirical  values  of  F(n,n,) 


DETERMINATION  OF  THE  DOWNWARD  IRRADIANCE 

A  similar  expression  can  be  derived  for  the  total 
downward  irradiance  below  the  cloud  layer.  The  flux 
transmittance  rT(l)  of  cloud  layer  /  is  given  by 

rT(D-  Eiz>,d)/Eizc,d)  (5) 

where  E(zt,d )  is  the  downward  irradiance  at  altitude zb  at 
the  base  of  the  cloud  layer.  Assuming 
Epc  (z„,d)  -  E0y  ( za,d )  -  ECr  iza,d)  we  have, 


EKi*b  *d)  “  Fin ,ii,)  Eayiit,  ,d) 

+  I  \-F\n,n,)]  ECr  izb,d) .  (6) 

The  total  downward  irradiance  can  be  expressed  as  the 
sum  of  the  diffuse,  pEizh^d),  and  direct  sun,  ,E(zb,d), 
components  hence, 


E<jt,d)-  DEizt,d)  +  ,Eizb,d) .  (7) 
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Let  us  assume  that  the  average  over  time  of  the  solar  irra- 
diance  component  in  partly  cloudy  atmospheres  is  deter* 
mined  by  the  cloud-free  solar  path  probability,  G(n,/*i) . 
such  that 

i Epcbh'd)  “  it, )  ,Ech  U^'d) 

+  l\-G(n,n,)\,Eov .  (8) 

Substituting  Eqs,  (7)  and  (8)  in  Eq.  (6)  and  assuming  that 
the  diffuse  irradiance  is  represented  by  the  time-averaged 
diffuse  irradiance,  we  have  for  partly  cloudy  skies 

C)Efc(ib,d)  ■  0Eoy(t),d) 

+  t  1-F(», nt)  1  pEct  (?b,d) 

+  l  \-F(n,tt,)-G(n^t)U,EcK^,d)-,Eoy%,d)]  .  (9) 

In  the  special  case  where  l-F—Q,  the  downward  diffuse 
irradiance  below  the  cloud  layer  is  given  by  the  weighted 
average  of  the  diffuse  irradiance  calculated  for  the 
corresponding  overcast  and  clear  atmospheres  as 
prescribed  by  Eq.  (4)  for  the  upward  irradiance  above  the 
cloud  layer.  However,  using  the  approximation  of  Allen 
and  Malick  (1983)  for  GOi.Mi)  and  Eq.  (3)  for  F(n,M »), 
we  And  that  the  quantity  (1-F)  is  invariably  larger  than  G 
ir  association  with  an  additive  component  to  the  downwel- 
ling  diffuse  irradiance  below  the  cloud  layer.  Likely 
sources  of  the  enhancement  of  difftise  irradiance  are  the 
forward  scatter  of  direct  solar  radiance  through  the  thin 
cloud  edges  and  the  reflection  downward  from  the  sides  of 
opaque  clouds.  The  additive  term  increases  with  increas¬ 
ing  solar  elevation  angle  and  increases  as  the  cloud  cover 
increases  from  scattered  to  broken.  Illustrations  of  these 
effects  through  sample  model  calculations  of  diffuse  radi¬ 
ance  are  given  later  in  this  section. 

DETERMINATION  OF  THE  AVERAGE  PATH  RADIANCE 

Conceding  for  reasons  discussed  above  that  detailed 
calculations  of  the  upper  hemisphere  radiance  fields  in 
partly  cloudy  atmospheres  are  not  practicable,  let  us  intro¬ 
duce  assumptions  that  greatly  simplify  calculation  of  the 
average  directional  path  radiance  distribution  that  is 
expected  to  occur  with  a  given  cloud  cover  and  aerosol 
distribution.  A  necessary  step  is  to  develop  approximate 
expressions  that  specify  the  expected  changes  in  the  path 
function  at  all  points  in  the  atmosphere  that  are  associated 
with  changes  in  cloud  amount  as  a  function  of  cloud  type 
and  solar  angle.  The  approach  lure  is  to  establish  general 
expressions  for  the  path  function  that  are  consistent  with 
the  average  behavior  of  the  diffuse  irradiance  fields  in 
partly  cloudy  conditions  as  represented  by  Eqs.  (4)  and 
(9). 

For  immediate  reference,  let  us  recall  the  basic 
expressions  for  path  radiance  and  path  function  given  in 
Appendix  A,  Section  2,  The  spectral  path  radiance  is 
given  by 


9 

where  £.(*',  8,*)  is  the  path  function,  defined  as  the 
point  function  component  of  the  path  radiance  that  is  gen¬ 
erated  by  the  scatteting  of  light  reaching  that  point  of  the 
path  from  all  directions.  The  path  function  can  be 
expressed  in  terms  of  its  contributions  from  the  scattering 
of  solar  (lunar)  scalar  irradiance  L-,(z,0,4)  and  from  the 
scattering  of  sky,  cloud  and  earth  radiances  £v(z,0,^> 


thus 

£•(1,0,4)  -  2,.,  (1,0,4)  +  Lv Or, 0,4) 

(11) 

where 

£.,(1,0.4) -,«(*)  F(z,0)i(z) 

(12) 

and 

£*(x,0,4>-  f  £(*,0\ 4')  FM)*(*Wn\ 

(13) 

/>(z,0)  is  the  normalized  single  scattering  phase  function 
at  angle  0  between  the  path  of  sight  and  the  direction  of 
the  source  light  and  i(z)  is  the  total  volume  scattering 
coefficient,  As  described  in  Section  2.1,  Appendix  A,  the 
delta-Eddlngton  approximation  is  employed  to  determine 
the  multiple  scattering  component  of  the  path  function. 
The  resultant  expression  for  the  path  function  (Eq.  (17), 
App.  A)  can  be  written, 

L.(z,0, 4)~  ■*(*)  I  ),*<*) 

+  LDii)  + 1  £.0>(z)cos6  + 

I  L0»(z)cos4  1 ,  (14) 

where  g  is  the  asymmetry  factor,  and  L0,  La  and  LD-  are 
the  average,  vertical  asymmetry,  and  azimuthal  asymmetry 
components  respectively  of  the  delta-Eddington  path  func¬ 
tion  representation. 

Given  approximate  expressions  for  the  calculation 
of  the  diffuse  component  of  the  path  function, 
Lv(z,0,d)  we  can  proceed  with  the  FASCAT  model  to 
calculate  the  spectral  path  radiance  distribution  in  partly 
cloudy  atmospheres  in  the  same  way  as  for  clear  and  over¬ 
cast  atmospheres.  Let  us  specify  that  the  diffuse  com¬ 
ponent  of  the  path  function  will  respond  to  changes  in 
cloud  amount  and  solar  zenith  angle  in  direct  correspon¬ 
dence  with  the  changes  in  the  diffuse  irradiance  above  and 
below  the  cloud  layer.  Consistent  with  Eq.  (4),  we 
assume  for  all  points  above  the  cloud  layer  that 

kLi  ,8, 4)  ”  F (n ,(i, )  ovLy  Or*  ,0,  4) 

+  1 1— Fin  ,M|) )  cr  f-v  .0.  4)  (15) 

ahd  by  analogy  with  Eq.  (9),  we  assume  for  all  points 
below  the  cloud  layer  that 
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pc (tj  , fl ,  <4 )  —  F(n,n,)  oyL^Ui,  ,9,  $)  (16) 

+  1 1— Fin  ,n, )  1  cp  L\i  ,6,  4) 

+  ~  1 1  -  F  in ,  m,  )■ -  G  (n ,  M , )  1 1 ,  £cs  <*„  ,d  >■ - ,  E0  v  % ,  rf )  1 

tf 

where  the  weighting  factors  F(n,n,)  are  given  by  Eq,  (3). 
The  probability  of  a  cloud  free  solar  path,  Gin, 9, ),  can 
be  expressed  (Allen  and  Malick,  1983) 


In  G(/i,e,)- (1  + c*tanfl,)  In  ,  (17) 

where  9,  is  the  solar  zenith  angle 

and 

pn  -  1-n  (l+3/i )/4  (18) 

and 

c„  -  0.55 — «/2  .  (19) 

For  completeness,  we  assume  that  Eq.  (15)  also  applies  at 
all  points  along  segments  of  cloud-free  viewing  paths  that 
are  between  the  base  and  top  of  the  partly  cloudy  layer. 

Sample  model  calculations  of  the  diffuse  scalar  Irra- 
diance  as  a  function  of  cloud  amount  and  solar  angle  for  a 
point  below  an  altocumulus  cloud  layer  of  average  optical 
thickness  are  shown  in  Fig.  3.1.  The  values  represent  the 
calculated  diffuse  radiant  flux  arriving  at  the  point  from  all 
directions  about  that  point  located  at  3  km,  with  the  base 
altitude  of  the  cloud  layer  at  3.5  km.  The  values  are  nor¬ 
malized  by  dividing  by  the  diffuse  scalar  irradiance  as  cal¬ 
culated  for  corresponding  the  clear-sky  condition.  They 
depict  only  the  normalized  diffuse  scalar  irradiance.  If  the 
point  is  In  direct  sunlight,  the  solar  irradiance  must  be 


NORMALIZED  DIFFUSE  SCALAR  IRRADIANCE 


Fig.  3-1.  Model  calculations  of  the  average  diffuse  radiant  flux 
arriving  at  a  point  0.5  km  below  the  base  of  an  altocumulus  cloud 
layer  as  a  function  of  cloud  amount.  The  Irradiance  values  are 
normalized  by  dividing  by  the  diffuse  scalar  Irradiance  for  the 
clear  atmosphere  with  corresponding  aerosol  optical  properties. 


added  to  obtain  the  total  scalar  irradiance  reaching  that 
point.  For  solar  zenith  angles  less  than  that  about  65 
degrees,  the  calculated  average  diffuse  scalar  Irradiance 
tends  to  be  larger  for  partly  cloudy  conditions  than  for 
either  clear  or  overcast  conditions.  As  illustrated  in 
Fig.  3.1,  the  diffuse  scalar  irradiance  increases  with 
increasing  cloud  amount  up  to  about  7/10  cloud  cover  and 
decreases  thereafter  as  the  shadowing  effects  of  adjacent 
cloud  elements  become  more  prominent.  Notice  that  the 
normalized  diffuse  irradiance  increases  markedly  with 
increasing  solar  elevation  angle. 

For  calculations  of  the  time  or  space  averaged 
values  of  path  function  for  points  below  the  cloud  layer 
with  Eq.  (11),  the  increase  In  the  diffuse  component  with 
increasing  cloud  amount  tends  to  be  offset  by  a  decrease 
in  the  direct  solar  component  of  the  path  function  due  to 
the  increased  probability  that  the  points  along  the  observa¬ 
tion  path  will  be  in  cloud  shadow.  The  average  or  most 
probable  value  of  the  solar  component  of  the  path  func¬ 
tion  for  partly  cloudy  atmospheres  is  given  by, 

PcLt,  (2 ,9,$) m  Gin ,nt )  cr />•» (2,9,4) 

•Ml  -<?0i,m,)]ovL*.  (*,*,*).  (20) 


TWO-LAYER  CLOUD  OPTION 

The  above  system  of  analytic  expressions  yields  esti¬ 
mates  of  the  average  path  radiance  distribution  for  atmo¬ 
spheres  with  a  single  cloud  layer  of  varying  cloud  amount. 
A  simple  strategy  was  used  to  extend  the  FASCAT  model 
calculations  to  include  a  second  layer  of  fractional  cloud 
cover,  although  the  physical  basis  for  the  determination  of 
the  average  diffuse  spectral  radiance  is  much  more  tenu¬ 
ous.  The  assumption  is  made  that  the  weighting  factors 
F(n,n,)  as  determined  for  a  single  cloud  layer  will  apply 
in  general  to  conditions  where  the  layer  exists  in  combina¬ 
tion  with  a  second  layer  of  variable  cloud  cover.  In  the 
absence  of  representative  experimental  data  for  definitive 
tests  of  this  assumption,  the  2-layer  cloud  option  is 
included  in  FASCAT  only  as  a  means  to  obtain  rough  esti¬ 
mates  of  the  multi-layer  cloud  effects  on  image  transmis¬ 
sion.  The  diffuse  path  function  component  is  calculated  as 
the  weighted  average  of  the  corresponding  values  as  calcu¬ 
lated  for  the  overcast  upper  layer,  U,  the  overcast  lower 
layer,  L,  overcast  both  layers,  UL,  and  for  the  clear  atmo¬ 
sphere,  CR,  For  the  region  above  the  base  of  the  upper 
cloud  layer,  we  assume 

p^L^ii ,9 ,4)  “  I  )~F(/ iny,i i, )  1  (  l-Fy  ( n y  ,n,  1  ch L'j  i2 ,9, 4 ) 
+  Fy iny,n, )  [  1 —Fl  iny tMi ) )  ot'.i/f'vlr.fl,  4 ) 

■t-  F^  iny  ,ti, )  [  I—  Fu  (ny,fi, )  1  ova  L^is  ,9, 4) 

+  Fy  («(/  ,/i,  )  Fy  {fly  ,/U,  )  QV  ML  (*  <9,  4 )  (21) 


For  the  region  below  the  base  of  the  upper  layer  we 
assume 
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kL^(z,6,  4)  “  1 1- Fy  Oiy,M»)  1  l  1 —Fi  (fly  ,Mt  )  cflLv  U,0, 


+  Fy(fly,M»)  [  1  I  ,4) 

+  fj,  (fl^  ,/t| )  (  1— Fy  (fly  ,M|)  1  L*i(z,Oy  ♦) 

+  Fy  (fly, Ml  )  Fy  (fly  ,M|)  ,9,  $) 

+  j  1  l-Fy(fly,M.)  ]  1  1-fi  <«i,M,)  ] 

~  Gy  (fly  ,Mi  )  Gy  (fly  ,M«)  | 

+  ~  - |  Fy  (fly  ,Mi )  (  l~Fy  (fly  .Mi )  1 

-  Gy  (fly  ,Mi )  I  l“Gy  0*0  >Mi )  1  J 

+  ~~~ —  |  fi  (»£  ,M» )  C  1— Fit/  (nU  iMj  )  1 

-  Gy  (fly ,/*, )  I  1- Gy  (fly  ,Mj)  1 J 

+  — ■  ’ -  |  Fy  (fly  ,Ml  )  fit  («£  >Mi ) 

-  [  1— Gy  (fly, Mi  )  1  t  1“  Gy  (fly  ,M|)  1  |  i  (22) 

where  the  factors  F{n,n,)  and  G(a,m>)  are  determined 
as  before  from  Eqs.  (3)  and  (17),  respectively,  except  that 
Gy(«y,M •)“  1  above  the  base  of  the  lower  cloud  layer. 
For  optically  thick  clouds,  ,£y,  ,EL  and  ,£yy  lend  to 
zero. 

Diagnostic  model  calculations  of  visual  detection 
range  as  a  function  of  cloud  amount  are  given  in  Sec¬ 
tion  3.3  below. 

3.2  Calculation  of  Apparent  Target  and 
Background  Radiances 

The  fundamental  equation  for  the  apparent  spectral 
radiance  of  distant  target  t  at  range  r,  as  given  by  Dunt- 
ley,  etal.  (1957),  [see  Eq.  (1)  Appendix  A]  is 

i Lr(z ,0 , 4>)  m  t LgiZf  ,0,'d)  Tr{z,9) 

+  l,*M,d),  (23) 

where  ,L0  is  the  inherent  target  radiance  at  altitude  z(,  and 

Tr  -  exp  -  J  a  (r  )dr  Is  the  path  transmittance  and  a  is 

the  volume  attenuation  coefficient.  The  apparent  radiance 
of  the  target  is  the  sum  of  the  residual  target  radiance  and 
the  path  radiance,  which  is  given  by  Eq.  (10).  We  note 
that  the  determination  of  the  apparent  target  radiance 
requires  resolution  of  3  inter-related  atmospneric  effects: 
the  influence  of  the  irradiance  distribution  on  the  inherent 
target  radiance,  the  attenuation  of  the  target  radMce,  and 
the  generation  of  the  spurious  path  radiance. 


«  v  v.  n-%  n »y  r-y »v  <r.  < j  rj  f 

k'»A  u.*-  k-av4*.»  j*  W»  >!>»>  »■«  b  ■»  W»-»  m  ■  1. 


For  the  spedflc  purpose  of  calculating  the  Inherent 
target  radiance,  it  is  assumed  that  both  the  downward  and 
upward  diffuse  radiance  fields  are  uniform.  Now  the 
diffuse  component  of  the  irradiance  on  the  target  surface 
can  be  expressed 


fl//,G,,0w)-o£Or,,d)  [  1+  cos«M  1/2 

+DE(*,,u)  [  1  -  cos 1/2  (24) 


where,  as  in  Eq.  (4)  and  Eq.  (7),  E(x,,u)  is  the  upward 
irradiance  on  a  horizontal  surface,  DE(*,,<t)  is  downward 
diffuse  irradiance  on  a  horizontal  surface  and  the  direction 
of  the  target  surface  normal  is  (0m>4m)-  The  total  (solar 
plus  diffuse)  irradiance  on  the  target  surface  is 


Af(  (*( ,0jf  m  <|  (if  )cos0n  +  d Hf  (?)  t9fft )  ,  ■  (25) 


where  the  angle  between  the  solar  direction  and  the  target 
surface  normal  direction  Is  given  by 


co*0*  -  COI0M  cosS, 

+  iln0M  sins,  cos(0M  -  4,) .  (26) 


The  cos#*  it  set  to  zero  if  cos9„<0  or  if  the  target  is  in 
local  shadow.  Finally,  it  is  assumed  that  the  target  sur¬ 
faces  reflect  In  accordance  with  Lambert’s  law  so  that  the 
expression  for  the  inherent  spectral  radiance  of  the  target 
becomes 


i L»(*i>®>^)  **  A,//, (x, ,9„ ,9^)/ir  (27) 

The  expressions  corresponding  to  Eqs.  (23),  (24), 
(25),  and  (27)  for  the  apparent  spectral  radiance  of  the 
local  background  are 


»Lr(i ,«,*)- */,(**, 9,*)  rf(z,tf)+  L;b,9,4>) ,  (28) 


yfyOfc.e,^)  -  DE(xt,d)  I  1  +  cos 8Nb  1/2 

+  c£(*»,m)  I  1  -  cosS/v*  1/2  ,  (29) 


//j(**,0,»,0M)-  «,( iA)cos0,4  +  DHb(zb.8Nb) ,  (30) 
and 

(zb  ,9,  #)  “  Eb Hb  (jb  ,9lb  ,9tfb  Vn  ,  (31) 

where  the  subscript  b  refers  to  the  local  background  sur¬ 
face  which  appears  in  immediate  contrast  with  the  desig¬ 
nated  target  surface.  As  in  the  case  of  the  target  surface, 
the  direction  of  the  normal  to  the  local  background  is  an 
independent  Input  variable  to  FASCAT.  For  objects 
viewed  against  a  sky  background,  is  zero.  If  one 
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selects  the  sea  surface  option  in  FASCAT,  Rb  is  ignored 
for  observing  paths  directed  toward  the  sea  surface,  and 
the  inherent  radiance  is  calculated  as  the  Fresnel 
reflectance  of  the  wind  ruffled  water  surface  assuming  the 
wave  slope  probabilities  of  Cox  and  Munk  (1954)  (see 
Section  3,  Appendix  A). 

3.3  FASCAT  Application  Methodology 

In  order  to  calculate  the  apparent  spectral  contrast  of 
distant  targets,  one  must  know  or  estimate  a  variety  of 
environmental  and  target  factors.  Needless  to  say,  the 
accuracy  and  representativeness  of  data  available  for  these 
determinations  varies  markedly  depending  upon  the  cir¬ 
cumstances  under  which  the  model  is  applied.  For  this 
reason  emphasis  has  been  placed  upon  the  development  of 
data  input  procedures  that  help  simplify  the  input  of  the 
required  information  and  are  adaptable  to  the  amount  and 
variety  of  available  data.  The  FASCAT  input  file  and  an 
item  by  item  set  of  user  instruction  notes  are  included  as 
Appendix  B  of  this  report.  Some  recent  changes  in  the 
input  format  and  some  special  features  of  the  modeling 
procedures  are  reviewed  in  the  following  paragraphs. 

The  basic  method  used  for  data  entry  has  been 
modified  to  improve  both  the  data  handling  and  the  com¬ 
putational  efficiency  of  the  model.  Under  the  new  format 
the  data  are  input  sequentially  in  four  functional 
categories:  general  parameters,  optica!  property  profiles,  cloud 
layer  lr\formatlon,  and  target  and  local  background  data. 
Selected  data  entry  items  in  each  category  will  be  discussed 
below  to  augment  the  user  instructions  given  in 
Appendix  B. 


GENERAL  PARAMETERS 


Wavelength  •  Although  the  model  calculations  are 
strictly  applicable  only  for  monochromatic  radiation,  the 
results  are  suitable  with  good  approximation  for  reason¬ 
ably  broad  spectral  bands  in  the  visible  and  near  infrared 
portions  of  the  spectrum.  The  designated  wavelength 
should  be  chosen  to  conform  as  closely  as  possible  with 


the  overall  spectral  characteristics  of  the  sensor.  The 
model  is  applicable  over  the  range  450  to  1060  nm. 


Average  surface  reflectance  -  The  weighted  spatial 
average  of  the  surface  reflectance  is  entered  to  calculate 
the  component  of  the  path  radiance  that  is  generated  by 
light  reaching  the  path  through  reflection  by  the  underly¬ 
ing  surface  and  subsequent  scattering  by  the  atmosphere 
in  the  direction  of  the  sensor.  Thus,  the  area  comprising 


the  spatial  average  may  be  small  as  in  the  case  of  slant 
observation  paths  near  the  earth's  surface  or  large  as  for 
paths  extending  from  high  altitude  to  the  horizon.  As 
described  under  target  data  below,  another  data  entry  is 
made  for  the  local  background  reflectance,  which  often 
differs  from  the  average  reflectance.  The  local  reflectance 
is  used  for  the  determination  of  the  inherent  radiance  of 


the  immediate  background  against  which  the  target 


appears. 


Base  altitude  of  the  top  layer  -  This  input  ite  n  sets 
the  lower  altitude  limit  of  the  primary  ozone  'ajer  and 


should,  in  meet  instances,  correspond  with  the  altitude  of 
the  tropopause.  Stored  climatological  values  of  total  ozone 
as  a  function  of  latitude  and  season  are  used  to  calculate 
the  ozone  optical  depth  of  the  top  layer  for  the  designated 
wavelength.  Desired  departures  from  the  climatological 
values  can  be  handled  through  appropriate  changes  in  the 
designated  single  scattering  albedo  for  the  top  layer.  The 
single  scattering  albedo  for  the  top  layer  normally  is  used 
to  define  only  the  absorption  by  molecular  or  aerosol  con¬ 
stituents  other  than  ozone. 

Observation  angles  -  Up  to  20  zenith  observing 
angles  may  be  selected.  Default  values  include  every  ten 
degrees  beginning  at  5  and  ending  at  175  degrees  plus  zen¬ 
ith  angles  of  100  and  110  degrees.  Three  azimuthal 
observing  angles  may  be  selected.  They  are  specified  in 
items  21  and  22  of  the  input  file  as  departures  from  the 
solar  azimuthal  angle.  The  default  values  are  0  (upsun), 
180  (downsun)  and  90  (cross  sun)  degrees. 


Under  the  revised  data  entry  format,  the  optical  pro¬ 
perty  profiles  for  the  cloud-free  atmosphere,  above,  below, 
and  between  cloud  elements,  are  established  first.  Next 
the  cloud  properties  are  superimposed  and  finally  informa¬ 
tion  with  respect  to  sensor,  target  and  background  factors 
are  entered. 

Optical  Properties  -  Depending  upon  the  complete¬ 
ness  of  available  measurements  and  observations,  one 
may  select  up  to  10  layers  to  define  the  vertical  distribu¬ 
tion  of  the  Individual  parameters.  For  each  designated 
layer,  the  optical  scattering  ratio,  the  single  scattering 
albedo,  and  the  single  scattering  phase  function  and  the 
associated  asymmetry  factor  must  be  entered.  Reference 
is  made  to  Sections  3.2  and  3.3  of  Appendix  A  for  a  dis¬ 
cussion  of  these  optical  parameters  and  for  a  description  of 
options  for  their  specification  and  entry  into  the  computer 
program. 

A  recent  extension  to  the  FASCAT  computer  code, 
developed  by  Lt.Col.  John  D.  Mill  and  Bric.P.  Shettle  of 
the  Air  Force  Geophysics  Laboratory,  enables  specification 
of  the  optical  properties  for  each  layer  through  use  of  the 
LOWTRAN  aerosol  models  (Shettle  and  Fenn,  1979). 
Representative  optical  properties  for  rural,  maritime, 
urban  and  tropospheric  atmospheres  as  well  as  for  advec- 
tion  and  radiation  fog  are  supplied  through  interpolation 
from  stored  look-up  tables.  The  aerosol  model  and  associ¬ 
ated  properties  are  designated  independently  for  each 
layer.  The  operator  enters  the  appropriate  aerosol  type, 
the  visible  (550  nm)  extinction  coefficient  and  the  relative 
humidity,  and  the  corresponding  values  of  optical  scatter¬ 
ing  rat  single  scattering  albedo  and  the  phase  function 
for  the  designated  wavelength  are  determined  and  supplied 
to  the  main  program  for  radiance  field  calculations. 

As  discussed  In  Section  3.3,  Appendix  A,  many 
options  are  available  In  particular  for  the  specification  of 
the  single  scattering  phase  function.  For  example,  the 
program  will  accept  the  coefficients  of  two-term  Henyey- 
Greenstein  functions  and  calculate  the  scattering  phase 
function.  When  other  information  is  lacking,  another 
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model  option  will  estimate  the  Henyey-Greenstein 
coefficients  from  the  scattering  ratio  using  algorithms 
derived  from  the  extensive  field  measurements  of  phase 
function  by  Barteneva  (1960). 

It  should  be  emphasized  that  profile  information  for 
alj  atmospheric  levels  must  be  entered  before  specific 
cloud  type  and  amount  information  is  entered  below.  For 
layers  of  fractional  cloud  cover,  the  data  entries  in  all 
cases  define  the  optical  properties  of  the  cloud-free  paths 
between  the  individual  clouds.  For  overcast  cloud  layers, 
the  operator  ha;  two  choices.  First  he  can  simply  make 
dummy  entries  for  the  optical  properties  for  the  overcast 
cloud  layer  for  the  initial  profiles  that  will  be  ignored  and 
new  data  will  be  supplied  through  entry  of  cloud  Informa¬ 
tion  in  the  input  file  items  discussed  below.  The  second 
option  for  overcast  layers  can  be  used  if  the  operator 
desires  to  enter  his  own  cloud  optical  properties  (scattering 
ratio,  single  scattering  albedo,  and  scattering  phase  func¬ 
tion)  rather  than  rely  on  the  average  values  prescribed  as  a 
function  of  cloud  type.  If  the  second  option  is  chosen  for 
overcast  cloud  representation,  the  cloud  optical  properties 
are  entered  along  with  and  in  the  same  manner  as  for  the 
cloud-free  layers,  and  no  entry  for  the  cloud  layer  is  made 
in  the  cloud  layer  input  file  below. 


The  input  for  the  individual  cloud  layers  (2-layer 
maximum)  consists  of  the  base  and  top  altitudes,  the 
cloud  amount  in  tenths,  the  cloud  type 
(cirrus/cirrostratus,  altostratus/altocumulus,  cumulus, 
stratus/stratocumulus,  or  nimbostratus/precipitation)  and 
the  relative  optical  depth  of  the  cloud.  The  relative  optical 
depth  has  3  options:  average,  thin  (minus  standard  devia¬ 
tion),  and  thick  (plus  standard  deviation).  In  a  separate 
phase  of  study  summarized  in  Section  4.1 ,  and  Table  3  of 
Appendix  A,  the  average  and  the  variability  of  cloud  opti¬ 
cal  depth  as  a  function  of  cloud  type  were  determined 
empirically  from  the  SOLMBT  data  base.  The  model  algo¬ 
rithms  return  the  optical  depth  for  the  specified  cloud  type 
and  relative  optical  thickness  (average,  thick  or  thin)  from 
the  table  without  direct  reference  to  the  specified  base  and 
top  altitude  of  the  layer.  In  other  words,  the  optical  thick¬ 
ness  is  determined  independently  of  the  specified 
geometric  thickness  of  the  cloud. 


Data  sets  consisting  of  sensor,  target,  and  local 
background  Information,  are  entered  in  the  last  section  of 
the  input  file.  Once  the  path  function  and  beam  transmit¬ 
tance  distributions  have  been  established  from  prior  entry 
data,  the  model  deals  in  turn  with  each  of  the 
target/background  data  sets  without  need  for  repeating  the 
basic  path  function  and  transmittance  calculations.  Data 
entry  for  all  cases  where  the  targets  are  viewed  from 
above  against  the  earth’s  surface  as  the  background  are 
listed  first,  followed  by  cases  where  the  objects  are  viewed 
from  below  against  a  sky  or  cloud  background.  The 
number  of  target/background  data  entries  is  unlimited. 
The  input  data  for  each  set  includes: 


Sensor  altitude  (km  mat) 

Target  altitude  (km  msl) 

Target  reflectivity 

Target  normal  zenith  angle  (deg.) 

Target  normal  azimuthal  angle  (deg.  departure 
from  solar  azimuth) 

Target  illumination  (3  options) 

(1)  full  sunlight 

(2)  cloud  shadow  •  includes  direct  solar  beam  and 
forward  peak  radiance  penetration  of  thin  clouds 

(3)  local  shadow 

Local  background  (against  which  target  appears) 
reflectivity 

Local  background  surface  normal  zenith  angle  (deg) 

Local '  background  surface  normal  azimuthal  angle 
(departure  from  solar  azimuth  angle) 

Local  background  illumination  (3  options) 

(1)  full  sunlight 

(2)  cloud  shadow 

(3)  local  shadow 


Data  entry  items  for  upward  paths  are  the  same  as 
for  downward  paths  of  sight  except  that  local  background 
input  data  are  not  required.  The  apparent  local  back¬ 
ground  radiance  for  a  clear  or  partly  cloudy  atmosphere  is 
the  calculated  sky  radiance  for  a  cloud-free  path.  For  the 
case  of  an  intervening  overcaat  layer,  apparent  radiance  of 
the  cloud  is  returned  as  the  local  background  radiance.  A 
separate  data  set  is  employed  to  calculate  the  apparent 
radiance  of  the  local  background  where  that  background  Is 
an  opaque  cloud  element  In  a  partly  cloudy  sky.  For  the 
second  data  set,  the  cloud  element  is  assumed  to  be  the 
target  for  the  upward  path  of  sight,  and  the  appropriate 
cloud  reflectivity,  cloud  surface  normal  direction,  and 
cloud  illumination  (full  sun  or  shadow)  are  entered  as  tar¬ 
get  data. 

It  should  be  noted  that  for  both  upward  and  down¬ 
ward  looking  paths,  the  choice  as  to  what  is  the  target  and 
what  is  the  local  background  against  which  the  target 
appears  can  be  arbitrary  In  many  situations,  Model  calcu¬ 
lations  of  the  apparent  radiance  of  opaque  surfaces  are  the 
same  regardless  of  the  target/background  assignment. 

3.4  Sensitivity  Analysts 

Special  attention  has  been  given  to  diagnostic  stu¬ 
dies  of  the  effects  of  uncertainties  in  the  input  data  on  the 
.  model  calculations  of  image  transmission  properties.  The 
results  of  a  broad  series  of  trial  sensitivity  calculations  are 
summarized  in  Section  O  of  Appendix  A.  Visual  detec¬ 
tion  range  was  used  as  the  yardstick  to  analyze  the  relative 
importance  of  the  changes  in  each  of  many  environmental 
and  target  factors  on  the  detection  of  distant  objects. 
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In  this  section  we  will  extend  the  analysis  to  include 
targets  in  shadow  as  well  as  in  sunlight  and  to  Include  tar* 
gets  and  local  backgrounds  with  non-horizontal  surfaces. 
The  following  calculations  are  based  on  the  same  average 
reference  atmosphere  (Ref.  A)  with  the  exception  that  the 
height  of  the  boundary  layer  extends  to  1.3  km  rather 
than  1  km.  As  before,  a  uniform  target,  30  m  in  diame¬ 
ter,  is  assumed,  which  is  viewed  from  aloft  against  an 
uncluttered  surface  background.  The  analytic  representa¬ 
tion  of  contrast  threshold  for  visual  detection  is  described 
in  Section  6.2  of  Appendix  A.  It  is  based  upon  the  data  of 
Taylor  (1964)  and  corresponds  to  a  99  percent  detection 
probability  for  a  glimpse  time  of  1/3  sec  and  a  lack  of 
knowledge  of  the  target  position  of  ±  4  deg  or  more. 
The  reflectivity  of  the  target  is  .15  and  the  background 
reflectivity  is  .075.  In  all  cases  the  target  is  at  sea  level 
and  the  sensor  altitude  is  3  km.  Unless  specified  other¬ 
wise,  the  wavelength  is  550  nm  and  the  solar  zenith  angle 
is  60  deg. 

The  results  of  the  supplementary  model  calculations 
are  shown  in  Figs.  3.2  to  3.5.  The  cartoon  inserts  in  each 
figure  depict  specific  conditions  with  respect  to  the  observ¬ 
ing  path  and  the  target  orientation  and  illumination.  In 
most  examples,  the  departure  of  the  observation  azimu¬ 
thal  angle  from  the  solar  azimuthal  direction  is  180  deg. 
(downsun).  The  target  surface  for  these  examples  is 
either  horizontal  or  vertical  and  the  target  and  local  back¬ 
ground  may  be  in  direct  sunlight,  in  local  shadow,  or  in 
cloud  shadow.  For  the  sensitivity  calculations,  all  input 
variables  other  than  the  factor  under  examination,  are 
held  constant  and  equal  to  the  reference  values  listed 
above. 

VARIATION  WITH  EQUIVALENT  SURFACE  VISIBILITY 

Fig.  3.2  shows  the  variation  in  horizontal  detection 
distance  as  measured  along  the  ground  surface  as  a  func¬ 
tion  of  the  scattering  ratio  in  the  boundary  layer.  The 
equivalent  surface  visibility  is  approximated  through 
Kochmieder's  law  under  the  assumption  of  complete  aero¬ 
sol  particle  mixing  within  the  boundary  layer  and  an 
assumed  threshold  contrast  of  5  percent.  Calculations  are 
shown  for  boundary  layer  scattering  ratios  of  52,  26,  13, 
and  7.5  which  correspond  to  equivalent  surface  visibilities 
of  5,  10,  20  and  40  km,  respectively. 

The  first  two  data  sets,  (a)  and  (b),  assume  that  the 
target  and  local  background  scene  are  in  direct  sunlight, 
The  detection  distances  are  somewhat  larger  for  the  verti¬ 
cal  target  surface  In  set  (b)  because  of  the  enhanced  solar 
irradlance  of  the  vertical  surface  (solar  zenith  angle  of  60 
deg.)  yields  a  larger  inherent  contrast  with  the  local  back¬ 
ground. 

Notice  the  substantial  reductions  In  detection  dis¬ 
tance  for  all  values  of  surface  visibility  for  data  set  (c) 
(scene  in  local  shadow)  and  data  set  (d)  (altocumulus 
overcast)  as  compared  with  the  sunlit  conditions  of  sets 
(a)  and  (b).  The  visual  detection  distances  are  a 
minimum  for  the  scene  in  local  shadow  (set  c)  where  the 
inherent  target/background  signal  is  small  in  the  shadowed 
area  and  the  path  radiance  remains  relatively  high  for  the 


sunlit  viewing  path,  The  detection  distances  are  larger  for 
the  intervening  overcast  condition  (set  d)  as  compared 
with  the  clear-sky,  local-shadow  case  (set  c)  for  two  rea¬ 
sons.  First,  the  downwelling  diffuse  irradlance  on  target 
and  background  is  larger  for  the  overcast  sky  condition, 
and  second,  the  deleterious  effects  of  path  radiance  are 
reduced  in  the  absence  of  the  direct  sun  component  of 
scattering  along  the  observing  path  under  the  overcast. 
The  latter  effect  is  offset  a  bit  by  the  increase  in  diffuse 
irradlance  of  the  path. 

VARIATIONS  WITH  SOLAR  ANGLE 

Shown  in  Fig.  3.3  are  examples  of  the  changes  in 
the  detection  distances  of  sunlit  targets  as  a  function  of 
changes  in  solar  zenith  angle.  The  model  calculations 
cover  the  range  of  solar  angle  corresponding  to  a  time 
span  from  early  morning  to  noon  in  middle  latitudes  in 
the  spring  or  fall  season.  Data  sets  (a)  and  (b)  are  Tor  a 
horizontal  and  vertical  target  surface,  respectively,  as 
viewed  against  a  horizontal  local  background.  Notice  that 
the  detection  distances  are  much  lower  for  the  low  solar 
elevation  angles  in  data  set  (c),  where  the  local  back¬ 
ground  is  Inclined  30  degrees  toward  the  sun.  Since  the 
reflectivity  of  the  local  background  for  all  model  calcula¬ 
tions  Is  assumed  to  be  about  one-half  of  the  target 
reflectivity,  the  prescribed  change  in  the  slope  and  orienta¬ 
tion  of  the  local  background  in  set  (c)  serves  to  reduce  the 
inherent  target/background  contrast,  particularly  at  high 
solar  zenith  angles. 

VARIATIONS  WITH  CLOUD  COVER 

As  emphasized  In  Section  3.1,  the  model  calcula¬ 
tions  deal  only  with  the  average  dependence  of  the 
apparent  radiance  distribution  on  cloud  amount  in  partly 
cloudy  atmospheres.  Variations  in  detection  distance  with 
changes  in  cloud  cover  from  clear,  to  scattered,  to  broken, 
to  overcast  sky  conditions  are  shown  in  Fig.  3.4  for  three 
target/background  scenes.  Data  set  (a)  for  a  horizontal 
target  surface  in  sunlight  and  set  (b)  for  a  vertical  target 
surface  in  sunlight  respond  in  a  similar  way  to  changes  in 
cloud  amount.  The  calculated  detection  distances  remain 
essentially  the  same  in  each  of  these  cases  for  sunlit  target 
scenes  with  clear,  0.5,  and  0.8  cloud  cover.  The  lack  of 
sensitivity  of  average  detection  distance  to  changes  in 
cloud  amount  is  the  result  of  compensation  in  the  atmos¬ 
pheric  effects.  The  increase  in  the  diffuse  atmospheric 
irradlance  with  increasing  cloud  cover  for  the  sunlit  target 
produces  an  increase  in  the  inherent  target  contrast  signal 
but  also  serves  to  increase  the  troublesome  path  radiance. 
On  the  other  hand,  as  the  cloud  amount  increases  the  pro¬ 
bability  that  points  along  the  viewing  path  beneath  the 
cloud  layer  will  be  in  sunlight  decreases  and  consequently 
the  average  observed  path  radiance  is  diminished.  Fluc¬ 
tuations  in  cloud  amount  also  affect  the  spatial  average  of 
the  downwelling  irradlance  at  the  surface  and  in  turn  affect 
the  surface  reflected  component  of  Irradlance  reaching  the 
target  surface,  the  local  background  surface,  and  the 
observation  path  below  the  cloud  layer.  The  results  for 
the  particular  reference  conditions  used  for  the  model  cal¬ 
culations  in  Fig.  3.4  show  a  modest  systematic  increase  in 


Variation  With  Surlaoa  Visibility  (Km) 

Fig.  3-2.  Model  calculations  of  visual  detection  distance  as  a  function  of  the  optical  depth  or  the  boundary  layer 
(expressed  as  equivalent  surface  visibility).  Oround  based  target  at  sea  level  is  viewed  downsun  (<A-  180 degrees)  from 
3  km  altitude.  Scene  (a)  •  horizontal  target  surface  in  sunlight.  Scene  (b)  •  vertical  target  surface  in  sunlight,  Scene 
(c)  -  horizontal  target  surface  in  local  shadow,  Scene  (d)  •  horizontal  target  surface  under  altocumulus  overcast  with 
base  altitude  of  3.5  km, 
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Variation  With  Solar  Angle 

Fig.  3-3.  Visual  detection  distance  vs  solar  zenith  angle  (deg),  Target  at  sea  level  Is  viewed  downsun  ( 4>  -  180“)  from 
an  altitude  of  3km.  Scene  (a)  ■  horizontal  target  surface  In  sunlight.  Scene  (b)  -  vertical  target  surface  in  sunlight. 
Scene  (c)  ■  same  as  (a)  except  the  local  background  surface  Is  inclined  30°  toward  the  solar  direction. 
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Variation  With  Cloud  Covar 

Fig.  >4.  Visual  detection  distance  vs  fractional  amount  of  altocumulus  cloud  with  base  altitude  of  3.5  km.  Sea  level 
target  Is  viewed  downsun  from  an  altitude  of  3km.  Scene  (a)  ■  horizontal  target  surface  In  sunlight.  Scene  (b)  •  verti¬ 
cal  target  surface  In  sunlight.  Scene  (c)  -  horizontal  target  surface  in  local  shadow. 
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Variation  With  Aerosol  Type 

Fig.  3-5.  Target  detection  distance  vs  LOWTRAN  aerosol  model  type  for  the  boundary  layer.  Ground  based  target  at 
sea  level  Is  viewed  from  3  km  altitude.  Scene  (a)  •  horizontal  target  surface  In  sunlight  Is  viewed  In  the  downsun 
180")  azimuthal  direction.  Scene  (b)  -  horizontal  target  surface  In  sunlight  Is  viewed  In  the  upsun  (0  —  0 “)  direc¬ 
tion.  Scene  (c)  -  horizontal  target  surface  In  local  shadow  Is  viewed  In  the  downsun  direction.  Scene  (d)  -  horizontal 
target  surface  In  local  shadow  is  viewed  In  the  upsun  direction. 
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the  average  estimated  detection  distance  with  increasing 
cloud  amount  for  target  scenes  in  local  shadow;  and  as 
stated  above,  there  is  small  net  sensitivity  to  cloud  amount 
changes  in  the  range  from  zero  to  0.8  fractional  cloud 
cover  for  sunlit  surface  targets. 


Hering,  W.S.  (1981),  “An  Operational  Technique  For 
Estimating  Visible  Spectrum  Contrast  Transmit¬ 
tance,”  SIO  Ref.  82-1,  AFGL-TR-81-0198,  ADA 
111  823,  University  of  California,  San  Diego,  Scripps 
Institution  of  Oceanography,  Visibility  Laboratory, 


VARIATION  WITH  AEROSOL  TYPE 

The  dependency  of  visual  detection  distance  on  the 
specific  optical  properties  associated  with  individual 
LOWTRAN  aerosol  models  is  illustrated  in  Fig.  3.5.  The 
model  calculations  were  made  for  a  wavelength  of  670 
nm.  For  each  of  the  4  target  scenes,  calculations  were  car¬ 
ried  out  assuming,  in  turn,  that  the  aerosol  characteristics 
of  the  boundary  layer  were  given  by  the  troposphere, 
rural,  urban,  and  maritime  aerosol  models.  All  other  fac¬ 
tors,  including  the  total  volume  extinction  coefficient  (550 
nm)  in  the  boundary  layer,  were  held  constant  and  equal 
to  the  average  reference  atmosphere  values.  A  70  percent 
relative  humidity  was  assumed  for  all  cases. 

For  the  prescribed  typical  hazy  atmosphere,  it  is  evi¬ 
dent  that  the  net  effect  of  the  differences  in  aerosol  model 
properties  on  the  calculated  detection  distance  at  670  nm 
is  small  for  each  of  the  target  scenes.  The  detection  dis¬ 
tance  is  systematically  above  average  for  the  urban  model 
and  below  average  for  the  maritime  model  atmosphere. 
This  reflects  the  differences  in  single  scattering  albedo 
which  result  in  a  higher  (lower)  path  radiance  for  the  mar¬ 
itime  (urban)  atmosphere  under  the  assumption  of  a  con¬ 
stant  total  volume  extinction  coefficient  in  the  boundary 
layer. 

Substantial  differences  are  noted  in  Fig.  3.4  between 
the  model  calculations  for  the  various  target  scenes.  For 
example,  the  visual  detection  distance  is  larger  by  a  factor 
of  4.5  for  case  (a)  where  a  sunlit  target  is  viewed  downsun 
as  compared  with  case  (d)  where  the  target  is  In  local  sha¬ 
dow  and  is  viewed  along  a  sunlit  path  In  the  upsun  direc¬ 
tion.  The  detection  distances  as  calculated  for  the  sunlit 
target  viewed  directly  upsun  (data  set  b)  and  the  shadowed 
target  viewed  downsun  (data  set  c)  are  about  the  same. 
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4.0  SUMMARY 

This  Final  Report  under  Contract  No.  F19628-82- 
C-0060  summarizes  the  results  of  several  research  studies 
directed  toward  the  development  of  models  of  the  optical 
properties  of  the  atmosphere  which  are  suitable  for  the 
operational  support  of  tactical  weapons  and  sensors  operat¬ 
ing  in  the  visible  and  Infrared  portions  of  the  electromag¬ 
netic  spectrum. 

The  overall  contractual  effort  has  been  an  integral 
part  of  the  overall  Air  Force  Geophysics  Laboratory 
(AFGL/AFSC)  program  to  develop  a  thorough  under¬ 
standing  of  atmospheric  effects  on  visible  and  Infrared 
radiation  through  the  atmosphere,  and  ultimately,  the  con¬ 
comitant  limitations  on  the  performance  of  electro-optical 
systems  such  as  target  acquisition  devices  and  precision 
guided  munitions.  , 


Cox,  C.  and  W.  Munk  (1954),  “Measurements  of  the 
Roughness  of  the  Sea  Surface  from  Photographs  of 
the  Sun’s  Glitter,”  J.  Opt .  Soc.  Am.,  44,  pp.  838-850. 

Duntley,  S.Q.,  A,R,  Bolleau,  and  R.W.  Preisendorfer 
(1957),  “Image  Transmission  by  the  Troposphere. 
I,"  /  Opt.  Soc.  Am.,41 ,  pp.  499-506. 


This  Final  Report  contains  a  brief  review  of  four 
previously  Issued  Interim  reports  which  were  submitted 
during  this  contractual  interval  In  support  of  the  primary 
model  development  task.  It  also  contains  the  complete 
documentation  for  the  FASCAT  code  and  a  discussion  of 
the  FASCAT  model  performance  under  fractional  cloud 
meteorological  conditions. 
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SUMMARY 


The  development  of  field  oriented  technique!  for  determining  the  atmoepherlc  effects  on  viiible  image  transmission 
are  deacribed  in  this  report.  Simplified  analytic  methods  are  presented  for  calculating  the  background  sky,  cloud  and  terrain 
radianoe  distribution,  and  in  turn  the  spectral  contrast  transmittance,  of  inclined  paths  of  sight  In  the  atmosphere.  The 
directional  path  radiance  it  determined  as  the  sum  of  the  singly  scattered  sunlight  component  and  a  multiply  scattered 
diltose  component  that  is  calculated  through  application  of  the  delta-Eddlngton  approximation.  Performance  testa  of  the 
model  radianoe  calculations  were  carried  out  through  comparisons  with  results  from  more  comprehensive  radiative  transfer 
models  and  by  direct  comparisons  with  hlgh-resolutlon  scanning  radiometer  measurements  gathered  by  Instrumented  air¬ 
craft  In  a  broad  range  of  environmental  conditions.  Techniques  for  the  representation  of  the  optical  parameters  required  as 
model  input  data  are  presented  and  discussed,  including  the  specification  of  cloud  optical  depth  aa  a  function  of  cloud  type. 

An  added  feature  of  the  contrast  transmittance  modeling  system  deals  with  the  combined  influences  of  atmoapherlc 
and  target  factors  in  order  to  calculate  visual  detection  range  for  objects  located  at  the  surface  and  viewed  from  aloft.  Ana¬ 
lytic  representations  of  vision  experiment  data  are  added  along  with  target  information  to  explore  the  relative  sensitivity  of 
detection  range  to  selected  changes  In  target  and  atmospheric  parameters.  Results  are  presented  for  several  examples 
assuming  simpllfled  targets  and  backgrounds. 
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Frontispiece 

The  analytic  techniques  described  in  this  report  are  a  natural  extension  to  a  broad  program  of  research  into  the  optical  properties  of  the 
atmosphere  carried  out  by  the  Visibility  Laboratory,  University  of  California,  San  Diego,  in  cooperation  with  and  under  the  sponsorship 
of  the  Air  Force  Geophysics  Laboiatory.  The  composite  photograph  shown  above,  illustrates  some  elements  of  the  field  program  and 
some  of  the  important  components  of  the  radiance  field  as  measured  by  the  instrumented  aircraft.  The  field  site  near  Mt.  Rainier,  Wash¬ 
ington,  is  one  of  many  locations  in  the  United  States  and  Europe  where  simultaneous  measurements  of  optical  and  meteorological  vari¬ 
ables  were  gathered.  Extensive  high  resolution  measurements  of  the  upper  and  lower  hemisphere  spectral  radiance  fields  for  scenes  such 
as  those  illustrated  by  the  flsheye  camera  photographs  shown  below  (flight  C-3S9)  were  of  prime  importance  both  for  algorithm  develop¬ 
ment  and  for  validation  of  the  modeling  techniques  over  a  wide  ranae  of  environmental  conditions. 
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ANALYTIC  TECHNIQUES  FOR  ESTIMATING 
IMAGE  TRANSMISSION  PROPERTIES  OF  THE  ATMOSPHERE 

Wayne  S.  Hering 


1.0  INTRODUCTION 

As  set  forth  by  Duntley  (1948),  the  degradation  of 
visible  image  contrast  along  a  viewing  path  of  increasing 
length  is  governed  by  (1)  the  attenuation  of  the  inherent 
background  and  object  radiances  by  air  molecule  and  aero¬ 
sol  particle  scattering  and  absorption,  and  (2)  the  genera¬ 
tion  of  path  radiance  by  molecular  and  aerosol  scattering 
of  incident  light  into  the  path  of  sight.  Direct  held  meas¬ 
urements  of  the  basic  optical  properties  required  for  the 
determination  of  these  factors  for  selected  path  segments 
cannot  be  made  simply  and  inexpensively.  Thus,  opera¬ 
tional  estimates  of  contrast  transmittance  must  be  derived 
from  model  simulations  that  relate  the  aerosol  directional 
scattering  and  absorption  properties  to  generally  available 
meteorological  observations  and  forecasts  (e.g.  Shettle  and 
Penn,  1979).  Additional  problems  stem  from  tho  com¬ 
plexity  of  multiple  scattering  processes  and  the  relatively 
large  amount  of  computer  time  and  capacity  required  for 
treating  the  radiative  transfer  calculations  completely  as  a 
physical  problem.  Effective  simplifications  must  be  made 
to  yield  fast  yet  consistent  calculations  of  the  radiance  dis¬ 
tribution.  The  successful  development  of  techniques  for 
real  time  estimates  of  the  atmospheric  effects  on  visible 
image  transmission  requires  a  reasonable  balance  between 
the  desired  accuracy,  the  complexity  of  the  modeling  pro¬ 
cedures,  and  the  reliability  and  representativeness  of  the 
available  input  data  Since  these  factors  may  vary  consid¬ 
erably  from  one  application  to  the  next,  it  is  important  to 
retain  reasonably  good  flexibility  and  completeness  in  the 
analytic  approximations  used  for  the  estimates. 

The  development  and  testing  of  field  oriented  tech¬ 
niques  for  calculating  directional  path  radiance  and  the 
object-background  contrast  transmittance  along  inclined 
paths  of  sight  In  clear-hazy  atmospheres  were  described  in 
a  previous  report  (Hering,  1981).  This  paper  addresses 
extensions  and  refinements  to  the  modeling  procedures 
and  their  validation.  In  particular,  techniques  are 
presented  for  (1)  estimating  the  background  radiance  and 
visible  image  transmission  characteristics  in  overcast  sky 
conditions  from  conventional  observations  or  forecasts  of 
cloud  type  and  altitude.  As  in  the  case  of  prior  studies, 
technique  development  was  based  extensively  on  experi¬ 
mental  data,  including  a  broad  series  of  simultaneous  opti¬ 
cal  and  meteorological  measurements  gathered  with  a  spe¬ 
cialty  instrumented  aircraft.  The  high  resolution  profile 
measurements  were  obtained  over  a  period  of  several 
years  in  a  wide  range  of  meteorological  and  geographical 


conditions  in  the  United  States  and  western  Europe. 

Further  extensions  of  the  modeling  concepts  deal 
with  the  combined  effects  of  both  environmental  and  tar¬ 
get  factors  on  visual  detection  range.  Although  indepen¬ 
dent  assessment  can  be  made  of  the  changes  in  contrast 
transmittance  that  are  associated  with  natural  changes  in 
the  relevant  atmospheric  variables,  the  impact  of  these 
changes  on  target  detection  distance  depends  markedly  on 
such  factors  as  inherent  target  contrast,  target  size,  and 
glimpse  time.  Analytic  representations  of  vision  experi¬ 
ment  data  (Blackwell,  1946  and  Taylor,  1964)  were  added 
to  the  contrast  transmittance  model  to  explore  the  relative 
sensitivity  of  visual  detection  distance  to  a  selected  range 
of  target  and  atmospheric  variables.  Results  are  presented 
for  several  examples  assuming  simplified  targets  and  back¬ 
grounds. 

A  review  of  the  basic  equations  and  a  summary  of 
the  analytic  approximations  used  for  calculating  path  radi¬ 
ance  and  for  estimating  visible  spectrum  contrast  transmit¬ 
tance  In  cloudless  atmospheres  are  given  in  Section  2. 
Also  presented  in  Section  2,  are  comparative  radiative 
transfer  calculations  betwoen  the  simplified  model  and 
more  complete  computational  methods.  In  Section  3, 
techniques  are  presented  for  the  representation  of  the 
input  variables  that  are  required  for  the  calculation  of 
background  radiance  and  spectral  contrast  transmittance. 
Extension  of  the  modeling  concepts  to  deal  with  image 
transmission  characteristics  for  viewing  paths  below  an 
overcast  cloud  layer  are  described  in  Section  4.  Results  of 
validation  tests  comparing  model  calculations  with  experi¬ 
mental  measurements  of  the  directional  radiance  distribu¬ 
tion  gathered  by  instrumented  aircraft  are  given  in  Sec¬ 
tion  5.  The  results  of  sensitivity  analyses  which  give  a 
measure  of  the  relative  Importance  of  individual  target  and 
environmental  factors  on  the  determination  of  target 
detection  range  are  summarized  in  Section  6. 

2.0  THEORY  AND  ANALYTIC  APPROXIMATIONS 

Although  some  refinements  and  additions  to 
cloudless-sky  modeling  procedures  are  presented,  this  sec¬ 
tion  is  primarily  a  review  of  theory  and  analytic  approxi¬ 
mations  given  In  Hering  (1981). 

Neglecting  turbulence  eflfects,  the  equation  for  the 
apparent  spectral  radiance  of  the  background  b  at  range  r 
along  the  path  of  sight  specified  by  zenith  angle  9  and 
azimuthal  angle  ^ ,  can  be  written  (Duntley  et  at.  1957) 
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,,Ltizfi4>)~  TfUfi)  bLtU„B,<t>)  +  L*U, 6,4>)  ,  (1) 

where  t L „  is  the  inherent  background  radiance  at  target 
altitude  z, , 

t 

Tr  -  exp-/  a(r)dr 


is  the  path  transmittance,  a  is  the  volume  attenuation 
coefficient,  and  L ,*  is  the  path  radiance  produced  by  the 
scattering  of  light  from  the  sun  and  from  the  surrounding 
sky,  clouds  and  terrain  into  the  path  of  sight. 

The  path  radiance  is  given  by 


L,‘(*,0, 4,)  -  /.!,(*', «, *)Tr  (z'fi)dr' , 


the  distribution  of  light  reaching  the  path.  The  first  term 
on  the  right  hand  side  of  Eq.  (3)  is  the  contribution  of  pri¬ 
mary  scattering  of  direct  solar  irradiance.  The  directional 
volume  scattering  function  may  be  expressed 

a  (r,0)  -  P(z,P)s(z)  .  (5) 

where  six)  is  the  total  volume  scattering  coefficient,  and 
P (x,0)  is  the  single  scattering  phase  function  for  com¬ 
bined  Rayleigh  plus  aerosol  particle  scattering  (see  discus¬ 
sion  in  Section  3.3),  which  defines  the  probability  that 
Incident  radiation  will  be  scattered  In  the  direction  given 
by  scattering  angle  /9.  The  phase  function  varies 
significantly  with  the  scattering  properties  of  the  aerosol 
particle  distribution  In  the  atmosphere. 

Substituting  Eq.  (5)  in  Eq.  (3),  the  expression  for 
L«Gr,0,4)  becomes 


where  L«(i',0,  <£)  is  the  path  function,  defined  as  the 
point  function  component  of  path  radiance  generated  by 
the  directional  scattering  of  light  reaching  that  point  of  the 
path.  The  expression  for  the  path  function  can  be  written 
in  terms  of  its  contributions  from  the  scattering  of  solar 
(lunar)  scalar  irradiance  ,«<z)  and  from  the  scattering  of 
sky  and  earth  radiances  L(z,0\  <f>')  as  follows  (e,g.  Gor¬ 
don  ,  1969b): 

L,(z,9,4)  -  ,« (r)cr (a./S, )  +  £  Liz  fi'&')criz  ,  (3) 

where  cr(z,/3)  is  the  directions:  volume  scattering  function 
at  angle  0  between  the  path  of  sight  and  direction  of  the 
source  light,  dtl  is  an  element  of  solid  angle,  and  /9,  is 
the  scattering  angle  with  respect  to  the  sun. 

In  turn,  the  spectral  contrast  transmittance  of  the 
path  of  sight  can  be  expressed  (see  Duntley  et  al.  1957) 
directly  as  the  product  of  the  path  transmittance,  T,(x&), 
and  the  ratio  of  the  inherent,  tL,,  and  apparent,  »£,,, 
background  radices  as  follows: 

C,M, 4)1  C0 (*,.«. 4> -  T, (*,«)* L#(z,,e,«)/tL, (i, 9, 4)(4) 

where  C,  -  i,L,  -  t,L,)l  bL,  is  the  apparent  target  con¬ 
trast  at  path  length  r,  C„  -  (,£,  -  *£»)/ is  the 
inherent  target  contrast  at  altitude  z,  and  ,4  it  the 
inherent  target  radiance.  Thus,  the  contrast  transmittance 
of  a  given  path  does  not  depend  upon  intrinsic  target 
characteristics  but  is  a  function  only  of  the  directional 
radiance  distribution  in  the  atmosphere  and  the  path 
transmittance.  The  expression  is  strictly  applicable  only 
for  monochromatic  radiation  but  may  be  applied  with  good 
approximation  to  reasonably  broad  spectral  bands  in  the 
visible  portion  of  the  spectrum. 

From  Eq.  (3)  we  note  that  the  path  function  and  in 
turn  the  contrast  transmittance  as  calculated  from  Eq.  (4) 
depend  upon  the  direction  of  the  viewing  path  relative  to 


L$<z,6,+)“  ,*(z)P(z,fi)s(z) 


+  f£Hz,0\'t>,)P^z,p)s^z)dn, 


Modelling  techniques  for  the  calculation  of  the  diffuse 
component  of  the  path  function  are  discussed  in  the  fol¬ 
lowing  paragraphs. 

2.1  Calculation  of  the  Path  Function 

The  second  term  on  the  right  hand  side  of  Eq.  (6)  Is 
the  component  of  the  path  function  resulting  from  the 
scattering  of  diffuse  radiance  reaching  the  path  from  the 
surrounding  sky  and  terrain.  It  has  a  directional  depen¬ 
dence  due  to  the  asymmetry  In  thq  background  sky  and 
earth  radiance  distribution  and  the  scattering  phase  func¬ 
tion  asymmetry.  Precise  numerical  calculation  of  the  path 
radiance  resulting  from  the  complex  multiple  scattering 
processes  requires  large  amounts  of  computer  time.  For 
this  reason,  rapid  approximate  methods  are  employed 
extensively  for  radiation  transfer  calculations.  The 
appropriate  choice  of  computational  method  from  among 
the  variety  of  available  methods  depends  upon  the  results 
desired  for  the  application  at  hand. 

While  it  is  important  to  retain  complete  directional¬ 
ity  for  calculation  of  the  path  radiance  component  due  to 
single  scattering  of  direct  solar  radiance,  approximate 
hemispherical  two  stream  methods  can  be  used  effectively 
for  fast  calculation  of  the  path  radiance  component  due  to 
scattering  of  the  background  sky  and  terrain  radiances  at 
any  point  and  direction  provided  that  the  asymmetric 
influence  of  the  prominent  forward  scatter  peak  and  the 
irradiance  profile  are  managed  adequately.  The  delta- 
Eddington  approximation  introduced  by  Joseph,  Wiacombe 
and  Weinman  (1976)  satisfies  the  requirement.  It  extends 
the  standard  Eddington  approximation,  which  assumes  a 
simple  cosine  dependence  of  the  single  scattering  phase 
function,  through  approximation  of  the  phase  function, 
by  a  truncated  forward  scatter  peak  and  a  two-term 
phase  function  expansion, 


APPENDIX  A 


An  P„  fc)  -  V'A  (1  -  cos/3)  +  (1  -  /'Ml  +  3f  'cos/8)  (7) 


where/'  is  the  fractional  scattering  represented  by  the  for¬ 
ward  peak  and  s'  is  the  asymmetry  factor  of  the  truncated 
phase  function.  In  effect,  the  delta-Eddington  approxima¬ 
tion  transforms  most  of  the  enhanced  radiance  in  the  solar 
aureole  into  the  direct  solar  flux  component,  and  it  is  gen¬ 
erally  assumed  that 


/'(*  J-s’tx). 


This  assumption  is  commensurate  with  representation  of 
the  actual  phase  function  with  a  single  term  Henyey- 
Greenstein  phase  function.  Alternate  expressions  are 
required  (see  McKellar  and  Box,  1981)  for  other  phase 
(Unction  representations  such  as  the  2-term  Henyey- 
Greenstein  functions.  Joseph  etal  (1976)  show  that  cal¬ 
culations  of  radiative  transfer  with  the  delta-Eddington 
approximation  can  be  carried  out  with  the  standard 
Eddington  computer  code  (Shettle  and  Weinman,  1970) 
with  the  following  changes  of  variable 


Ar'-  (!-«/')  Ar  , 


«'- «(!-/')/  (1-ui/'), 


(1-/0, 


where  Ar  is  the  optical  thickness  of  the  layer  and  at  is  the 
single  scattering  albedo. 

As  an  integral  part  of  the  technique  for  estimating 
directional  contrast  transmittance,  the  products  of  the 
Eddington  computer  program  (Shettle  and  Weinman, 
1970)  are  used  directly  to  calculate  the  second  term  on  the 
right  hand  side  of  Eq.  (6). 

From  the  standard  Eddington  approximation,  the 
dlffUse  radiance  is  assumed  to  be  given  by 


£,(z,0\  A')  “  Ld(i)  +  coa0' . 


As  shown  by  Shettle  (1981),  if  we  substitute  Eq.  (12)  and 
the  delta-Eddington  approximations  given  by  Eq.  (7),  (8) 
and  (11)  into  the  last  term  of  Eq.  (6)  and  Integrate  over  O' 
and  <f>\  we  have 


JJT  <A')/*(z,/3)s (z)Jn 


-  i(i)  l Ld(j)  +  jlfl  CtlcosO)  . 


A  recent  refinement  In  the  above  expression  for  the 
diffuse  component  of  the  path  function  consists  of  an 


additional  term  to  help  account  for  the  effects  of  the 
azimuthal  asymmetry  in  the  background  sky-terrain  radi¬ 
ance  distribution.  The  term  is  introduced  as  a  rough 
approximation  to  an  iteration  of  the  singly  scattered  sun¬ 
light  component  of  the  path  (Unction.  Its  magnitude  is 
assumed  proportional  to  the  single  scattering  component, 
,« (r)F(z,/9),  and  proportional  to  the  fractional  contribution 
of  the  total  diffuse  energy,  An  LD(jt),  to  the  sum  of  the 
total  diffuse  energy  and  singly  scattered  sunlight  com¬ 
ponent  at  the  level  of  computation.  In  other  respects,  the 
modification  conforms  with  the  extension  of  the  delta- 
Eddington  model  introduced  by  Davies  (1980);  in  particu¬ 
lar,  the  contribution  of  the  asymmetry  term  is  proportional 
to  aind,  slno,  and  cosd.  The  modified  expression  for  the 
diffuse  component  of  the  path  function  is 


L(x,o\  *')P(j (*)dn 


-  s(r)  Up  Or)  +  f£0.(i)cos0  +  si,fi»Of)cos0)  ,  (14) 


t0»(z)  ■ 


An  LD(x)P(x>p)  ,*b)a\n9,i\n0 


Commensurate  with  the  delta-Eddington  approximations 
and  Eq.  (9),  the  expression  for  the  solar  scalar  Irradiance 
at  altitude  z  is 


,«(*)- l«(oo)exp(-Tx/eoi0l) , 


where  ,«(<»)  is  the  extraterrestrial  solar  scalar  irradiance. 
The  optical  depth  rx  is  assumed  equal  to  the  delta- 
Eddington  optical  depth,  r',  everywhere  except  within  the 
forward  peak  03  <  25  deg)  where  tx  is  equal  to  the 
unmodified  value  of  r.  Substituting  Eqs.  (13)  and  (14) 
into  Eq.  (6),  the  expression  for  the  path  function  becomes 


i,,Ot,0,(p)  -  j(z)  (F(z,0,),«(<»)exp(-rx/coi0l) 


+  L0(z)  +  gZ.fl(z)cose  +  #L0«.U)coi^) . 


2.2  Calculation  of  Sky  and  Terrain  Radiance 

For  an  assumed  plane  parallel  and  horizontclly 
homogeneous  atmosphere,  the  directional  path  radiance, 
l/(z,0,  d>),  and  the  radiance,  l(z,0,  </>),  can  now  be  cal¬ 
culated  from  Eqs,  (1),  (2)  and  (17)  through  finite  summa¬ 
tion  over  adlacent  atmospheric  layers  using  the  trapezoidal 
rule.  For  upward  paths  of  sight,  the  Inherent  background 
radiance  of  clear  sky  at  target  altitude  is  given  by 


4  to  (r,  .0,  0)  •"  Lm  (z,  ,8 , 0)  , 


where  Ll  (z,,0,6)  is  the  path  radiance  as  determined  at 
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altitude  z,  for  the  slant  path  from  the  top  of  the  atmo* 
sphere  to  z, .  Looking  downward,  the  inherent  background 
radiance  for  the  case  of  uniform  Lambertian  reflectance 
and  a  horizontal  surface  is  given  by 

+  L*(r,  (19) 

where  R(0,<l>)  is  the  local  surface  reflectance, 
£(0,rf)  -  «,  (O)cos0,  +  nlLD  (0)  +  V 3/.o  (0)l  is  the 
downwelling  irradiance  at  the  surface,  and  T„  is  the 
transmittance  of  the  slant  path  from  the  surface  to  target 
altitude  z, . 

2.3  Comparisons  with  Other  Model  Calculations 

Tests  were  carried  out  to  evaluate  the  performance 
of  the  simplified  model  for  the  radiance  calculations  as 
described  above.  Model  calculations  of  radiance  distribu¬ 
tions  were  compared  with  the  results  calculated  from  more 
comprehensive  radiative  transfer  models,  which  are 
mathematically  precise  but  in  general  require  relatively 
large  computer  capacity  and  time  for  the  radiative  transfer 
calculations.  For  these  tests,  the  comparative  accuracy 


was  determined  under  the  assumption  that  the  physical 
system  is  perfectly  observed,  and  the  input  data  are  identi¬ 
cal  for  the  approximate  operational  model  and  the 
comprehensive  numerical  models.  The  set  of  examples 
used  for  the  analysis  were  selected  from  those  proposed  by 
the  ad  hoc  Working  Group  of  the  Radiation  Commission, 
International  Association  of  Meteorology  and  Atmospheric 

Table  1.  Summary  of  numerical  comparison  cues  listed 
in  the  Working  Group  Report  (Lenoble,  1977)  of  the 
Radiation  Commission,  IAMAP,  for  5  reference  atmo¬ 
spheres. 


Case 

Atmosphere 

Optical 

Thickness 

Single. 

Scattering 

Albedo 

Solar 

Zenith 

Angle 

1 

Hus  L 

1 

1.0 

0* 

2 

Hus  L 

1 

0.9 

0* 

3A 

Haul  (6-0*) 

1 

0.9 

60* 

3B 

HUS  L(6-  90‘) 

1 

0.9 

60* 

JC 

Hue  L(4*  lW*) 

1 

0.9 

60* 

4 

Cloud 

64 

1.0 

0* 

s 

Cloud 

64 

0.9 

'  0* 

to0 


10° 


10-'  + 


lO-'-f  10-J 

CASE  I 


to' 


10"*-L 

CASE  3 


10° 


to-'-f 


CASE3 

a-o* 


io-i 


10-J-j-  IO-'t 


CASE) 

a-w 


to-1 


10-J-t- 

CASE  3 
A-  ISO* 


Soltr 
Zenith 
Anglo  -  0* 


Solir 
Zenith 
Angle  -  60* 


Fig.  1.  Comparative  radiance  calculations  for  optical  depth  0.1  for  reference  esses  1,2  and  3.  Model  results  (dashed  lines)  sre  compared 
with  the  spherical  harmonics  method,  (Lenoble,  1977),  (solid  lines).  Optical  parameter  data  for  each  case  are  given  in  Table  1. 
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Phyiics  (Lenoble,  1977,  see  Table  1).  The  results  of  the 
spherical  harmonic  method  listed  in  the  Working  Group 
report  were  chosen  for  a  reference  standard  since  com* 
parative  calculations  using  this  method  were  presented  for 
all  live  test  cases.  A  single  layer,  plane  parallel  and  homo¬ 
geneous  atmosphere  with  a  black  underlying  surface  was 
prescribed  in  each  case.  Three  cases  have  a  Haze  L  (Deir- 
mendjian,  1969)  aerosol  particle  distribution  with  an  opti¬ 
cal  thickness  of  1,  and  the  other  two  have  assumed  cloud 
particle  distributions  with  an  optical  thickness  of  64.  The 
phase  function  for  single  scattering  and  the  single  scatter¬ 
ing  albedo  are  specified  for  each  case.  Comparative  results 
of  the  radiance  calculations  are  illustrated  in  Figs.  1-4. 

It  is  important  to  note  that  the  agreement  among 
the  more  rigorous  numerical  methods,  as  given  in  the 
Working  Group  report,  is  generally  within  a  few  percent 
depending  upon  the  completeness  of  scattering  history 
retained  in  the  calculations.  We  see  from  Fig.  1  thnt  the 
departures  of  the  simplified  model  radiance  calculations 
from  the  reference  values  for  the  haze  cases  are 
significantly  larger,  yet  the  indicated  accuracy  is  good  when 
considered  in  context  with  the  uncertainties  and  complete¬ 
ness  of  the  observational  data  base  from  which  the  esti¬ 
mates  must  be  derived  in  operational  practice.  The 


simplified  model  requires  only  a  few  seconds  of  large  scale 
computer  time  for  each  case,  and  the  model  is  designed 
for  easy  microprocessor  application  for  real  time  require¬ 
ments.  The  approximations  inherent  In  the  computation¬ 
ally  fast  model  serve  to  smooth  out  the  asymmetry  in  the 
calculated  radiance  component  due  to  the  scattering  of 
incident  diffuse  light;  hence  it  is  expected  that  the  error  of 
estimate  would  continue  to  increase  somewhat  in  the  solar 
azimuthal  plane  (<f>  -  0-180  deg)  as  the  solar  zenith  angle 
increases  beyond  60”. 

The  relative  effectiveness  of  the  azimuthal  asym¬ 
metry  term,  g  La* cos0,  is  Illustrated  in  Table  2.  Model 
calculations  for  Case  3  (6,-60  deg,  <w-0.9)  are  listed 
with  and  without  the  azimuthal  component  contribution 
for  comparison  with  the  reference  model  results.  It  is  evi¬ 
dent  that  the  accuracy  of  Eq.  (17)  is  much  better  with  the 
azimuthal  correction.  However,  systematic  errors  remain 
such  that  further  study  of  this  aspect  of  the  modeling  pro¬ 
cedures  is  particularly  important. 

The  relative  performance  of  the  model  calculations 
of  radiative  transfer  in  a  high  extinction  environment  is 
illustrated  in  Figs.  3  and  4.  Both  of  the  cloud  examples  in 
the  reference  series  assume  an  overhead  sun  position. 
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Ft|.  2.  Syne  as  Fig.  1  except  for  0.7S  optical  depth. 
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TaUt  2.  Comparison  of  Cue  3  radiance  distribution!  u  calculated  by  the  Spherical  Harmonics  method  (column  SH),  the  simplified  model 
without  azimuthal  uymmetry  term  (column  A),  and  eimpllfled  model  with  the  uimulhal  uymmetry  term  (column  B)  The  tolar  zenith 
«n|le  ii  60  deg,  Calculation!  are  ahown  for  a  level  near  the  top  of  the  hue  layer  (r  -  0.1)  and  a  level  near  the  base  (r  -  0.7S).  Note  that 
the  azimuthal  term,  t  £/>»  coed,  I*  zero  for  6"  90  deg,  hence  columns  A  and  B  are  identical* 


T-0.1 


e-9o* 

4-0* 

d-lto* 

SH 

A 

SH 

A 

B 

SH 

A  ' 

• 

1 

1.01-2 

9.31-3 

1.01-2 

931-3 

9.31-3 

1.01-2 

9.31-3 

9.31-3 

9,30-3 

9.43-3 

1,63-t 

1,66-1 

1.12-2 

3.30-3 

1.001 

2,79-3 

■  W 
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7,71-1 

7,50-1 

7.14-2 

4.30-1 
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1.38-2 

1.17-0 

1.11-0 

1.17-0 

7.46-3 

9.14-3 

6.39-3 

2.73-2 
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9.07-1 

7,90-1 

1.19-1 
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2.00-2 

1.43-2 

•.2 

9.11-2 

7.72-2 

6,30-1 

3,10-1 

5.30-1 

6.92-2 

t.SI-2 

6,47-2 

■i4 

6.36-2 

3,92-2 

3.01-1 

1,59-1 

2.21-1 

5.11-2 

6.40-2 

S.IO-2 

-.6 

4.32-2 

4,36-2 

1.36-1 

7.32-2 

9.71-2 

3.72-2 

4.14-2 

1.92-2 

2.17-2 

3,24-2 

3.74-2 

4,04-2 

4.73-2 

2.13-2 

3.S4-2 

3,21-2 

•1,0 

2,00-2 

2.47-2 

2.00-2 

2.47-2 

2.47-2 

2.00-2 

2.47-2 

2,47-1 
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SH 

A 

SH 

A 

B 

SH 

A 

• 

1 
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7.64-2 

6.73-2 

7.64-2 

7.66-2 
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7*4-1 

7,64-2 

.1 

6.69-2 

7.53-2 

6.79-1 

5.25-1 

6.10-1 

2.61-2 

3.16-2 

4,71-2 

.6 

7.31-2 

7.72-2 

2.10-0 

2.0IMI 

2.21-0 

3.002 

5.77-1 

5,19-2 

.4 

1.12-2 

1.44-2 

2.71-0 

2,39-0 

2.61-0 

4.13-2 

6,97-2 

6.13-2 

.2 

1.03-1 

9.44-2 

1.61-0 

1.03-0 

1.300 

5.56-2 

1.99-2 

7,34-2 

4,41-2 

4,36-2 

2.66-1 

1.66-1 

2.39-1 

2,92-2 

4.61-1 

141-2 

1.94-2 

2.03-2 

1.41-2 

5.12-2 

7.202 

1.41-2 

2.11-2 

1.77-2 

1.00-2 

1.02-2 

2.93-2 

1.19-2 

2,53-2 

1.96-3 

I.IS-2 

1,97-1 

3.91-2 

4.67-2 

1.10-2 

6.71-3 

1.603 

6.24-3 

6.23-1 

4,57-3 

-1.0 

4.02-3 

2.14-3 

4.02-3 

2.14-3 

2.14-3 

4.02-3 

2.14-1 

214-1 

I  3.  S'  ne  u  Fig.  1  except  for  cloud  case  4  at  optical  depths  Fig.  4.  Same  u  Fig.  3  except  for  cloud  case  S  with  single  scatter- 

3  i,  12.8,  and  48  and  no  absorption.  ing  albedo  of  0.9, 
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The  agreement  with  the  reference  model  calculations 
holds  well  for  both  cloud  cases,  although  the  disparity 
increases  with  increasing  depth  for  the  case  with  strong 
absorption.  The  simplified  model  calculations  of  in-cloud 
radiance  were  carried  out  through  division  into  sub-layers 
with  boundaries  corresponding  to  levels  where  reference 
model  calculations  were  available.  These  subdivisions 
represent  large  increments  of  optical  depth  in  the  cloud 
cases.  However,  the  estimates  if  in-cloud  radiance  are 
good  in  general  over  the  full  range  of  optical  depth  and 
observing  angles  shown  in  Fig.  3  and  4,  A  notable  excep¬ 
tion  is  for  the  data  point  directly  upsun,  0-0,-iHO,  at 
optical  depth  3.2.  The  calculated  values  of  in-cloud  radi¬ 
ance  for  this  data  point  result  in  overestimates  by  at  least 
an  order  of  magnitude  and  plot  offscale  in  Figs.  3  and  4. 
The  modeling  procedures  require  important  refinement  in 
layer  by  layer  resolution  in  order  to  handle  combined 
effects  of  the  rapid  attenuation  of  the  direct  solar  beam  in 
the  cloud  near  the  upper  limit  of  the  cloud  layer  and  the 
strong  change  in  the  single  scattering  phase  function  with 
scattering  angle  for  /3  less  than  3  deg,  The  large  errors  are 
confined  to  viewing  paths  within  a  few  degrees  of  the  sun 
and  to  the  region  near  upper  limit  of  the  cloud  layer,  In 
this  region,  the  singly  scattered  sunlight  component  makes 
a  significant  contribution  to  the  path  function  and  it 
changes  rapidly  along  the  path. 

3.0  APPLICATION  OF  ANALYTIC  TECHNIQUES 

Environmental  data  and  forecasts  relevant  to  the 
determination  of  the  background  spectral  radiance  distri¬ 
bution  and  image  transmission  characteristics  may  be 
available  in  many  forms  and  with  varying  degrees  of  com¬ 
pleteness  and  representativeness,  Emphasis  has  been 
placed  on  the  development  of  a  general  modular  format 
for  requisite  data  entry  seeking  to  take  full  advantage  of 
all  available  Information  important  for  a  particular  applica¬ 
tion  of  the  radiance  model.  A  summary  of  input  variables 
and  alternate  techniques  for  the  specification  of  the  atmos¬ 
pheric  properties  are  Included  in  Appendix  A.  Basic  data 
entries  include  the  wavelength  representative  of  the  sensor 
spectral  characteristics,  solar  zenith  angle,  and  the  obser¬ 
vational  paths  of  interest  for  the  problem  at  hand. 
Specification  of  the  underlying  surface  reflectivity  is  dis¬ 
cussed  below  In  Section  3.1. 

Several  options  are  available  for  the  Input  of  specific 
atmospheric  variables.  One  may  introduce  as  many 
atmospheric  layers  as  warranted  by  the  accuracy  and  com¬ 
pleteness  of  the  meteorological  observations  or  forecasts 
available  for  the  specification  of  the  optical  parameter 
profiles.  Minimum  information  for  each  designated  layer 
i  icludes  the  altitude  limits,  the  average  scattering 
coefficient,  the  average  absorption  coefficient  and  the  sin¬ 
gle  scattering  phase  function  for  aerosol  particle  scattering. 
Techniques  for  the  representation  of  these  variables  are 
discussed  in  Sections  3.2  and  3.3. 

3.1  Representation  of  Surface  Reflectance 

The  background  surface  reflectance  contributes  to 
the  calculated  radiance  distribution  in  two  ways,  First,  the 


general  or  broad  area  average  surface  reflectance  is 
entered  to  calculate  the  component  of  path  radiance  that  is 
generated  by  light  reaching  the  path  through  scattering 
from  the  underlying  surface  and  in  turn  scattered  in  the 
direction  of  the  sensor  as  given  by  Eq.  (6).  For  a  land 
surface,  the  reflectance  is  assumed  to  conform  with 
Lambert's  law,  so  that  the  resultant  radiance  is  indepen¬ 
dent  of  observation  angle  and  depends  only  on  the 
downwelling  irradlance  and  the  average  surface 
reflectance. 

Second,  the  local  background  reflectance,  which  may 
differ  from  the  area  averaged  reflectance  discussed  above, 
is  entered  for  the  determination  of  the  Inherent  back¬ 
ground  radiance  and  its  contribution  to  the  apparent  spec¬ 
tral  radiance  as  given  by  Eqs.  (I)  and  (19). 

The  above  algorithms  deal  with  surface  reflectance 
over  land  areas.  Another  option,  for  use  primarily  over 
an  open  water  surface,  assumes  specular  surface  reflection. 
Here  the  surface  radiance  is  a  function  of  both  the  viewing 
angle  and  the  downwelling  radiance  distribution.  For  the 
water  option,  the  Fresnel  reflectance  of  the  water  surface 
as  a  function  of  observation  angle  is  set  equal  to  average 
values  for  a  sea  surface  roughened  by  light  surface  winds 
(4m  see-1).  These  values  as  calculated  by  Gordon 
(1969a),  assume  that  the  crosswind  and  upwind  wave 
slope  probabilities  of  Cox  and  Munk  (1934)  may  be 
approximated  by  a  single  circular  distribution  of  wave 
slopes. 

3.2  Representation  of  Total  Volume 
Scattering  Coefficient 

The  spatial  distribution  of  total  volume  scattering 
coefficient  deserves  prime  consideration  since  it  is  the 
major  determinant  of  visible  spectrum  contrast  transmit¬ 
tance.  Techniques  for  ther  specification  of  the  scattering 
coefficient  profile  have  been  Investipted  during  the  course 
of  the  aircraft  measurement  and  analysis  program.  Some 
of  the  results  of  these  studies  and  their  application  to 
operational  modeling  procedures  were  summarized  by 
Herlng  (1981).  A  brief  review  is  given  here  for  immedi¬ 
ate  reference. 

For  profile  modeling  purposes,  it  is  important  to 
consider  a  conservative  measure  of  scattering  coefficient 
that  in  the  absence  of  local  aerosol  particle  sources  or 
sinks  does  not  change  appreciably  following  the  air 
motion.  The  optical  scattering  mixing  ratio,  Q Or),  is  such 
a  parameter,  As  the  vertical  mixing  within  an  Identifiable 
atmospheric  layer  becomes  more  complete,  Q(z )  becomes 
more  constant  with  height  within  the  layer. 

The  optical  scattering  ratio  Is  defined 

C?U)  -  s(i)/r*  U) ,  (20) 

where  sR(z)  is  the  total  volume  coefficient  for  Rayleigh 
scattering  at  altitude  z.  Note  that 

* (r)  -  I*  (z)  +  sM(z)  -  ut(z)a(z)  ,  (21) 
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where  a  (z)  is  the  extinction  coefficient  and  s«(z)  is  the 
aerosol  scattering  coefficient.  It  follows  that  the  aerosol 
scattering  ratio  is  given  by 

%(*)/*«(*)-  Q<*>-  l.  (22) 

The  aerosol  scattering  ratio  also  would  be  constant  under 
conditions  of  complete  aerosol  mixing.  An  additional 
computational  advantage  of  a  scattering  ratio  representa¬ 
tion  is  that  it  provides  normalization  with  respect  to  both 
density  altitude  and  wavelength. 

Profiles  of  Q(x)  derived  from  the  extensive  series 
of  airborne  optical  measurements  made  by  the  Visibility 
Laboratory,  reveal  large  variability  depending  upon  the 
aerosol  particle  source  strength  and  the  nature  of  the  con¬ 
vective  and  turbulent  mixing  processes.  The  problem  Is  to 
model  the  essential  characteristics  of  the  Q(z)  profiles  In  a 
way  that  recognizes  operational  observing  and  forecasting 
limitations  yet  takes  maximum  advantage  of  existing  capa¬ 
bilities.  A  prominent  feature  of  the  daytime  aircraft 
soundings  over  inland  areas  was  the  marked  tendency  for 
Q(z)  to  remain  essentially  constant  with  height  within  the 
boundary  layer,  and  also  in  the  relatively  haze-free  region 
of  the  upper  troposphere  above  the  primary  haze  layer.  It 
should  be  emphasized  that  the  assumption  of  constant 
scattering  ratio  with  height  does  not  hold  well  for  ground- 
based  stable  layers  with  little  vertical  mixing  such  as  those 
associated  with  the  nocturnal  formation  of  fog.  However, 
for  application  to  problems  of  contrast  transmittance  in 
hazy  atmospheres  In  the  daytime  following  the  dispersion 
of  any  surface  inversion  existing  at  sunrise,  a  simple  3- 
layer  troposphere  model  with  constant  Q(z)  In  each  layer 
provides  In  most  cases  a  good  first  approximation  of  the 
aerosol  scattering  profile.  Thus,  the  forecasting  problem  is 
reduced  to  the  prediction  of  the  upper  altitude  limit  of 
mixed  boundary  layer  and  the  average  scattering  ratio 


within  each  layer.  An  evaluation  of  model  representations 
of  the  high  resolution  scattering  ratio  profiles  as  measured 
during  several  deployments  of  the  instrumented  aircraft  in 
western  Europe  is  given  in  Hering  (1981). 

Three  examples  of  scattering  coefficient  profiles  as 
measured  by  the  integrating  nephetometer  are  shown  in 
Fig.  S.  Each  experimental  flight  usually  included  profile 
measurements  with  4  spectral  filters  having  peak 
wavelengths  of  475,  550,  660,  and  750  nm.  The  4  profiles 
were  measured  sequentially  over  approximately  a  2-hour 
period.  Data  for  individual  profiles  were  extracted  at  30-m 
intervals.  Figure  5  illustrates  some  recurrent  features  of 
the  lower  tropospheric  profiles;  the  well  mixed  boundary 
layer  In  late  morning  and  afternoon  over  inland  areas  of 
northern  Europe  (Soesterberg  profile),  the  multi- 
structured  profile  in  low-level  fog  conditions  offshore  near 
Rodby,  Denmark,  and  the  deep  haze  layer  typical  of  the 
summer  profiles  measured  in  the  Mediterranean  Sea  area 
(Trapani,  Sicily). 

Cursory  analysis  of  a  large  number  of  simultaneous 
aircraft  measurements  of  temperature  and  scattering 
coefficient  profile  structure  indicate  that  a  temperature 
lapse  rate  equal  to  the  standard  lapse  rate  of  6.5*C  per  km 
or  greater  will  support  rather  complete  vertical  mixing  of 
aerosol  particles  in  the  size  range  contributing'  predpm- 
inantly  to  extinction  In  the  visible  portion  of  the  spectrum. 
On  the  other  hand,  atmospheric  layers  of  thickness  greater 
than  about  300  m  and  lapse  rate  less  than  about  4.3*C  per 
km  significantly  Inhibit  vertical  mixing  and  in  general 
appear  sufficient  to  cap  the  mixed  boundary  layer. 

3.3  Representation  ef  Single  Scattering 
Phase  Function  and  Single  Scattering  Albedo 

The  single  scattering  phase  function  Piffi)  as 
employed  in  Eq.  (6)  Is  a  bulk  parameter  of  the  atmos- 


OPTICAL  SCATTERING  RATIO,  Q(Z)  OPTICAL  SCATTERING  RATIO,  Q(Z)  OPTICAL  SC  A!  TERINO  RATIO,  Q(Z) 

FIs.  5.  Scattering  ratio  profiles  as  measured  for  flight  C-372  (Soesterberg,  Netherlands),  flight  C-390  (Rodby,  Denmark)  and  flight  C-462 
(Trapani,  Sicily).  The  heavy  dashed  line  Is  an  objective  model  repreaentatlon  of  the  profile.  The  spectral  nephelometer  measurements 
were  corrected  for  the  effects  of  Incomplete  purging  of  the  Instrument. 
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pheric  layer,  representing  the  combined  aerosol  and  Ray* 
leigh  phase  functions.  It  is  given  by 

MM)- 08) 

where  Pu((l,i)  is  the  phase  function  for  aerosol  scattering 
and  the  theoretical  Rayleigh  phase  function  is 


So  that  the  fraction  of  radiation  absorbed  for  each  single 
photon  collision  Is  1  -  w(z).  Commensurate  with  other 
variables,  »(r)  Is  introduced  as  a  bulk  parameter  and  is  a 
measure  of  the  combined  absorption  effects  of  both  air 
molecules  and  aerosol  particles.  In  the  visible  spectrum, 
w(z)  in  the  boundary  layer  ranges  from  1.0  (no  absorp¬ 
tion)  to  less  than  0.7  in  urban  atmospheres  with  appreci¬ 
able  carbon  concentration. 


PR(p)~  3(1  + cos*#  /  16* .  (24) 

The  phase  function  for  single  scattering  has  the  normal¬ 
ized  form, 

£/>(*,»,  *)rfn  -  1  (25) 

The  computational  scheme  offers  several  options  for 
specification  of  the  average  single  scattering  phase  func¬ 
tion  Ptofl),  for  each  designated  atmospheric  layer. 
Detailed  estimates  of  Py(fi)  available  from  prescribed 
aerosol  models  such  as  those  associated  with  LOWTRAN 
(Shettle  and  Penn,  1979),  may  be  entered  directly  and 
used  in  a  table  look-up  format.  A  second  option  is  the 
representation  of  the  aerosol  phase  functions  by  two-term 
Henyey-Greenstein  functions  (Irvine,  1968)  as  follows, 

Puffi<ii\iisc) m  c PnoVl,i\)  +  (l  -  c)/Vo(0,f]) ,  (26) 


Pm(fi, t) -  0  -  t1)/  ( 4w (1  —  2rcoi0  +  ,  (27) 

and  the  asymmetry  factor,  g ,  is  given  by 


P(fi)oo»p  sing  d  p 


To  the  extent  that  the  Henyey-Greenstein  asymmetry  fac¬ 
tors,  fi  and  g],  and  the  partitioning  factor,  c,  can  be 
estimated  from  reference  atmosphere  calculations  they  can 
be  entered  for  the  specification  of  the  required  scattering 
phase  functions.  In  the  absence  of  information  that  may 
serve  to  identify  the  directional  scattering  properties  of 
individual  aerosol  layers,  yet  another  option  may  be  used 
as  described  in  Hering  (1981).  Empirical  functions  were 
developed  which  prescribe  the  Heny<  y-Oreensttin  func¬ 
tion  parameters  fi,  gj,  and  c,  as  a  function  of  the  scatter¬ 
ing  ratio  (?(*>.  The  derived  expressions  are  based  to  a 
large  extent  on  the  average  phase  functions  for  selected 
ranges  of  scattering  coefficient  as  measured  by  liarteneva 
(I960). 

The  other  direct-entry  variable  for  each  designated 
atmospheric  layer  is  the  single  scattering  ^bedo,  wG).  It 
la  defined  by  the  expression 

«(i )  —  n(t)/a(j) .  (29) 


4.0  APPLICATION  TO  OVERCAST  SKY  CONDITION 

The  clear-sky  modeling  concepts  and  computational 
scheme  can  be  applied  in  the  same  form  to  calculate  the 
visible  spectrum  radiance  fields  and  contrast  transmittance 
for  continuous  cloud  cover  conditions.  Given  reliable 
empirical  estimates  of  cloud  optical  properties,  the  delta- 
Eddington  method  yields  good  estimates  of  the  Irradiance 
profile!  in  overcast  atmospheres.  In  turn,  the  directional 
path  radiance  fields  and  the  contrast  transmittance  for 
selected  slant  paths  may  be  calculated  from  the  expres¬ 
sion!  given  in  Section  2.  In  practice,  one  muat  rely  on 
cloud  model  estimates  of  the  optical  depth,  the  single 
scattering  phase  function,  and  the  single  scattering  albedo 
as  a  function  of  standard  cloud  type. 


4.1  Eitlmatee  of  Optical  Depth  from 
Conventional  Cloud  Observations 

Rough  empirical  approximations  of  the  average  opti¬ 
cal  depth  as  a  function  of  cloud  type  were  derived  through 
referonoe  to  a  climatological  summary  of  the  SOLMET 
(1977)  data  base  u  included  In  a  report  by  Shapiro 
(1982).  The  SOLMET  data  base  wu  oomplled  by  the 
National  Oimatlc  Center  under  the  sponsorship  of  the 
Department  of  Energy.  It  consists  of  a  long  record  of 
edited  and  adjusted  tolar  radiation  data  (hourly  averagoa) 
for  26  National  Weather  Servioe  stations  spread 
throughout  the  United  State*  coupled  with  the  conven¬ 
tional  hourly  surface  observations  for  these  stations,  The 
climatological  tabulations  of  the  broad-band  nolar  Irradi¬ 
ance  were  available  in  the  form  of  averages  for  all  stations 
and  all  seasons  for  standard  cloud  types  (10/10  cloud 
cover)  for  selected  solar  zenith  angles.  These  data  were 
combined  with  the  average  irradiance  data  for  cloudless 
skies  to  calculate  the  ratio  of  the  average  broad-band  Irra¬ 
diance  for  the  sky  covered  by  clouds  of  specific  types  to 
the  average  broad-band  irradiance  for  cloudless  atmo¬ 
spheres. 

It  a  important  to  note  that  the  calculated 
overcast/clear  irradiance  ratios  for  different  cloud  types 
and  solar  angles  derived  from  the  SOLMET  data  base 
agree  very  closely  with  the  average  overcast/clear  irradi¬ 
ance  ratios  calculated  by  (Haurwltz  (1948))  for  compur- 
abie  cloud  classifications  using  8  years  of  observations  at 
the  Blue  Hill  Observatory  in  Massachusetts. 

To  the  extent  the  fractional  transmittance  for  the 
different  cloud  types  as  determined  from  the  broad-band 
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(.293-2.ty**>)  pyronometer  measurements  ire  repreaenta- 
tive  In  general  of  the  viiible  spectrum,  one  can  um  the 
delta-Eddiniton  method  to  estimate  the  effective  average 
optical  depth  associated  with  the  individual  cloud  types. 
The  results  of  trial  calculations  made  to  obtain  rough  esti¬ 
mates  of  the  optical  depth  of  standard  cloud  types  are 
shown  in  Pig.  6.  Through  successive  approximation, 
determination  was  made  of  the  best  least-squares  fit 
between  the  average  overcast/clear  irradiance  ratios  as  a 
function  of  solar  zenith  angle  derived  from  the  SOLMET 
data  base  and  the  corresponding  values  calculated  with  the 
delta-Eddlngton  model.  A  central  wavelength  of  SSO  nm 
and  an  average  surface  reflectance  of  .07  were  used  for  the 
model  calculations  and  absorption  within  the  cloud  layer 
was  assumed  negligible.  The  altitude  limits  of  the  cloud 
layers  were  chosen  so  as  to  conform  in  general  with  typical 
values  for  each  cloud  type. 

Except  as  noted  below,  the  asymmetry  factor,  y, 
was  assumed  equal  to  .85  in  calculations  for  all  cloud 
types.  It  should  be  emphasized  that  the  model  calcula¬ 
tions  are  sensitive  to  variations  in  g,  However,  in  the 
.  case  of  water  clouds,  the  variation  in  g  is  small  (Twomey, 
1980)  and  the  resultant  uncertainties  are  negligible  in 
comparison  with  the  natural  variations  In  optical  depth  for 
each  cloud  type.  In  the  case  of  Ice  clouds,  very  significant 
fluctuations  in  g  are  associated  with  variations  in  ice  cry¬ 
stal  size,  shape,  and  orientation.  For  this  reason,  u 
emphasized  by  Welch,  era/.  1980,  attempts  to  estimate 
the  optical  depth  of  clouds  composed  of  non-spherlcal  par¬ 
ticles  from  observed  cloud  reflectance  (or  transmittance) 
are  subject  to  considerable  error.  For  comparison  pur¬ 
poses,  estimates  were  derived  for  both  a  nominal  value  of 
.85  for  the  asymmetry  factor  in  cirrus  clouds  and  a  marked 
reduction  to  .75.  In  agreement  with  the  results  of  Welch, 
it  al.  1980,  the  calculated  optical  depths  corresponding  to 
the  same  average  cirrus  cloud  transmittance  differ  by  more 


COSINE  SOLAR  ZENITH  ANGLE 


Hi.  •.  Average  overcut/clear  Irradiance  ratios  from  SOLMET 
data  base  (solid  linns)  and  estimated  with  delta-Eddlngton  approx¬ 
imation  (dashed  lines).  Listed  are  the  values  of  optical  depth,  r, 
giving  the  beet  correspondence. 


Table  3. 

The  approximate  average  and  ±  standard  deviation 
values  of  optical  thickness  for  various  doud  forms  aa 
calculated  using  the  delta-Eddlngton  epproxlmation  and 
the  SOLMET  data  baat.  The  assumed  value  of  asym¬ 
metry  factor  la  given  in  column  g .  For  reference  pur- 
pom  a  value  of  scattering  ratio,  Q(5S0nm),  is  shown 
for  each  value  of  At  aa  determined  for  a  doud  of 
geometric  cloud  thickness  of  1  km  and  t  base  altitude  aa 
indicated. 


than  a  factor  of  2  for  the  two  assumed  values  of  g  as 
shown  In  Table  3. 

As  Indicated  by  the  examples  In  Fig.  6,  one  can 
obtain  a  rather  close  At  of  the  average  irradiance  ratios  vr 
solar  zenith  angle  for  the  standard  cloud  types  using  the 
delta- Eddington  model.  The  optical  depths  corresponding 
to  the  closest  approximations  of  the  SOLMET  data  aver¬ 
ages  for  each  cloud  type  are  listed  in  Table  2.  Rough  estl- 
matec  of  the  variability  In  the  optical  depth  were  obtained 
in  i  similar  manner  from  the  SOLMET  data.  The  tabu¬ 
lated  values  of  optical  depth  labeled  "±  rms  variation" 
were  calculated  from  SOLMET  Irradiance  ratios  equal  to 
the  average  ratio  for  a  given  cloud  type  and  solar  zenith 
angle  plus  and  minus  the  root-mean-square  departure  of 
the  Individual  SOLMET  values  of  Irradiance  ratio  'rom  the 
averages  for  that  cloud  type. 


APPENDIX  A 


Thus,  Tor  overcast  sky  conditions,  the  radiance  Reid 
at  any  altitude  may  be  calculated  using  the  optical  proper* 
ties  given  in  Table  2  corresponding  to  the  observed  cloud 
type  and  altitude  and  the  estimated  input  variables  for  the 
additional  atmospheric  layers  above  and  below  the  cloud 
layer.  The  optical  scattering  ratio  for  the  cloud  layer  Is 
given  by  the  cloud  optical  depth  divided  by  the  Rayleigh 
optical  depth  for  the  corresponding  altitude  interval  of  ‘he 
cloud.  In  view  of  the  uncertainty  in  the  asymmetry  factor 
for  cirrus  cloud  particle  distributions,  estimates  for  g  equal 
to  .75  are  shown  as  well  as  for  the  more  representative 
value  of  .85.  In  the  absence  of  other  Information,  opera¬ 
tional  use  of  the  results  for  g-,85  are  suggested, 

Several  comparisons  of  the  calculated  spectral  radi¬ 
ance  fields  with  the  high  resolution  aircraft  measurements 
for  both  clear  and  cloudy  atmospheres  are  given  in  the  fol¬ 
lowing  section, 

8.0  COMPARISONS  OF  MODEL  RESULTS 
WITH  MEASURED  DATA 

From  the  fundamental  Eq,  4,  we  note  that  the  relia¬ 
bility  of  model  calculations  of  slant-path  contrast  transmit¬ 
tance  are  dependent  directly  upon  the  accuracy  with  which 
the  radiance  fields  £,(*,«,<*»)  are  determined,  For  this 
reason,  continuing  emphasis  has  been  given  to  direct  com¬ 
parisons  of  the  radiance  fields  as  calculated  by  the  model¬ 
ing  procedures  and  the  corresponding  radiance  fields 
observed  as  part  of  the  extensive  airborne  measurement 
program.  The  upper  and  lower  scanning  radiometers 
measured  the  apparent  radiance  fields  surrounding  the  air¬ 
craft  as  it  probed  the  lower  troposphere  up  to  an  altitude 
of  6  km.  The  processed  radiometer  measurements  pro¬ 
vided  data  with  5-degree  resolution  In  the  4  specific 
wavebands  commensurate  with  the  integrating  nephelome- 
ter  measurements  of  total  volume  scattering  coefficient 
(peak  wavelengths  of  475,  550,  670  and  750  nm). 

Several  of  the  comparisons  between  model  calcula¬ 
tions  and  measured  radiance  fields  are  shown  in  Figs.  8  to 
12.  The  plotted  data  are  presented  Tor  either  one  of  two 
vertical  planes  (0-180  deg  or  90-270  deg)  with  respect  to 
the  solar  azimuth  and  depict  data  sweeping  from  the  nadir 
(«  -  180  deg)  up  through  the  horizon  and  zenith  and 
then  do  vnward  through  the  opposite  horizon  to  nadir,  as 
shown  In  Fig.  7.  A  series  of  examples  were  selected  in 
order  to  illustrate  model  performance  over  a  wide  range  of 
haze-layer  extinction,  solar  zenith  angle,  and  surface 
reflectance  (land  and  sea)  for  both  clear  and  overcast  sky 
conditions. 

The  specific  input  data  for  each  of  the  illustrated 
examples  are  listed  in  Appendix  B,  Except  for  the  scatter¬ 
ing  mixing  ratio  profile  in  the  lower  troposphere,  the  other 
input  variables  were  approximated  by  their  average  or 
climatological  values.  The  scattering  mixing  ratio  for  the 
boundary  layer  was  entered  as  determined  directly  from 
the  nephelometer  measurements  for  that  particular  flight 
profile.  The  scattering  ratio  for  the  upper  troposphere  and 
stratosphere  during  this  period  of  minimal  volcanic  activity 
was  held  constant  for  all  cases  and  set  equal  to  a  charac¬ 
teristic  value  of  1.2  at  550  nm,  The  mixing  ratio  for  other 


Fig,  7.  Coordinate  system  for  the  calculation  of  sky  and  terrain 
radiance  and  slant-path  contrast  transmittance. 


wavelengths  at  upper  levels  was  determined  assuming  the 
wavelength  dependence  of  the  aerosol  scattering 
coefficient  equal  to  (A/Xjjo)"1.  The  aeroaol  single  scatter¬ 
ing  phase  function  for  all  haze  layers  and  all  eases  was 
specified  by  a  two-term  Henyey-Oreensteln  function  with 
parameters;  g \  -  .714,  gj  -  -.613,  and  c  -  .963.  The 
corresponding  values  for  the  cloud  layers  are  ,890,  -.660 
and  .980, 

Optical  depths  for  the  overcut  sky  layers  were 
specified  equal  to  the  average  value  for  the  observed  cloud 
type  as  developed  in  Section  4  (see  Table  2).  The  single 
scattering  albedo  and  surface  reflectance  were  auumed 
equal  to  average  values  for  the  prevailing  conditions.  For 
flights  conducted  over  open  water,  the  average  Fresnel 
reflectance  distribution  for  wlnd-ruffied  sea  (see  Sec¬ 
tion  3.1)  was  used. 


5.1  Comparative  Results  for 
Individual  Flights 

Flight  C-379,  Fig.  8,  was  a  midday  flight  off  the 
south  coast  of  Lolland  Island,  Denmark.  The  example 
depicts  the  characteristic  changes  with  altitude  or  the 
upper  and  lower  hemisphere  radiance  fields  under  cloud¬ 
less  skies  and  light  haze  conditions.  Note  the  characteris¬ 
tic  uniformity  and  symmetry  of  the  sea  surface  radiance 
fields  in  the  cross-sun  plane. 

Flight  C-466,  Fig.  9,  was  conducted  over  open  farm¬ 
land  and  wooded  areas  In  northwestern  Germany.  Scat¬ 
tered  cirrus  clouds  and  Isolated  patches  of  altocumulus 
were  observed.  A  layer  of  moderate  haze  extended  to  an 
altitude  of  1350m.  Notice  the  marked  fluctuations  in  the 
observed  apparent  terrain  radiance  associated  with  varia¬ 
tions  in  ground  cover  (green  und  brown  fields  and  dark 
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woods),  as  opposed  to  the  smooth  terrain  radiance  distri¬ 
bution  calculated  with  the  assumption  of  uniform  Lamber¬ 
tian  surface  reflectance. 

Flight  C-469,  Figs.  10  and  11,  was  an  afternoon 
flight  over  open  farmland  in  central  Netherlands.  The 
examples  illustrate  the  variations  in  radiance  Aelds  with 
altitude  in  dense  haze  conditions.  Depicted  also  are  the 
characteristic  differences  in  the  radiance  distributions  as 
observed  in  the  narrow  blue  passband  (475  nm)  as  shown 
in  Fig.  10  and  in  the  near  infrared  passband  (750  nm) 
shown  in  Fig.  11.  Note  in  particular  the  reduction  in  over¬ 
head  sky  radiance  at  750  nm  associated  with  diminished 
molecular  and  aerosol  particle  scattering  and  the  enhanced 
apparent  terrain  racliance  due  to  the  greater  surface 
reflectance  in  the  near  infrared  portion  of  the  spectrum. 

Flight  C-401,  Fig.  12,  was  a  winter  day  flight  with 
low  sun  elevation  over  a  rural  area  in  northwestern 
France.  Scattered  cumulus  and  scattered  thin  cirrus  clouds 
were  observed  with  a  layer  of  light  haze  extending  to  an 
altitude  of  about  1  km.  The  relatively  bright  horizon  in 
contrast  with  low  overhead  sky  and  terrain  radiance  is  a 
prominent  feature  of  low-sun  radiance  fields.  As  Indicated 
by  this  example,  the  good  agreement  between  model  cal¬ 
culations  and  measured  radiance  fields  characteristically 
found  in  the  comparisons  extends  to  the  cases  with  low 
solar  elevation  angle. 

Flight  C-422,  Fig.  13,  is  also  an  over  water  flight  off 
the  coast  of  Denmark,  but  on  this  day  the  sky  was  over¬ 
cast  with  stratocumulus  clouds  having  a  base  altitude  of 
1800m.  Heavy  haze  was  present  at  all  levels  below  the 
cloud  base.  The  example  illustrates  the  highly  symmetric 
and  uniform  distribution  of  cloud  and  sea  radiance  typical 
of  measurements  made  with  the  flight  path  over  open 
water  and  under  an  overcast  cloud  layer. 

Flight  C-465,  Fig.  14,  was  an  afternoon  flight  In  mid 
summer  over  northwestern  Germany.  The  sky  was  over¬ 
cast  with  stratocumulus  clouds  having  a  base  altitude  of 
1600m.  A  layer  of  light  haze  was  observed  below  the 
cloud  base.  In  contrast  to  the  previous  case,  the  radiance 
field  was  more  irregular  in  structure  due  to  the  variations 
in  cloud  optical  thickness  and  to  the  variations  in  ground 
cover  reflectivity. 

Flight  C-43S,  Fig.  15,  was  a  winter  flight  over  a 
snow  covered  plateau  In  southwestern  Germany,  The  rk y 
was  overcast  with  altostratus  clouds  with  a  base  altitude  of 
3  km.  The  apparent  radiance  of  the  snow  surface  as 
measured  at  an  altitude  of  6  km  is  roughly  equal  to  the 
apparent  radiance  of  the  overcast  sky.  The  minimum 
values  of  lower  hemisphere  radiance  are  associated  with 
underlying  areas  of  dark  woods  having  an  effective  aver¬ 
age  reflectance  of  about  .04  in  this  case.  Notice  that  the 
effective  reflectance  of  the  snow  covered  areas  is  some¬ 
what  higher  than  the  average  value  of  0.75  that  was 
assumed  for  the  model  calculations. 

The  correspondence  between  the  Held  measure¬ 
ments  and  the  model  calculations  illustrated  above  are 
representative  in  general  of  the  many  such  comparisons 
made  during  the  course  of  model  development  and  valida¬ 


tion.  It  should  be  emphasized  that  the  accuracy  of  the 
illustrated  model  calculations  is  enhanced  considerably  by 
the  availability  of  the  detailed  vertical  profile  of  scattering 
coefficient  in  the  boundary  layer,  which  is  the  mqjor  deter¬ 
minant  of  the  fluctuations  in  path  radiance  in  the  tropo¬ 
sphere.  On  the  other  hand,  the  correspondence  between 
the  actual  and  calculated  radiance  fields  could  be  Improved 
through  relaxation  of  the  constraint  that  other  input  vari¬ 
ables  be  held  equal  to  their  typical  or  average  values.  For 
example,  Identification  of  an  aerosol  model  appropriate  for 
the  boundary  layer  in  a  given  synoptic  situation  would  pro¬ 
vide  estimates  of  the  single  scattering  phase  function  and 
single  scattering  albedo  that  are  more  realistic  than  the 
prescribed  overall  average  valuea  used  In  illustrated  com¬ 
parisons.  A  brief  assessment  of  the  sensitivity  of  esti¬ 
mates  of  the  maximum  detection  range  of  distant  objects 
to  uncertainties  in  both  atmospheric  and  target  factors  Is 
given  in  the  following  Section. 


6.0  SENSITIVITY  OF  TARGET  DETECTION  RANGE 
TO  ENVIRONMENTAL  AND  TARGET  FACTORS 

Diagnostic  calculations  have  been  made  routinely 
during  the  course  of  the  investigations  to  determine  the 
sensitivity  of  the  model  calculations  to  variations  In  the 
input  data.  In  prior  investigations  (Herlng,  1981),  atten¬ 
tion  was  directed  primarily  toward  analysis  of  the  sensi¬ 
tivity  of  contrast  transmittance  to  changes  or  uncertainties 
In  the  relevant  optical  properties  of  the  atmosphere.  The 
interpretations  of  these  results  are  straightforward  in  that, 
except  for  scintillation  effects,  Eq.  4  describes  completely 
the  atmospheric  effects  upon  contrast  transmittance.  On 
the  other  hand,  contrast  transmittance  must  be  coupled 
with  the  effects  of  specific  target  factors  to  estimate  the 
maximum  distance  an  object  can  be  detected  and 
identified.  A  given  change  in  the  slant-path,  contrast 
transmittance  may  or  may  not  be  critical  for  target  detec¬ 
tion  depending  the  minimum  contrast  required  for  detec¬ 
tion,  Thus,  the  assessment  of  the  combined  Influences  of 
atmospheric  and  target  factors  is  Important, 

For  this  purpose,  analytic  techniques  were  intro¬ 
duced  into  the  computor  program  to  calculate  visual  detec¬ 
tion  range  as  an  additional  output  variable,  The  represen¬ 
tation  of  target  information  and  visual  search  factors  in 
this  preliminary  analysis  are  greatly  simplified.  Neverthe¬ 
less,  the  results  provide  Important  insight  into  the  sensi¬ 
tivity  of  visual  detection  distance  to  changes  in  specific 
atmospheric  variables  for  selected  combinations  of  target 
characteristics. 

6.1  Determination  of  Visual 
Detection  Range 

The  analytic  representations  of  the  contrast  thres¬ 
holds  for  visual  detection  were  derived  from  basic  visual 
data  from  the  Tiffany  experiments  (Blackwell,  1946)  and 
the  Visibility  Laboratory  experiments  (Taylor,  1964)  which 
were  adapted  and  summarized  by  Gordon  (1979),  The 
graphical  representations  of  threshold  contrast  os  a  func¬ 
tion  of  angular  target  size  shown  in  Fig,  16  are  taken 
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Fig.  10.  Comparison  of  sky  and  terrain  radiance  as  measured  by  airborne  scanning  radiometer  (Johnson  198Id)  (solid  line)  and  model  cal¬ 
culations  (dashed  tine). 


Fig.  11.  Comparison  or  sky  and  terrain  radiance  as  measured  hy  airborne  scanning  radiometer  (Johnson  I98td)  (solid  line)  and  model  cal¬ 
culations  (dashed  line). 


A-19 


V 


■v 

1  ■  *m  .  ►,  •  (*.'(< 


APPENDIX  A 


,  a| 

ql 

lGl 


\ 

1 


C-401  MUZ,  PRANCE 

CLEAR  SKY 

LOW  SUN  BLEVATION . 

S50nm 

#,  -  76  d«t. 


|  10J: 

Ti  : 

c 

* 


*-1*0* 


•90  6  90 

ZENITH  OBSERVATION  ANGLE  (Degrees) 


Fig,  12.  Comparison  of  sky  and  terrain  radiance  as  measured  by  airborne  scanning  radiometer  (Johnson  1981d)  (solid  line)  and  model  cal¬ 
culations  (dashed  line). 
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Fig.  13,  Comparison  of  sky  and  terrain  radiance  as  measured  by  airborne  scanning  radiometer  (Johnson  1981d)  (solid  line)  and  model  cal 

culations  (dashed  line). 
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n«.  1C.  Angular  subtense  of  Uriel  u  a  function  or  threshold 
contrut  for  a  99  pet  cent  probability  (confidence)  of  detection  and 
a  lack  of  knowledge  of  target  pod t Ion  of  ±  4,  deg  or  more  (Gor¬ 
don,  1979).  iBf  li  the  adaptation  brightneai  and  At  Is  the 
tiaumed  visual  fixation  (dwell)  lima. 


directly  from  Fig.  1  of  Gordon  (1979).  Two  search  modes 
are  represented.  The  first  corresponds  to  Ions  duration 
viewing  (>24  seconds,  TlfTany  data).  The  second 
corresponds  to  a  viewing  time  of  one-third  of  a  second 
(Visibility  Laboratory  data),  which  is  the  approximate 
fixation  time  associated  with  search  performed  by  the 
human  eye  (White,  1964). 

From  Fig.  16  we  see  that,  for  .the  normal  range  of 
daylight  conditions,  the  visual  contrut  threshold  does  not 
vary  significantly  With  background  luminance.  However, 
the  contrut  threshold  does  vary  markedly  with  the  angular 
subtense  of  the  target.  The  relationships  shown  in 
Fig.  16,  represent  the  contrast  threshold  associated  with 
confident  detection  (99  percent  probability)  under  condi¬ 
tions  where  the  lack  of  knowledge  of  target  position  is  ±  4 
degrees  or  more. 

For  application  to  the  problem  at  hand,  the 
prescribed  relationships  given  in  Fig.  16  were  approxi¬ 
mated  analytically  by  the  expression, 


In  y  -  /f+i»(-ln  Cr)  +  C(-ln  Cr)2  +  0(-ln  Cr)3  (30) 


where  y  is  the  angular  size  of  the  target,  Cr  is  the  con¬ 
trast  detection  threshold  and  A,B,C  and  D  are  the  paran  j- 
terization  coefficients  as  determined  for  the  two  search 
modes. 


Given  the  above  method  for  the  representation  of 
target  factors,  let  us  now  consider  analytic  representation 


of  the  model  calculations  of  slant  path  contrast  transmit¬ 
tance.  The  computational  scheme  provides  discrete  values 
of  contrut  transmittance  for  designated  paths  of  sight.  As 
an  additional  option,  the  contrast  transmittance  of  air  to 
ground  observation  paths  for  fixed  values  of  azimuthal 
viewing  angle  and  sensor  altitude  is  represented  through 
the  method  of  least  squares  as  a  continuous  function  of 
path  length  (range)  by  the  expression, 


t  -  o+Mln  7 V  (*£))  +  clln  Cr,®))2 


(31) 


Since  the  angular  subtense  of  the  target  is  1ms  than  a 
degree  or  so,  the  slant  range,  r,  can  be  approximated  by 
r  -  d/y  where  d  is  the  effective  diameter  of  the  target. 
Given  the  inherent  target  contrut,  C» ,  the  range  auoci- 
ated  with  a  selected  value  of  apparent  contrut  C, ,  is 
determined  by  Eq.  31.  The  maximum  detection  range  for 
specified  sensor  altitude  and  azimuthal  viewing  angle 
corresponds  to  the  range  where  C,  -  C,.  The  method  of 
successive  approximation  yields  rapidly  converging  solu¬ 
tions  for  maximum  detection  range  from  trial  values  of  C, 
and  Q  using  Eqs.  30  and  31. 


6,2  Diagnostic  Calculations  of 
Visual  Detection  Distance 

A  seriu  of  additional  model  calculations  were  car¬ 
ried  out  to  illustrate  the  relative  sensitivity  of  target  detec¬ 
tion  range  to  changes  or  uncertainties  in  some  of  the  key 
variables.  For  this  purpose,  reference  conditions  were 
auumed  which  reflect  a  typical  visual  detection  scenario 
for  hazy  atmospheres  in  the  daytime.  In  an  attempt  to 
Isolate  the  effects  of  fluctuations  in  individual  factors,  the 
procedure  wu  to  vary  a  given  parameter  over  its  normal 
range  while  holding  all  other  reference  conditions  con¬ 
stant.  The  reference  atmospheru  used  in  this  analysis  are 
summarized  in  Appendix  B  (Ref.  A  for  clear  sky  and 
Ref.  B  for  overcut  sky),  The  reference  atmospheres  have 
a  surface  equivalent  visibility  of  20  km  and  a  boundary- 
layer  hazy  depth  of  1  km.  Except  u  noted,  an  ideal  circu¬ 
lar  target  is  assumed  having  a  diameter  of  30m  and  an 
inherent  contrut  of  one.  The  functional  representation  of 
psychophysical  data  of  Taylor  (1964)  corresponding  to  the 
0.33  sec  fixation  time  (see  Fig.  16)  was  used  for  the  trial 
calculations  except  for  the  one  comparative  example  (para¬ 
graph  6.2.5)  based  upon  the  Tiffany  visual  threshold  data 
for  a  fixation  time  greater  than  24  seconds. 


Visual  Horizontal  Detection  Distance 
versus  Observer  Altitude  and  Azimuthal 
Viewing  Angle  (Fla.  17), 


For  Intermediate  size  targets  of  good  inherent  con¬ 
trast,  the  maximum  horizontal  detection  distance  (ground 
range)  Increases  significantly  with  increasing  observer  alti¬ 
tude  in  the  lower  troposphere  above  the  boundary  layer. 
The  variations  in  ground  range  with  azimuthal  viewing 
angle  are  determined  largely  by  the  directional  dependence 
of  the  path  radiance  as  defined  by  the  combined  phase 
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Fl«.  17.  Vliual  horizontal  detection  dliUnoe  u  a  function  of  Fig.  11.  Vtaual  horizontal  detection  distance  at  a  function  of 
observer  altitude  and  azimuthal  observation  angle.  Inherent  contrait  for  30-m  target. 


function  for  molecular  and  aerosol  scattering.  For  the  typ¬ 
ical  condition  in  the  free  atmosphere  above  the  haze  layer 
where  the  aerosol  particles  have  low  concentration,  the 
difference  in  detection  distance  in  the  upsun  (4  -  0  deg) 
and  downsun  (<t>  m  180  deg)  viewing  directions  remains 
rather  constant  with  increasing  observer  altitude  above  the 
boundary  layer  and  increases  with  increasing  solar  zenith 
angle.  Of  special  interest  In  Fig.  17  is  the  relatively  small 
variation  in  detection  range  with  azimuthal  angle  in  the 
antisolar  hemisphere  as  indicated  by  the  rather  close 
correspondence  in  the  180  degree  (downsun)  and  90 
degree  (cross  sun)  profiles.  It  should  be  noted  that  the 
trial  calculations  which  follow  in  subsequent  examples 
assume  a  90-degree  azimuthal  viewing  direction. 


Visual  Detection  Distance  versus 
Inherent  Contrast 


J3T1EJ 


The  rate  of  increase  in  visual  detection  range  with 
increasing  observer  altitude  depends  significantly  upon  the 
initial  target  contrast.  Notice  in  Fig.  18  the  dramatic 
Increase  in  detection  range  with  Increasing  viewing  altitude 
for  bright  targets  viewed  against  a  dark  background 
(C0  -  5)  as  contrasted  with  the  restricted  increase  In 
detection  distance  for  low  contrast  targets  (C«  -  0,2).  As 
a  point  of  reference  information,  the  detection  range  is 
independent  of  the  sign  of  the  target  contrast.  For  targets 
darker  than  the  background,  the  lower  contrast  limit  Is 
C,  -  - 1  (black  target).  The  inherent  contrast  of  targets 
brighter  than  background  is  unbounded.  The  assumed 
reference  value  of  target  contrast  of  C0  -  1  (dashed 
curve)  corresponds  to  the  case  where,  under  the  same 
irradiance,  the  reflectivity  of  the  target  is  twice  the  back¬ 
ground  reflectivity. 


Fig.  19.  Visual  horizontal  detection  distance  as  a  (Unction  of 
inherent  contrast  for  S-m  target. 


in  visual  detection  distance  as  compared  with  the  results 
for  the  30-m  target  (Fig.  18).  Notice,  for  example,  In  the 
case  of  a  target  of  0.2  inherent  contrast  under  the  assumed 
conditions,  there  is  little  advantage  in  increased  detection 
range  to  be  gained  by  climbing  to  an  observation  altitude 
above  the  haze  layer  top  and  that  the  reference  target  hav¬ 
ing  0,2  contrast  cannot  be  detected  with  99  percent  proba¬ 
bility  above  a  viewing  altitude  of  6  km.  The  calculated 
visual  detection  range  for  5-m  targets  for  this  example  is 
everywhere  less  than  12  km  for  targets  having  an  inherent 
contrast  less  than  one. 


Visual  Detection  Distance  versus 
Inherent  Contrast  for  Reduced 
Target  Size  (Fig.  19), 

Diagnostic  model  calculations  for  a  smaller  S-m  tar¬ 
get  diameter  (Fig.  19)  show  the  corresponding  reduction 


Visual  Detection  Distance  as 
a  Function  of  Environmental 
Factors  (Fia,  20t 

Sensitivity  crtculations  of  the  response  of  ground 
detection  range  to  selected  changes  in  Individual  factors 
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Fig.  20.  Sensitivity  of  horizontal  detection  distance  to  changes  In  relevant  factors.  Solid  bar  In  each  case  (22  km)  corresponds  to  the  calcu¬ 
lation  for  reference  target  (30-m  diameter)  and  reference  atmosphere  parameters.  Surface  target  with  a  contrast  of  one  Is  viewed  against  at 
a  uniform  background  from  an  altitude  of  3  km,  For  each  case,  all  variables  are  fixed  equal  to  reference  values  except  for  the  variable 
under  examination. 


are  shown  in  Fig.  20  for  a  ground  target  viewed  from  an 
altitude  of  3  km.  All  variables  except  the  factor  under 
examination  are  held  constant  and  equal  to  the  reference 
values  listed  in  paragraph  6.1  and  Appendix  B  (Ref,  A). 
The  scattering  and  absorption  coefficients  are  allowed  to 
change  only  within  the  primary  haze  layer.  The  calcula¬ 
tions  assume  a  A.-1  wavelength  dependence  of  the  aerosol 
scattering  coefficient  and  no  changes  In  surface  reflectance 
with  wavelength.  The  illustrated  changes  in  detection  dis¬ 
tance  relative  to  changes  in  single  scattering  albedo  were 
calculated  assuming  that  the  optical  thickness  of  the  haze 
layer  remains  constant,  i.e,  the  scattering  ratio  was 
changed  by  an  amount  sufficient  to  compensate  for  the 
postulated  change  In  absorption  coefficient. 

The  solid  bar  In  each  of  the  graphs  shown  In  Fig,  20 
depicts  the  ground  range  (22  km)  associated  with  the 
reference  values  for  the  Individual  variables.  Notice  for 
the  prescribed  conditions  that  a  decrease  in  target  distance 
of  about  25  percent  is  associated  with:  (Da  reduction  In 
target  size  by  50%;  (2)  a  reduction  In  surface  equivalent 
visibility  by  35%;  (3)  an  increase  in  solar  zenith  angle 
from  60  to  78  degrees;  (4)  a  decrease  in  surface 
reflectance  by  50  percent,  and  (5)  a  change  In  wavelength 
from  550  to  450  nm,  The  detection  distance  is  less  sensi¬ 
tive  to  changes  in  the  absorption  contribution  to  a  fixed 
value  of  total  volume  extinction  coefficient  In  the  boun¬ 
dary  layer.  It  should  be  emphasized  that  the  results 
shown  In  Fig,  20  are  specific  only  for  the  conditions 
assumed  for  the  diagnostic  calculations.  Simultaneous 
changes  in  two  or  more  variables  could  of  course  result  in 


a  much  broader  range  in  the  variation  of  target  detection 
distance. 

Visual  Detection  Distance  versus 
Visual  FIxation  lltTijlFIg,  2.1 )~ 

The  variations  in  target  detection  distance  that  are 
associated  with  changes  in  stimulus  duration  are  shown  In 
Fig.  21,  The  two  functional  relationships  for  contrast 
threshold  data  used  for  the  calculations  are  the  same  as 
those  depicted  in  Fig.  16.  As  discussed  above,  the  first 
applies  to  unlimited  viewing  time  (>24  sec)  and  the 
second  applies  to  the  average  fixation  time  (1/3  sec)  for 
the  performance  of  visual  search  by  the  human  eye. 
Comparative  calculations  for  the  two  viewing  conditions 
are  shown  for  two  atmospheres.  Both  atmospheres  and 
target  factors  are  the  same  as  reference  conditions  cited 
above  except  that  the  haze  layer  of  depth  1.3  km  is 
represented  by  the  LOWTRAN  (Shettle  and  Fenn,  1979) 
rural  model  atmosphere  for  relative  humidity  equal  to  70 
percent  for  the  one  case  and  for  relative  humidity  equal  to 
99%  In  the  other  case.  Notice  that  the  relative  difference 
in  detection  distance  versus  glimpse  time  for  the  example 
with  high  humidity  and  high  extinction  coefficient  Is 
roughly  the  same  as  that  for  the  lower  humidity  and  lower 
extinction  coefficient  example, 

Visual  Detection  Distance 
for  Overcast  Sky  (Fig,  22), 

A  cloudless  atmosphere  was  assumed  for  the  sensi¬ 
tivity  analyses  described  above,  Let  us  now  consider 
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Fl|.  21.  Vliuil  horizontal  detection  diitance  versus  altitude  for 
fixation  time*  (A/)of  0.33  aec  and  Infinity.  Calculation!  assume 
LOWTRAN  rural  model  for  boundary  layer  with  70  percent  humi¬ 
dity  (equivalent  aurfaca  vlilbllity  of  about  20  km)  and  with  99  per¬ 
cent  humidity  (equivalent  surface  visibility  of  about  S  km). 


briefly  the  behavior  of  contrast  transmittance  and  the  vari¬ 
ations  in  target  detection  diitance  that  are  associated  with 
overcast. sky  conditions.  In  the  region  below  the  cloud 
layer  in  the  daytime,  both  the  inherent  background  radi¬ 
ance  *4 (z, ,0, 0)  .and  the  apparent  radiance  tL,  (z,0,  <f>) 
are  essentially  proportional  to  the  downwelling  radiance  at 
the  cloud  base  when  the  direct  sun  component  of  the  irra- 
dlance  is  small  in  comparison  with  the  diffuse  component. 
Thus,  except  for  thin  clouds  with  significant  direct  sun 
penetration,  the  contrast  transmittance  of  slant  paths 
belovy  the  cloud  deck  Is  to  a  good  approximation  indepen¬ 
dent  of  the  downwelling  irradiance  from  the  cloud  layer 
(see  Eq.  4).  However,  the  path  radiance  and  radiance 
transmittance  (hence  the  contrast  transmittance)  are  of 
course  strongly  dependent  upon  the  optical  properties 
along  the  cloud  free  paths  of  sight  below  the  overcast. 

A  comparison  of  target  detection  distance  as  calcu¬ 
lated  with  and  without  an  overcast  cloud  are  shown  in  Fig. 
22.  The  diagnostic  calculations  again  were  made  using 
both  the  70-perr‘ent  and  99-percent  humidity  LOWTRAN 
rural  models  for  the  boundary  layer.  For  the  cloud  free 
cases,  calculations  are  shown  for  solar  zenith  angles  of  25, 
65  and  85  degrees.  A  solar  zenith  angle  of  55  degrees  and 
stratocumulus  overcast  of  average  optical  depth  were 
assumed  for  the  calculation  of  the  cloud  case.  As  dis¬ 
cussed  above,  the  slant  path  contrast  transmittance  and 
detection  distance  for  the  overcast  condition  are  not  sensi¬ 
tive  to  cloud  optical  depth  (or  solar  zenith  angle) .  Notice 
in  Fig.  22  that  the  calculated  detection  distance  for  the 
cloud  case  corresponds  to  the  case  where  the  cloud  is 
removed  and  the  solar  zenith  angle  Is  near  75  degrees. 
This  Is  a  typical  result  for  conditions  where  a  Lambertian 
surface  reflectance  is  assumed  and  the  optical  properties 
along  the  viewing  path  are  the  same  for  both  the  clear  and 
cloud  examples. 

7.0  SUMMARY  COMMENTS 

The  analytic  techniques  that  were  developed  as  part 
of  the  experimental  optical  measurement  and  analysis  pro¬ 


ng.  22.  Comparison  of  calculated  targtt  detection  distance  for  a 
reference  atmosphere  with  and  without  m  overcast  cloud  layer. 
Calculations  for  dear  sky  case  (solid  lines)  are  shown  for  aoler 
zenlth  angles  of  25,  (5,  and  IS  deg.  Calculations  for  esse  with 
stratocumulus  overcast  (dashed  line)  are  for  soler  zenith  angle  of 
35  deg. 

gram  provide  relatively  fast  end  consistent  estimates  df  the 
visible  spectrum  contrast  trsnemituutce  along  vertically 
Inclined  paths  of  eight.  Essential  elements  of  the  scheme 
include  modeling  techniques  suitable  for  field  use  for  the 
determination  of  the  apparent  background  sky,  cloud  and 
terrain  radiance  distribution  as  well  as  the  radishes 
transmittance  along  the  viewing  path.  Emphasis  was 
placed  upon  the  development  of  algor  ithms  which  are  sim¬ 
ple  end  easy  to  apply  yet  are  adaptable  so  aa  to  take  advan¬ 
tage  of  all  relevant  observations  and  forecast  information. 
The  procedural  developed  for  the  calculation  of  clear  sky 
radiance  distribution  are  directly  applicable  to  the  overcast 
sky  condition  through  introduction  of  the  appropriate  opti¬ 
cal  propartiaa  of  the  cloud  layer.  Comparisons  of  model 
calculations  with  experimental  aircraft  measurements  of 
cloud  radiance  indicate  that  use  of  average  optical  depths 
for  various  cloud  typaa  as  derived  from  the  SOLMET  data 
base  yields  a  good  drat  approximation  to  the  actual  over¬ 
cut  cloud  radiance  distribution. 

Another  recant  extension  to  the  image  transmiuion 
modeling  procedures  is  the  addition  of  techniques  to  deal 
with  the  combined  influences  of  itmospheric  and  target 
factors  in  order  to  determine  target  detection  range  for 
objects  located  at  the  surface  and  viewed  from  eloft. 
Although  the  reiulte  described  in  this  report  pertain  to 
simple  targets  with  uniform  backgrounds,  the  diagnostic 
model  calculations  provide  insight  into  the  relative  Impor¬ 
tance  of  the  uncertainties  Involved  in  the  estimates  of  tar¬ 
get  detection  distance.  Not  surprisingly,  the  determina¬ 
tions  are  most  sensitive  to  uncertainties  in  the  optical 
depth  of  the  atmosphere  which  in  turn  impact  directly  all 
three  components  of  contrut  transmittance,  it  the 
inherent  and  apparent  background  radiance  and  the  path 
transmittance.  The  sensitivity  calculations  also  confirm 
the  strong  dependence  of  target  detection  distance  on 
knowledge  of  target  and  background  reflectivity.  The 
model  capability  for  joint  assessment  of  target  and 
environmental  effects  on  target  detection  range  also 
reveals  the  interactive  effects  of  such  parameters  u  solsr 
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zenith  angle,  surface  reflectivity,  target  size,  and  visual 
fixation  time  for  overcast  and  clear  sky  conditions. 

Efforts  are  underway  to  improve  the  methodology 
for  estimating  image  transmission  characteristics  by  relaxa¬ 
tion  of  some  of  die  constraints  which  limit  the  applicability 
snd  accuracy  of  the  techniques.  It  is  important  in  particu¬ 
lar,  to  develop  techniques  to  exploit  expanding  knowledge 
of  the  directional  dependence  of  the  reflectivity  of  various 
types  of  ground  cover.  Work  is  also  underway  to  incor¬ 
porate  analytic  techniques  to  estimate  slant-path  contrast 
transmittance  and  target  detection  distance  for  inclined 
paths  of  sight  in  fractional  cloud  cover  conditions. 
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SUMMARY  OF  INPUT  DATA, 
COMPUTATION  STEPS  AND 
RESULTANT  DATA 


A.1  General  Inpat  Data 

A.  1 . 1  Number  of  atmoapheric  layer* 

A.  1.2  Representative  wavelength  (pm) 

A.  1.3  Extraterrestrial  solar  (lunar)  Irradiance 
0 V/mfym) 

A.  1 .4  Rayleigh  optical  depth  of  top  layer 

A.  1 .3  Solar  zenith  angle*  (deg.) 

A.  1 .6  Zenith  viewing  angles  (deg.) 

A.  1.7  Average  reflectance  of  underlying  surface  and 
reflectance  of  background  In  the  vicinity  of  target 
if  different  than  the  average  reflectance. 

A.2  Input  Data  Each  Layer 

A.2.1  Base  altitude  (km) 

A, 2.2  Scattering  mixing  ratio 

A.2.3  Single  scattering  albedo 

A.2.4  Single  scattering  phase  function  (normalized  to 

integral  of  one) 

A.3  Computation  Steps 

A.3.1  Compute  Rayleigh  atmosphere  optical  thickness 
for  each  layer 

A.3.2  Compute  total  optical  thickness  for  each  layer. 

A.3.3  Compute  optica!  depth  for  the  base  of  each  layer. 

A.3.4"  Compute  'single  scattering  phase  function  and 
asymmetry  factors  for  each  layer.  (These  param¬ 
eters  may  be  specified  if  desired  in  lieu  of  the,  use 
of  model  algorithms). 

A.3.S  Transform  optical  depth,  single  scattering  albedo, 
and  asymmetry  factor  in  accordance  with  delta- 
Eddingion  approximation, 

A.3.6  Compute  path  function  distribution  for  the  base 
of  each  sub-layer  (each  layer  is  divided  equally 
into  2  sub-layers). 

A. 3.7  Compute  path  radiance  distribution  for  each  sub¬ 
layer. 

A.3.8  Compute  path  transmittance  distribution  for  each 
sub-layer. 

A.3. 9  Compute  sky/terrain  radiance  distribution  at 

sub-layer  boundaries. 


A.3, 10  Compute  slant  path  contrast  transmittance  of 
selected  paths. 


A.4  Output  Tables  aai  Data 

A.4.1  Path  radiance,  for  O*<0<B5*  (upward 

looking  paths)  for  selected  levels  (max  3)  for 
each  0,  (i.e.  Sky  Radiance)  (Wf mfam  ar) 

A.4, 2  Path  radiance,  for  95*<0<  180* 

(downward  looking  paths)  for  selected  levels 
(max  S)  for  each  0,  (W/m^m  ar) 

A, 4.3  Path  plus  terrain  radiance  L,  (xftji)  for 
95*<0<  ISO*  for  selected  levels  (max  5)  for  each 
0,  ar) 

A.4.4  Control  transmittance,  Tc(*fi&)  between  object 
and  sensor  altitudo  (max.  5  altitude  intervals) 
for  O*<0<  85*  and  for  each  0,. 

A.4, 5  Control  transmittance  Tc(x#&)  between  ohiect 
and  sensor  altitudes  (max.  5)  for  95*<0<  180* 
and  for  each  0, . 


A.9  Additional  Output  From  The  Supplemental 
Eddington  Computer  Program 

The  standard  Eddington  computer  program  is 
employed  o  an  integrated  part  of  the  model  calculations 
of  contrast  transmittance.  The  following  output  from  the 
program  as  modified  by  the  delta-Eddington  transforma¬ 
tion  of  input  parameters  is  available  in  tabular  form  for 
the  base  of  each  sub-layer  and  for  each  solar  zenith  angle: 

'  optical  depths,  t  and  r' 

•  components  of  Eddington  diffuse  radiance,  Lp Or)  and 

Lp'U)  ar) 

•  total  scalar  irradiance  (W/m^m) 

•  downward  diffuse  irradiance  (W/m^m) 

•  upward  diffuse  irradiance  (W/m^m) 

•  solar  irradiance  (W/m^m) 

'•  total  downward  irradiance  (W/m^m) 
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SPECIFIC  INPUT  DATA  FOR  MODEL  CALCULATIONS 
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GLOSSARY  AND  NOTATION 


The  notation  uaed  in  reporta  and  Journal  article!  produced  by  the  Vlaibility  Laboratory  atalT  followa,  in  general,  the 
rulee  set  forth  In  pages  499  and  S00,  Duntley  tl  al.  (19S7).  These  rules  are: 

•  Each  optical  property  is  indicated  by  a  baaic  (parent)  symbol. 

•  A  presubecript  may  be  used  with  the  parent  symbol  as  an  identifier,  «.#„  b  indicates  background  while  t  denotes 
an  otysct,/.*.  target. 

•  A  postsubscript  may  be  used  to  Indicate  the  length  of  a  path  of  sight,  r  denotes  an  apparent  property  as 
measured  at  the  end  of  a  path  of  sight  of  length  r,  while  o  denotes  an  inherent  property  based  on  the  hypothetical 
concept  of  a  photometer  located  at  zero  distance  from  an  object,  /.#,  target. 

•  A  poetsuperscript*  or  postsubscript*,  Is  employed  as  a  mnemonic  symbol  signifying  that  the  radiometric  quantity 
has  been  generated  within  the  path  or  path  segment  by  the  scattering  of  ambient  light  reaching  the  path  from  all 
directions  and/or  by  emission. 

•  The  parenthetical  attachments  to  the  parent  symbol  denote  altitude  and  direction.  The  letter  z  Indicates  altitude  in 
general;  t,  Is  uaed  to  specify  the  altitude  of  a  target.  The  direction  of  a  path  of  sight  Is  speoiAed  by  the  zenith 
angle  t  and  the  azimuth  4.  In  the  case  of  irradiances,  the  downwelling  Irradianoe  Is  designated  by  d,  the 
upwolling  by  u. 

•  The  radiometric  symbols  used  herein  now  correspond  to  the  OSA  recommendations  in  Section  1  of  Driscoll  and 
Vaughn  (197S).  Prior  to  June  1980,  the  symbol  used  for  radiance  L  was  Af,  for  Irradianoe  E  was  /f,  and  for 
attenuation  length  £  was  L, 


Symbol  Units  Quantity 


A(r) 

none 

Albedo  at  altitude  z,  defined 

A(z)  -  £(z,u)/  £(z,rf), 

ait) 

Absorption  coefficient. 

Cit ) 

sr~l  m~l 

Diffuse  component  of  the  solar  almucantar  sky  radiance  ratio  to  the  sun  scalar 
Irradiance  and  relative  airmass 

C(z)  -  J Lizt$',4>')<r(i,p')d(l  /  ,a(z)  +  L»  it)  /  ,*it), 

c 

m/s 

Speed  of  light  c  -  2.997924S8£8. 

C\ 

Wm 1 

First  radiation  constant  ct  «  lithe1  -  3.741 832 £-16, 

Ct 

mK 

Second  radiation  constant  C]  -  ch/k  -  1, 438786£-2. 

Dit) 

none 

Radiance  distribution  function  Dit)  -  j LitflW)stc8‘d fl  /  tit), 

d(k) 

none 

Limb  darkening  factor  relating  the  average  sun  radiance  and  the  center  sun 
radiance  Dik)  -  ,L0  /  ,L0, 

Ek 

W/m1 

Spectral  irradiance  (formerly  symbol  H)  defined  as  Ek  «  J Lk(z,M)  cos 9  dCl, 
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Symbol  Uniti 


Quantity 


E 

W/m'nm 

Broad  bind  sensor  Irradiance,  defined  as  £  -  J  £\ SKTk  rfX  /  tX. 

0 

E(x,d) 

W/m 1 

1  ' 

Irvadiancc  produced  by  downwelling  flux  as  determined  on  a  horizontal  flat  plate 
at  altitude  z  (formerly  /f(z,rf».  In  thia  report  d  ia  used  in  plaoo  of  the  minus 
sign  in  the  notation  lH(x,-)]  which  appura  in  Duntley  (1969).  This  property 
may  be  deAned  by  the  equation  £(z,d)  -  J L(z,0',<p')  cosd'dO', 

E{x,u) 

W/m1 

Irradiance  produced  by  the  upwelling  flux  as  determined  on  a  horizontal  Aat  plane 
at  altitude  i  (formerly  tf(r,u)).  Here  u  is  substituted  for  the  plus  sign  formerly 
used  in  the  notation  (//<*,+)]. 

//(*) 

m 

Scale  height  at  altitude  z,  the  height  of  a  homogeneous  atmosphere  having  the 
density  of  the  layer  at  altitude  z. 

h 

J/s 

Planck  conitant  h  «  6,626l76£-34. 

k 


J/K 


Boltzmann  conitant  k  -  1,3806621-23, 


Lk 


W/sr 


Spectral  radiance  (former  symbol  N). 


L 


W/sr  minm 


Broad  band  aenior  radiance  la  defined  as  L 


L0(x„9,4>)  W/sr  m1  Inherent  radiance  based  on  the  hypothetical  concept  of  a  photometer  located  at 

zero  distance  from  an  object  at  attitude  z  In  the  direction  ipedAad  by  zenith  angle 
9  and  azimuth  d>, 


Lf(i,9<d>)  Whr  m 1  Apparent  radiance  as  determined  at  altitude  z,  from  tho  end  of  a  path  of  sight  of 

length  r  at  zenith  angle  9  and  azimuth  d>.  This  property  may  be  defined  by 
L,(z,0,0)  -  L0(i,fid)Tt(t,9)  +  L'{x,9,d>). 


L(k,T) 


W/sr  mJ 


Black  body  radiance  at  wavelength  X  and  temperature  T, 


Lq{z,9,<t>)  W/sr  m1  Equilibrium  radiance  at  altitude  x  with  the  direction  of  the  path  of  sight  specified 

by  zenith  angle  9  and  azimuth  <t>.  Hits  property  Is  a  point  function  of  position 
and  direction, 


wJaiwAseiweaaJI'wi’i^Ji  i«i»«  l.'i?**!. ", 


As  discussed  by  Duntley  et  al.  (1957),  many  Image  transmission  phenomena  are 
most  clearly  understood  in  terms  of  the  concept  of  equilibrium  radiance.  This 
concept  is  a  natural  consequence  of  the  equation  of  transfer  which  states 
analytically  that  in  any  path  segment  the  difference  between  the  output  and  Input 
radiances  is  attributable  to  a  gain  term  and  a  toss  term,  such  that  some  unique 
equilibrium  radiance  L4(z,0,^)  must  exist  at  each  point  such  that  the  loss  of 
radiance  within  the  path  segment  is  exactly  balanced  by  the  gain,  l,e., 
AZ,(z,fl,<p)  -  0. 
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Symbol 


Quantity 


By  virtue  of  this  concept  and  the  equation  of  transfer,  one  can  ihow  that  each 
segment  of  every  path  of  sight  has  associated  with  it  an  equilibrium  radiance,  and 
that  the  space  rate  of  change  in  image  forming  radiance  caused  by  the  path 
segment  is  in  such  a  direction  as  to  cause  the  output  radiance  to  be  closer  to  the 
equilibrium  radiance  than  is  the  input  radiance.  This  segment  by  segment 
convergence  of  the  apparent  radiance  of  the  object  held  to  the  dynamic 
equilibrium  radiance  was  clearly  illustrated  by  the  data  In  the  1957  paper 
referenced  above. 


£•<*,#.«) 


W/sr  ml 


L'/z.tf,*)  W/sr  mi1 


Path  function  at  altitude  z  with  the  direction  of  the  path  of  sight  speuified  by 
zenith  angle  9  and  azimuth  d>  This  property  is  defined  by  the  equation 
!,(*,«, si)  -  lf|(z,0,^)  +  L,a(x). 

Path  function  due  to  shattering,  defined  by  £.*(:,», d)  -  ^  <r(z,/3'>L(*,0',d')tffl'. 


L»a  Or)  W/sr  m1 


Path  function  due  to  eminence  (absorption)  L*(z)  -  u(z)  i(X,7). 

m 


L',{ts9A) 


W/sr  m1 


Path  radiance  as  determined  at'  altitude  *  at  the  end  of  a  path  of  sight  of  length  r 
In  the  direction  of  zenith  angle  6  and  azimuth  d- 


Ll  (z,0,d) 


W/sr  m1 


Sky  radiance  at  altitude  z,  zenith  angle  9  and  azimuth  Also  the  path  radiance 
for  the  path  of  sight  of  length  a«  from  out  of  the  atmosphere  to  altitude  z. 


fL«(z,0,,d)  W/sr  m J 


Apparent  radiance  ot  the  center  of  the  solar  disk  as  determined  from  the  end  of  a 
path  of  sight  of  length  *>  from  out  of  the  atmosphere  to  altitude  z  at  the  zenith 
angle  of  the  sun  0,, 


Attenuation  length  at  altitude  z,  The  attenuation  length  is  the  distance  at  which 
the  signal  is  attenuated  to  l/e. 


M(k%T) 


black  body  exiiance  (emittance)  M(k,T)  »  nL(k,T), 


m(h,T) 


Black  body  scalar  exltance  (emittanco)  •»  4wZ,(\,D. 


Scalar  exltance  (emittance)  per  length  m*(z)  -  J  £,(z,0,d)</r>. 

Scalar  exltance  (emittance)  per  length  due  to  scattering 

m ,t{z)  -  J L.,(z,0,d)dn  -  «(z)j(z). 


m.a(z) 


Scalar  exltance  (emittance)  per  length  due  to  absorption  m.  (J  -  a(,z)m(k,T). 


m„(z,0) 


Absolute  air  mass  at  angle  9  m»(i,0)  -  f  p(z)<fr. 

I 
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Symbol 


Quantity 


none 


Relative  optical  air  mass. 


Relative  optical  airmass  at  zenith  angle  of  the  sun,  shorthand  for 

Refractive  modulus  V(r)  -  [11(2)  -  l]£6. 

Refractive  index. 

Optical  scattering  mixing  ratio  at  altitude  2.  This  quantity  is  defined  as  the  ratio 
of  the  total  volume  scattering  coefficient  at  attitude  z*  to  the  molecular  (or 
Rayleigh)  volume  scattering  coefficient  at  the  same  altitude  2. 
Q(z)  -  s(z)/*s(z). 


2<z,0) 


Volume  scattering  function  ratio  it  attitude  2.  This  quantity  is  defined  as  the  ratio 
of  the  total  volume  scattering  function  at  altitude  2  and  scattering  angle  0,  to  the 
molecular  (or  Rayleigh)  volume  scattering  function  at  the  same  altitude  and 
scattering  angle.  9(2,0)  -  0(2,0)  /  *0(2,0), 


Path  length,  for  paths  of  sight  at  zenith  angles  0  to  70  degrees,  r  -  seed  As. 


Standardized  relative  spectral  response  of  fllter/cathode  combination  where  5*  is 
spectral  sensitivity  of  the  multiplier  phototube  cathode  and  7V  is  spectral 
transmittance  of  optical  fitter.  The  relative  spectral  response  values  are 
normalized  to  the  peak  value. 


Total  volume  scattering  coefficient  as  determined  at  altitude  2.  This  property  may 
be  defined  by  the  equations 


2(2)  -  f<r(z,0)d O  -  *2(2)  +  ys(z), 


In  the  absence  of  atmospheric  absorption,  the  total  volume  scattering  coefficient  is 
numerically  equal  to  the  attenuation  coefficient. 

Volume  scattering  coefficient  for  Mie  l.e,  particulate  or  droplet,  scattering  at 
altitude  z, 


Volume  scattering  coefficient  for  Rayleigh  l.e,  molecular,  scattering  nt  altitude  z. 


Absolute  temperature  in  degrees  Kelvin. 


f,(z,0) 


Radiance  transmittance  as  determined  at  altitude  z  for  a  path  of  sight  of  length  r 
at  zenith  angle  9  (formerly  referred  to  as  "beam"  transmittance).  This  property  is 
Independent  of  azimuth  in  atmospheres  having  horizontal  uniformity.  It  is  always 
the  same  for  the  designated  path  of  sight  or  its  reciprocal. 


Tm(z,9) 


Radiance  transmitance  for  the  path  of  sight  at  zenith  angle  9  from  out  of  the 
atmosphere  to  the  altitude  z. 
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Symbol 

Units 

Quantity 

f4,(z) 

none 

Optical  thickness  r4I(z)  —  f  a(z)dz. 

•i 

none 

Absorption  optical  thickness  arA/(z)  -  f  a(z)dz. 

‘ , 

v'4x  (-) 

none 

Mie  optical  thickness  W/Al(z)  -  J  ys(z)dz. 

•i 

H '«<*) 

none 

Rayleigh  or  molecular  optical  thickness  «rA{(z)  -  f  Ks(z)dz. 

' , 

W 

7/m1 

Radiant  density. 

W 

W 

Watt  W  -  J/s. 

z 

m 

Altitude,  usually  used  as  above  ground  level. 

h 

m 

Altitude  of  any  applicable  object  or  target. 

a(z) 

Volume  attenuation  coefficient  as  determined  at  altitude  z.  a(z)  -  a(z)  +  s(z). 

0 


deg 


Symbol  for  scattering  angle  of  flux  from  a  light  source.  It  Is  equal  to  the  angle 
between  the  line  from  the  source  to  any  unit  scattering  volume  and  the  path  of  a 
ray  scattered  off  this  direct  line.  See  illustration. 


A  none  Symbol  to  indicate  incremental  quantity  and  used  with  r  and  z  to  indicate  small, 

discrete  increments  in  path  length  r  and  altitude  z. 


8\  nm  Effective  passband  (formerly  designated  "response  area”)  for  a  Altered  sensor  is 

M 

defined  as  8*  -  j  $KTK  dk. 

"o 
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Symbol 

Units 

Quantity 

«(z) 

ffVm 1 

Scalar  irradiance.  This  may  be  defined  as  the  radiant  flux  arriving  at  a  point,  from 
all  directions  about  that  point,  at  altitude  z  (Tyter  and  Preisendorfer  (1962)). 
«(z)-  JL  (z,»', <t>')dn'. 

dt(z) 

W/m* 

Diffuse  scalar  irradiance  dt(z)  -  «(z)  -  ,«(z>. 

/t(z) 

W/m1 

Sun  scalar  irradiance  at  altitude  z  t«(z)  -  ,«(<»)  r„(z  ,(>,), 

C 

m 

Radius  of  the  earth. 

e 

deg 

Symbol  for  zenith  angle.  This  symbol  Is  usually  used  as  one  of  two  coordinates  to 
specify  the  direction  of  a  path  of  sight. 

9' 

deg 

Symbol  for  zenith  angle  usually  used  as  one  of  two  coordinates  to  specify  the 
direction  of  a  discrete  portion  of  the  sky. 

«, 

deg 

Zenith  angle  of  the  sun. 

X 

nm 

Symbol  for  wavelength. 

X 

nm 

Mean  wavelength  Is  defined  as  X  -  J  X  SI 77  d\  /  CX. 

0 

sr1 

Ratio  of  solar  almucantar  sky  radiance  to  scalar  sun  irradiance  and  airmass. 
/*<*,/ 3)  -  Ll(z,6„p)/  [  ,«(z)mM. 

t(l) 

IV/m5 

Net  irradiance.  ((z)  -  E(z,d)  -  E{z,u). 

p(z) 

kg/m* 

Density  at  altitude  z. 

tr1*  t 3) 

m~lsr~l 

Symbol  for  volume  scattering  function.  Parenthetical  symbols  are  z  to  designate 
altitude  and  /3  to  designate  the  scattering  angle  from  a  source. 

tr(z,/3)/s(z) 

sr-1 

Normalized  volume  scattering  function.  This  may  be  defined  by  the  equation 
|l  <t(z,i 3)  /  j(i)]«fn  -  i. 

sr 

Optical  thickness  function.  tAi(z,|3)  -  JV(z,0)rfc. 
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Symbol 

Units 

Quantity 

* 

deg 

Symbol  for  azimuth.  The  azimuth  is  the  angle  In  the  horizontal  plane  of  the 
observer  between  a  fixed  point  and  the  path  of  sight.  The  fixed  point  may  be  for 
example,  true  North,  the  bearing  of  the  sun,  or  the  bearing  of  the  moon.  This 
symbol  is  usually  used  as  one  of  two  coordinates  to  specify  the  direction  of  a  path 
of  sight. 

deg 

This  symbol  for  azimuth  is  usually  used  as  one  of  two  coordinates  to  specify  the 
direction  of  a  discrete  portion  of  the  sky. 

* 

deg 

Angular  solar  radius  at  true  earth-to-sun  distance. 

$ 

deg 

Angular  solar  radius  at  mean  solar  distance. 

n 

sr 

Symbol  for  solid  angle,  l!or  a  hemisphere,  A  -  2ir  steradians.  For  a  sphere, 
a -Aw  steradians. 

Ctf(z) 

none 

Single  scattering  albedo  u(z)  —  a (z)  /  a(z). 
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Table  C.l.  Notation*!  equivalencies. 


Chandrasekhar  (1960) 


Visibility  Laboratory 


Symbol 


Definition 
Radiant  energy 
Net  flux 

Specific  intensity 
Average  intensity 


Mass  emission  coefficient 


Symbol 


Hz, 9 ,+) 

t/4n 


Definition 
Irradiance 
Net  irradiance 
Radiance 
Average  radiance 


Path  function/density 


Mass  emission  coefficient  x  density 

Mass  emission  coefficient  due 
to  scattering 

Mass  absorption  coefficient 
Mass  absorption  coefficient  x  density 


L*(z,9,$)  Path  function 


p(cor0)  Phase  function 

$  Thickness 

3  Source  function 


Integrated  energy  density 

Single  scattering  albedo 

Scattering  angle 

Polar  angle 

cosine  of  polar  angle 

Frequency 

Density 

Normal  optical  thickness 

Azimuth 

Solid  angle 


U(z,6,4>) 


4tror(/S)/a 


Path  function  due  to  scattering 
divided  by  the  density 

Attenuation  coefficient  divided 
by  the  density 

Attenuation  coefficient 

Four  rr  times  volume  scattering 
function/attenuation  coefficient 

Path  length 

Equilibrium  radiance 

Radiant  density 

Single  scattering  albedo 

Scattering  angle 

Zenith  angle 

cosine  of  zenith  angle 

Inverse  wavelength 

Density 

Optical  thickness 
Azimuth 
Solid  angle 
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APPENDIX  D 


FASCAT  PROGRAM 


STRUCTURE  CHART 
SUBROUTINE  SUMMARY 


INPUT/OUTPUT  FILES 


USER  NOTES 


FLOW  DIAGRAM 


■  •'-YAV.v 
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C*«***********************  TITLE  *****************%*********** 
C  FASCAT » . (FORTRAN  PROGRAM) 


A  FAST  ATMOSPHERIC  SCATTERING  MODEL  FOR  CALCULATING 
APPARENT  BACKGROUND  AND  TARGET  RADIANCE  FIELDS. 
RESEARCH  BY  WAYNE  S.  HERING 
VISIBILITY  LABORATORY 
SCRIPPS  INSTITUTION  OF  OCEANOGRAPHY 
UNIVERSITY  OF  CALIFORNIA  AT  SAN  DIEOO 
LA  JOLLA*  CALIFORNIA  92093 


METHOD 

C  THE  FASCAT  COMPUTER  CODE  CALCULATES  THE  SPECTRAL  RADIANCE 

C  OF  DISTANT  OBJECTS  AND  BACKGROUNDS  FOR  DESIGNATED  SLANT 

C  PATH  IN  THE  ATMOSPHERE . THE  APPARENT  RADIANCE  18  THE  SUM 

C  OF  RESIDUAL  BACKGROUND  OR  TARGET  RADIANCE  PLUS  THE  PATH 

C  RADIANCE  GENERATED  BY  SCATTERING  OF  LIGHT  REACHING  THE 

C  PATH  FROM  THE  SUN  AND  THE  SURROUNDING  SKY  AND  TERRAIN. 

C  SEVERAL  INPUT  OPTIONS  ARE  AVAILABLE  FOR  THE  SPECIFICATION 

C  OF  THE  VERTICAL  PROFILES  OF  AEROSOL  PROPERTIES,  CLOUD 

C  INFORMATION  IS  ENTERED  BY  SPECIFICATION  OF  THE  CLOUD  TYPE 

C  AMOUNT .ALTITUDE  OF  BASE  AND  TOP  AND  THE  RELATIVE  OPTICAL 

C  THICKNESS  ( AVERAGEiTHIN  OR  THICK).  FOR  DESIGNATED  SENSOR 

C  AND  TARGET  ALTITUDES.THE  APPARENT  CONTRAST  IS  CALCULATED 

C  FOR  CLOUD  FREE  SLANT  PATHS.UHICH  MAY  DR  MAY  NOT  TRAN3- 

C  VERSE  A  FRACTIONAL  CLOUD  LAYER.  OBJECTS  IN  SUNLIGHT  OR 

C  SHADOW  MAY  BE  VIEWED  AGAINST  SKY.  CLOUD  OR  TERRAIN  BACK- 

C  GROUNDS. 
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TRIDAO 


] 


SKYRAD 


NSTAR 
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C*********************  PARAMETERS  ***************************** 
C 

c  SUBROUTINE  ADJUST  < I  OPT i INDEXA» INDEXB ) 

C  ASSIGN  ARRAYS  OF  ATMOSPHERIC  LAYERS 

C  IOPT  -  ASSIGNMENT  CHOICES 

C  INDEXAr INDEXB  -  MOVE  INDEXB  INTO  INDEXA 

C 

C  SUBROUTINE  CALC  <TAUiM» IPRINT»N»IBTORE> 

C  DRIVES  FLUXES  SUBROUTINE 

C  TAU-  OPTICAL  DEPTH 

C  N«  INDEX  ON  LAYER/SUBLAYER 

C  IPRINT  *  PRINT  OPTION  <0-N0>l-YEB) 

C  N  -  TOTAL  ATMOSPHERIC  LAYERS 

C  ISTORE...O-NO'1-YES 

C 

C  SUBROUTINE  CLDFR< CLR t A » B t AB» ANS» IOPT » M > 


C  COMBINE  DELTA  EDDINGTON  VARIABLES  FOR  DIFFUSE 

C  RADIANCE  CALCULATION 

C  CLR  -  CLEAR  SKY  CONTRIBUTION 

C  A  -  FIRST  CLOUD  LAYER  CONTRIBUTION 

C  B  ■  SECOND  CLOUD  LAYER  CONTRIBUTION 

C  AB  «  ALL  CLOUD  CONTRIBUTIONS 

C  ANS  »  COMBINED  RADIANCE 

C  IOPT. . .1-ABOVE  TOP  CLOUD>  2-  BELOW  TOP  CLOUD 

C  M  -  INDEX  ON  LAYER 

C 

C  SUBROUTINE  DELTED<N» IPRINT r I STORE) 

C  EXPAND  AND  PRINT  DELI A-EDDINGTON  PRODUCTS 


N  «  TOTAL  ATMOSPHERIC  LAYERS 
IPRINT  ■  PRINT  OPTION  <  0«NO»  1-YES) 
ISTORE  -  STORE  OPTION  <  0«N0»  1-YES) 
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C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  ERRMOGUERR) 

HANDLES  FATAL  ERROR  MESSAGES 
JERR  -  POINTS  TO  FATAL  ERROR 

SUBROUTINE  FLUXE8<N»TAU> 

DETERMINE  IRRADIANCE  PROF ILE8  FOR  N  LAYER  ATMOSPHERE 
N-  NUMBER  OF  ATM08PHERIC  LAYERS 

TAU-  OPTICAL  DEPTH  FOR  WHICH  FLUXES  ARE  CALCULATED 
ERIC  SHETTLE  AFGL/OPA  HANSCOM  AFB  MA  01731 

SUBROUTINE  FUNCTION  QL<Q) 

GL  *  ASYMMETRY  FACTOR  FOR  COMBINED  RAYLEIGH  AND 
AEROSOL  SCATTERING 
WHERE  Q  IB  THE  SCATTERING  RATIO 

(RATIO  OF  TOTAL  1  RAYLEIGH  SCATTERING  COEFFICIENTS) 

r 

SUBROUTINE  INDATA< ISEAi LAMDA » IPRINT »N» ISTORE » BRO > 
SUBROUTINE  TO  INPUT  USER  DATA 
ISEA  «  REFLECTANCE  OPTION  <0-N0»l-YE8> 

LAMDA  «  WAVELENGTH 

IPRINT  -  PRINT  OPTION  <0-N0»  1-YES) 

N  »  TOTAL  ATMOSPHERIC  LAYERS 
ISTORF.  -  STORE  OPTION  (0-N0»  1-YES) 

BRO-SURFACE  RALEIGH  SCATTERING  COEFF  (550nm) 
SUBROUTINE  NSTAR ( I  TIN » IBALT * IUD» ANS» T> IT r TT * IPHI ) 

COMPUTES  RADIANCE  GENERATED  OVER  PATH  FROH  T  *  N  OF 


C 

C 


LAYERS 

ITIN-  TOP  ALTITUDE  LEVEL 
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C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


IBALT-BOTTOM  ALTITUDE  LEVEL  AND  LATER 
IUD  ■  1  FOR  UP  THETA  (0-83) 

2  FOR  DOWN  THETA  (90-180) 

ANS  "  RESULT 

T  -  TRANSMITTANCE  FOR  RANGE 
IT  >  INDEX  ON  THETA  (VIEW  ANGLE) 

TT  ■  DIFFUSE  TRANSMITTANCE 
IPHI  *  INDEX  ON  AZIMUTHAL  ANGLE 

SUBROUTINE  OPTIC(LAMDA» N. BRO » KCLOUD. IPRINT » ISTORE ) 
SUBROUTINE  TO  COMPUTE  OPTICAL  THICKNESS 
LAMDA  -  WAVELENGTH 
N  -  TOTAL  ATMOSPHERIC  LAYERS 
PRO-  SURFACE  RALEIGH  SCATTERING  COEFF  (SSOnm) 
KCLOUD  >  CLOUD  LAYER  INDEX 

IPRINT  »  PRINT  OPTION. ♦.i«YESrOxNO 

* 

I8T0RE  -  STORE  OPTION  ... 1»YE3.0»NQ 


SUBROUTINE  PHASEF (VEXT . WAVLEN. NMODEL » RHUMID » ALT .MPHASE. 
EXTG.8CATTC.SSALB.SRATIQ.PHANGL.PHFUNC .NPHASE. 
ASSYM* IERR ) 

SUPPLIES.  VIA  LOOK-UP  TABLES.  THE  EXTINCTION  AND 
SCATTERING  COEFFICIENTS.  SINGLE  SCATTER  ALBEDO. 
SCATTERING  RATIO.  SINGLE  SCATTER  PHASE  FUNCTION 
AND  A3SYMETRY  PARAMETER  FOR  ANY  OF  THE  LOWTRAN  3 
AEROSOL  MODELS.  FOR  WAVELENGTHS  BETWEEN  0.35  AND 


B-6 


APPENDIX  B 


C 

C 

C 

C 

C 

C 

C 

C 
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c 
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c 

c 

c 

c 

c 
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1.06  MICR0-HGTER8. 

INPUT  <**MAY  BE  CHANGED  BY  SUBROUTINE  IF  OUT  OF  RANGE) 
VEXT  -  VISIBLE  EXTINCTION  COEFFICIENT  IN  1/KM. 

TOTAL  EXTINCTION*  INCLUDING  RAYLEIQH  AND 
AEROSOL  SCATTERING  AND  ABSORPTION* 

BUT  NOT  OAS  ABSORPTION  ( SEE  ABOVE). 
REQUIRED*  NO  DEFAULT*  RANGE*  IE-2  TO  1E3 
WAVLEN  -  WAVELENGTH  IN  MICRO-METERS »  RANGE*  0.55 
TO  1.06*  DEFAULTS  TO  CLOSEST  IF  OUT  OF 
RANGE. 

NMODEL-LOUTRAN  5  AEROSOL  MODEL  NUMBER* DEFAULT  *  1 

1  -  RURAL 

2  -  MARITIME 

3  -  URBAN 

4  -  TROPOSPHERIC 

3-7  -  NOT  USED*  DEFAULT  TO  1 

8  -  ADVECTION  FOG 

9  -  RADIATION  FQQ 

RHUMID  -  RELATIVE  HUMIDITY  IN  PERCENT.  0.0-99X* 
DEFAULT-  70X  IF  TOO  SMALL  OR  99X  IF  TOO 
LARGE. 

FOR  FOG  MODELS*  RH  INPUT  IS  IGNORED. 

ALT  ALTITUDE  IN  KM  <FOR  RAYLEIGH  CALCULATION) 
0  TO  10  KM*  DEFAULTS  TO  NEAREST  IF  OUT 
OF  RANGE. 

MPHASE  -  DIMENSION  OF  ARRAYS  PHANOLE  AND  PHFUNC* 
BELOW. 

OUTPUTS 

EXTC  -  EXTINCTION  COEFFICIENT  IN  1/KM. 

SCATTC  -  SCATTERING  COEFFICIENT  IN  1/KM. 
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8SALB  -  SINGLE-SCATTERING  ALBEDO  (0-1) 

8RATI0  -  RATIO  OF  3CATT  TO  RAYLEIOH  SCATTERING 

COEFFICIENT  »  DIMENSIONLESS  (’SCATTERING 
RATIO* ) 

PHFUNC  -  ARRAY  (1  TO  NPHASE)  OF  DIRECTIONAL  SCAT¬ 
TERING  COEFFICIENTS  IN  l/(KM*8TER). 

(SEE  IERR  > 

PHANOL  -  ANGLE8  ASSOCIATED  WITH  PHFUNC  IN  DEGREE8 
SAME  FOR  ALL  MODELS 

NPHASE  -  NUMBER  OF  ELENENTS  IN  PHANOL  AND  PHFUNC 
“M IN ( MPHASE »NZANG) (SEE  AERMOD  SUBROUTINE 

ASSYM  -  PHASE  FUNCTION  ASSYMMETRY  PARAMETER ( 0-1 ) 
IERR  -  ERROR  FLAG 

-2  -  FATAL »  MPHASE  IS  .LE «  0 

-1  -  FATAL ,  VISIBLE  EXTINCTION  NOT  SUPPLIED  OR 

OUT  OF  RANGE » 

0  -  NORMAL  EXECUTION 

1  -  WARNING »  DEFAULTS  USED  FOR  ONE  UR  MORE 

INPUTS. 

2  -  WARNING, FULL  PHASE  FUNCTION  NOT  RETURNED 

(MPHASE  TOO  SMALL) 

3  -  BOTH  ONE  AND  TWO. 

4  -  EXTINCTION.LT.PAYIEIGH,  RESET  TO  RAYLEIGH 
5-  7  -  COMBINATIONS  OF  1,2,4 

LTC  JOHN  D,  MILL,  AFGL/OPA ,  HANSCOM  AFB »  MA  01731 

SUBROUTINE  PLTCLD(N) 

DETERMINE  MULTIPLE  SCATTERED  RADIANCE  COMPONENT  FOR 
PARTLY  CLOUDY  SKIES 
N  ■  TOTAL  ATMOSPHERIC  LAYERS 
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SUBROUTINE  SCATT < SVAR1 » SVAR2 » BRO p I > 

DETERMINE  8CATTERIN0  RATIO  OF  OVERCAST  CLOUD  LAYER 
SVAR1  -  SCATTERING  VARIANCE 
SVAR2  *  SCATTERING  VARIANCE 
BRO  -  RAYLEIGH  SCATTERING  COEFFICIENT 
I  »  INDEX  ON  CLOUD  LAYER 

SUBROUTINE  SENSITtJUD'N' JRT  > JTOT) 

FINDS  SENSOR  MATCH  WITH  BASE  ALTITUDE 
JUD  . .»  1“  DOWN  PATH'  2-UP  PATH 
N  ■  TOTAL  ATMOSPHERIC  LAYERS 
JRT  ...  1-RADIANCE'  2-TARGET  CONTRAST 
JTOT  -  TOTAL  SETS 

SUBROUTINE  SKYIT <N» IS£A» JUD » JTOT » JRT > 

SUBROUTINE  FOR  AZIMUTHAL  ANOLE  INDEXING  AND  FOR 
CALCULATION  OF  TRANSMITTANCE  AND  SEA  SURFACE 
REFLECTANCE 

N  »  TOTAL  ATMOSPHERIC  LAYERS 
I8EA  «  REFLECTANCE  OPTION  <0*N0»  1-YES) 

JUD...  1-DOWN  PATH'  2  »  UP  PATH 
JTOT...  TOTAL  SETS  OF  COMPUTATIONS 
JRT...  1  -  RADIANCE'  2  *  TARGET  CONTRAST 

SUBROUTINE  SKYPRT < I  STORE r JUD ' JTOT ' JRT» ISEA) 

PRINTS  TABLES  FOR  RADIANCE 
AND  APPARENT  TARET  CONTRAST 
ISTORE'  O-NO »  1-YES 
JUD...  1-  DOWN  PATH'  2-UP  PATH 
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JTOT « » »  TOTAL  SETS 

JRT ♦ » •  1-  RADIANCE*  2-  TARGET  CONTRAST 
ISEA. » . 1-SEA  REFLECTANCE. O-NO  SEA  REFLECTANCE 

SUBROUTINE  SKYRAD<N* ALB*  ALBP * ISEA *  JUD*  JTOT  *  JRT ) 

CALCULATES  APPARENT  BACKGROUND  AND  TARGET  RADIANCE 
FIELDS  AND  APPARENT  TARGET  C0NTRA3T 
N-TOTAL  NUMBER  OF  ATMOSPHERIC  LAYERS 
ALB  ■  AVERAGE  SURFACE  REFLECTANCE 
ALBP-LOCAL  OBJECT/BACKGROUND  REFLECTANCE 
< IF  NOT  ALB) 

ISEA  «  REFLECTANCE  OPTION  <0»N0* l-YES)/ 

JUD...  1-DOWN  PATH  2-UP  PATH 
JTOT...  TOTAL  8ET8 

JRT...  1  -  RADIANCE*  2  -  TARGET  CONTRAST 

SUBROUTINE  BPF<BETA*Q*P. IOP*Gi *G2»C> 

CALCULATES  6IN0LE  SCATTERING  PHASE  FUNCTION 
BETA  •  SCATTERING  ANGLES 
Q  «  SCATTERING  RATIO 
P  -  PHASE  FUNCTION 
IOP 

1  P-F  <  U) 

2  FOR  FOO  OR  CLOUD 

3  INPUT  TABLE  AND  G  TO  BE  SPECIFIED  BY  USER 

4  USER  SPECIFIED  01,G2*C  FOR  H.  GREENSTEIN 
CALCULATION*  AER080L  SCATTERING  PHASE  FUNCTION 
01*02*0  -  USER  I  JT  ASYMMETRY  FACTORS 

AND  WEIGHTING  FUNCTION 
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SUBROUTINE!  SPFTABIIOP.BETA.P) 

INTERPOLATE  P  FROM  IOP  OPTION  TABLE  (BETAfP) 

IOP  «  INDEX  ON  BETA  TABLE 
BETA  «  SCATTERING  ANGLES 
P  -  INTERPOLATED  PHASE  FUNCTION 

SUBROUTINE  TANIT<T»SN»L»M»  THETA* A» PHI » 

RTS.RTrlTf ISEA.JRT) 

CALCULATES  T  < TRANSMITTANCE)  AND  SNIPATH  RADIANCE) 

FOR  EACH  SUBLAYER 

T  «  TRANSMITTANCE 

SN  »  LAYER  PATH  RADIANCE 

L  »  INDEX  FOR  LAYER 

M  -  INDEX  ON  LAYER/SUBLAYER 

THETA  -  LOOK  ANGLE 

A  -  SINGLE  SCATTERING  ALBEDO  FOR  LAYER 

PHI  ■  AZIMUTHAL  ANGLE 

RTS  «  SOLAR  ZENITH  ANGLE  IN  RADIANS 

RT  »  VIEW  ANGLE  IN  RADIANS 

IT  =  INDEX  ON  VIEW  ANGLES 

ISEAr  O-NO  REFLECTANCE.  1  «  REFLECTANCE 

SUBROUTINE  TBPR < Nr  I  STORE > 

READ  IN  TABLES  TO  PRINT 
N  *  TOTAL  ATMOSPHERIC  LAYERS 
ISTORE* * .STORE  OPTION* .0-NC»  1-YES 
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C  H  “  TOTAL  ATMOSPHERIC  LAYERS 

C  LAHDA  »  WAVELENGTH 

C  BRO-SURFACE  RALEIGH  SCATTERING  COEF 

C 

C  SUBROUTINE  TRI DAG < N» SUB » DIA r SUP> CON > V > 

C  SOLVES  SYSTEM  OF  N  SIMULATNEOUS  LINES 

C  N  «  NUMBER  OF  ATMOSPHERIC  LAYER 

C  SUB r DIA r SUP » CON»SYSTEM  OF  EQUATIONS  FROM  SUB. FLUXES 

C  V  -  OUTPUT  VECTOR 

C 

C  SUBROUTINE  XRANG<X»XMIN»XMAX> 

C  DETERMINE  IF  X  IS  IN  RANGE  OF  XMIN  AND  XMAX 

C 

C  SUBROUTINE  ZENITH 

C  COMPUTATION  OF  SOLAR  ZENITH  ANGLE 

C 

c**»*  *#*n******)M  ********** * ************  ********** ******  ****** 
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Ct*tt*****t*****t  EXTERNAL  DATA  FILE**************************** 
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INPUT  FILE 

INPUT  DATA  IS  LIST  DIRECTED  DATA 


1  ♦  <  INFO  (  JA>  ?  JA-1  . . .  FORMAT  <  20A4  > 

ALPHA  TITLE 

2.  I  PRINT*  I  STORE*  I  SEA.  «  . «  . . . * « *  *  FORMAT  <31 5) 

IPRINT  ■  PRINT  OPTION  (O-NOrl-YES) 

ISTORE  -  STORE  OPTION  (O-NOi  1-YES) 

ISEA  *  REFLECTANCE  OPTION  (0«N0»  1-YES) 

3.  ALB* L AMD A r BAT  ......  .  .  .  . . . . FORMAT < 3F10 . 4  ) 

ALB  -  AVERAGE  SURFACE  REFLECTANCE 

LAMDA  ■  REPRESENTATIVE  WAVELENGTH  *  MICROMETERS ) 

BAT  »  BASE  ALTITUDE  OF  TOP  LAYER  (XM  MSL) 

4  «  FAC  . FORMAT  (FIO  *  4 ) 

EXTRATERRESTRIAL  SOLAR  IRRADIANCE  (UM .CM-2.NM-1 > 

5.  IZEN . . .  .FORMAT  <  IS  ) 

OPTION  TO  INPUT  SOLAR  ZENITH  ANGLE 
<0»NO»  1-YES) 

6.  THETAS* JULIAN r XL AT  <IF  IZEN-1 > *FQRHAT(F10.3r I5rF10,3> 

THETAS  -  SOLAR  ZENITH  ANGLE  < DEGREES  AND  TENTHS) 

JULIAN  -  JULIAN  DATE 

XLAT  -  LATITUDE  (DEGREES  AND  TENTHS) 

7.  HOUR*  JULIAN  *  XLAT ( IF  IZEN-Q) . . . ,F0RMAT(F10.3» I3»F10.3> 

HOUR  -  LOCAL  TIME  (HOURS  AND  TENTHS) 

JULIAN  -  JULIAN  DATE 

XLAT  -  LATITUDE  (DEGREES  AND  TENTHS) 
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8.  NTHETA . . . . . FORMAT  (15) 

TOTAL  NUMBER  OF  VIEW  ANGLES 

IF  NTHETA  *  Of  USE  DEFAULT  VIEU  ANGLES 

9.  < THETA <JA)f  JA®1 f NTHETA) . FORMAT ( 1  OF  10 , 3 ) 

IF(NTHETA.GT.O) 

INDIVIDUAL  VIEW  ANGLES 

10.  NLAY . . . FORMAT  <  15) 

TOTAL  NUMBER  OF  INITIAL  LAYERS  OF  ATMOSPHERE 

1 1 .  ISPF  ( ILAY) . . . . . . .  FORMAT  ( 15 > 


NOTE.. REPEAT  INPUT  RECORDS  11  t  12  NLAY  TIME8 
AEROSOL  TYPE  OF  ATMOSPHERIC  LAYER 

1  «  INPUT  ZLfQfU  < CALCULATES  PHASE  FUNCTIONS) 

2» INPUT  ZLfQfW(SPECIFIES  CLOUD  l  FOO  PHASE  FUNCTION) 
31 f 32 f 33  -  ZLfOfUfG  (ENTER  PHASE  FUNCTION  TABLES) 

A  «  ZLfOfWfQI fG2fC  (H.Q.  PHASE  FUNCTION  CALCULATION) 
5  »  ZLfVEXTfRHUMIDiNMODEL  (LOWTRAN  MODELS) 

12.  RECORD  12«INPUT  ACCORDING  TO  AEROSOL  TYPE ( RECORD  11) 
TYPE  1  *  ZL< ILAY ) f  Q< ILAY) fU( ILAY) . • » .FORMAT ( 3F10.3) 
TYPE  2  ■  ZL< ILAY) fQ< ILAY) fW( ILAY) . . . .FORMAT ( 3F10.3) 
TYPE  31 f 32 f ETC. . , ZL < ILAY) f 0< ILAY > f W< ILAY) r 3< ILAY) 
FORMAT  <  4F 1 0 ♦ 3 ) 

TYPE  4  •  ZL  ( ILAY)  fO(  ILAY)  fU(ILAY)fOKILAY)  »02(ILAY)f 
C(  ILAY ) * . .F0RMAT(7F10.3) 

TYPE  5*ZL(ILAY) fVEXT  fRHUMIDfNMODELfFORMAT(3F10.3fI5) 
ZL  «  BASE  ALTITUDE  OF  LAYER  <KM  MSL) 

U  -  SCATTERING  RATIO 
A  -  SINGLE  SCATTERING  ALBEDO 
0  *  AEROSOL  A8YMMETRY  PARAMETER 

G1  AEROSOL  ASYMMETRY  FACTOR  ( H . 0  .FUNCTION ) 
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02  -  AEROSOL  ASYMMETRY  FACTOR  (H.G .FUNCTION ) 
C  -  PARTITIONING  FACTOR  <H.G.  FUNCTION) 


NMODEL  -  LOWTRAN  AEROSIL  MODEL 
1  -  RURAL 


2  ■  MARITIME 


3  ■  URBAN 

4  ■  TROPOSPHERE 


5  -7  ■  NOT  USEDi  DEFAULT  TO  RURAL 


8  -  ADDICTION  FOG 

9  -  RADIATION  FOG 

VEXT  *  VISIBILITY  EXTINCTION  <KH-1> 
RHUMID  -  RELATIVE  HUMIDITY  (PERCENT) 


1^; 


13.  IENDT ( IF  AEROSOL  TYPE  31-33  IS  CHOSEN) . FORMAT ( 15) 

TOTAL  BETA/P  ENTRIES  FOR  TABLE 

REPEAT  RECORDS  13  AND  14  FOR  EACH  31-33  TYPE  CHOSEN 

14.  TBETArTP  < IF  AEROSOL  TYPE  31-33  CHOSEN ) rFORMAT < 2F10 . 4 > 

TBETA  “BETA  ANGLE 
TP  »  PHASE 

(REPEAT  IENDT  TIMES  FOR  EACH  31-33  AEROSOL  TYPE  ) 


IS.  NCLOUD . . . . . FORMAT ( IS ) 

TOTAL  NUMBER  OF  CLOUD  LAYERS  (MAXIMUM  OF  2) 


i\>v 


16.  ALTOP  > ALTBT » AMTFRr ICLOUD  f ICLFAC  (IF  NCLOUD . NE . 0) . . 
FORMAT  <3F10»2r2I5) 

ALTOP  •  TOP  ALTITUDE  OF  CLOUD  LAYER  (KM  MSL) 

ALTBT  -  BASE  ALTITUDE  OF  CLOUD  LAYER  (KM  MSL) 
AMTFR  ■  AMOUNT  FRACTION  OF  CLOUD  LAYER 
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ICLOUD  -  CLOUD  TYPE* 

1  ■  CIRRUS/C I RROSTRA7US 

2  »  ALTOSTRATUS/ALTOCUMULUS 

3  •  CUMULUS 

A  -  STRATUS/STRATOCUMULUS 
5  «  NIMBQSTRATUS/PRECIP IT  AT  ION 
ICLFAC. .. .RELATIVE  OPTICAL  THICKNESS 

1  ■  AVERAGE  OPTICAL  THICKNESS 

2  ■  THICK 

3  "THIN 

(REPEAT  RECORD  NCLOUD  TIMES) 

17.  NDLOOK .  . . . ......FORMAT  IS 

TOTAL  TABLES  FOR  DOWN  PATH  OF  SIGHT 

18  .  DSENS ( I ) i DTARG ( I ) * DTAREF (I)r DZNORH ( I >  * DANORM ( I > 

ID7RG< I >  ,DBREF(I ) *DZSLOP( I ) * DASLOP< I ) * IDLUM( I > 

FORMAT(3F10.2»IS*3F10.2» IS) 

IF ( NDLOOK . NE « 0) 

< I«1 i NDLOOK) 

DSENS-SENSOR  ALTITUDE  DOWNWARD  PATH  OF  SIGHT (KM  MSL) 
DTARQ»TARGET  ALTITUDE  DOWNWARD  PATH  OF  SIGHT ( KM  MSL) 
DTAREF  -  TARGET  REFLECTIVITY 

DZNORM  «  TARGET  NORMAL  ZENITH  ANGLE 
DANORM  -  TARGET  NORMAL  AZIMUTH  ANGLE 
(DEGREES  DEPARTURE  FROM  SOLAR  AZIMUTH) 

1DTRG  ■  TARGET  ILLUMINATION* 

1  -  SUNLIGHT 

2  -  cloud  shadow 

3  •  LOCAL  SHADOW 

DBREF  ■  LOCAL  BACKGROUND  REFLECTIVITY 


DZSLOP  -  BACKGROUND  NORMAL  ZENITH  ANGLE  (DEGREES) 
DA3L0P  -  BACKGROUND  NORMAL  AZIMUTH  ANGLE 
(DEGREES  DEPARTURE  FROM  SOLAR  AZIMUTH) 

IDLUM  -  BACKGROUND  ILLUMINATION?  1-  SUNLIGHT 

2*  CLOUD  SHADOW 
3  -  LOCAL  SHADOW 

19.  NULOOK . . . FORMAT  IS 

TOTAL  TABLES  FOR  UP  PATH  OF  SIGHT 

20.  USENS  ( I  >  f  UTARG  ( I )  r  UT AREF <  I  ) . UZNORM (Dr UANORM  ( I )  r 

lUTRG(I) 

IF(NULOOK.NE.O) 

FORHAT (SF10.2r2IS) 

(I  «  lr  NULOOK) 

USENS-SENSOR  ALTITUDE  UPWARD  PATH  OF  SIGHT  (Kh  MSL) 
UTARG-TARGET  ALTITUDE  UPWARD  PATH  OF  SIGHT  (KM  MSL) 

UTAREF  -  TARGET  REFLECTIVITY 

UZNORM  -  TARGET  NORMAL  ZENITH  ANGLE  (DEGREES) 

UANORM  -  TARGET  NORMAL  AZIMUTH  ANGLE 
(DEGREES  DEPARTURE  FROM  FROM  SOLAR  AZIMUTH) 

IUTRG  ■  TARGET  ILLUMINATION 

1  «  SUNLIGHT 

2  -  CLOUD  SHADOW 

3  -  LOCAL  SHADOW 

21.  . . . . FORMAT  15 

0  -  DEFAULT  AZIMUTHS  -  0rl80.90 
1  -  USER  INPUTS  3  AZIMUTH  ANGLE8 
(DEGREES  DEPARTURE  FROM  SOLAR  AZIMUTH) 

22.  PHINT . . FORMAT  (  3F10. 4 ) 

IF  NPHI  .NE.O  ENTER  3  AZIMUTH  ANGLES 
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Cttttt ********************************************************** 

C  NOTES  ON  INPUT/OUTPUT  DATA  FILES 

C 

C  ITEM 

C  2  IPRINT-1  PRINTS  A  SERIES  OF  TABLES  ASSOCIATED 
C  WITH  THE  CALCULATION  OF  THE  IRRADIANCE  AND  THE 

C  DIRECT  AND  DIFFUSE  RADIANCE  COMPONENTS  REACHING 

C  EACH  ATMOSPHERIC  LAYER.  THE  NUMBER  OF  TABLE!  SETS 

C  DEPENDS  ON  THE  NUMBER  OF  CLOUD  LAYERS.  FOR  EXAMPLE » 

C  IF  THE  MAXIMUM  OF  TWO  LAYERS  OF  FRACTIONAL  OR  OVER- 

C  CAST  CLOUD  ARE  PRESCRIBED!  FOUR  SETS  OF  TABLES  OF 

C  4  TABLES  EACH  ARE  PRINTED.  THE  FIRST  SET  INCLUDES 

C  ONLY  THE  UPPER  CLOUD  LAYER  IN  THE  CALCULATIONS!  THE 

C  SECOND  SET  THE  LOWER  CLOUD  LAYER!  AND  THE  THIRD  SET 

C  INCLUDES  BOTH  CLOUD  LAYERS,  THE  LAST  SET!  OR  FOURTH 

C  SET  IN  THIS  CASEi  REFERS  TO  THE  CLOUD  FREE  PORTION 

C  OF  THE  ATMOSPHERE  EXCEPT  FOR  INSTANCES  WHERE  OVER- 

C  CAST  IS  PRESCRIBED  FOR  A  GIVEN  CLOUD  LAYER. 

C 

C  THE  FIRST  TABLE  IN  EACH  SET  LISTS  THE  RAYLEIGH  AND 

C  TOTAL  OPTICAL  THICKNESS  FOR  EACH  SUBLAYER.  THE 

C  SECOND  TABLE  LISTS  THE  VERTICAL  BEAM  TRANSMITTANCE! 

C  THE  SINGLE  SCATTERING  ALBEDO  AND  THE  COMBINED 

C  RAYLEIGH  PLUS  AEROSOL  ASYMMETRY  FACTOR  FOR  EACH 

C  SUB-LAYER.  THE  THIRD  TABLE  LISTS  THE  TOTAL  OPTICAL 

C  DEPTH  PROFILE  FOR  THE  ACTUAL  AND  MODIFIED 

C  DELTA-EDDINOTON  ATMOSPHERE,  THE  FOURTH  TABLE  IN 

C  EACH  SET  LISTS  THE  COMPONENTS  OF  THE  IRRADIANCE 

C  PROFILES  A8  DETERMINED  BY  THE  DELTA-EDDINOTON 

C  COMPUTER  CODE.  B.lf, 


I8T0RE  -  1  STORES  DATA  IN  UNFORMATTED  BINARY  FORM 
OUTPUT  HANDLING  IS  SYSTEM  DEPENDENT 
LOGICAL  UNIT  NUMBER  IN  PROORAM  FASCAT  IS  4 

ISEA  -  1  PRESCRIBES  LOWER  BOUNDARY  CONDITIONS 
WHICH  CONFORM  TO  A  WIND  RUFFLED  SEA  SURFACE  WITH 
LIGHT  TO  MODERATE  WIND  SPEED.  8INCE  FRESNEL 
REFLECTANCE  IS  ASSUMED  AND  A  MATCHING  DOWNUELLING 
RADIANCE  AND  THE  REFLECTED  SOLAR  BEAM  COMPONENT  ARE 
REQUIRED*  THIS  OPTION  IS  SET  UP  TO  OPERATE  ONLY  IF 
THE  Or  90 r  180  DEGREE  AZIMUTHAL  ANGLE  OPTION  IS 
USED.  (SEE  ITEM  21). 

ALB  IS  USED  TO  CALCULATE  THE  CONTRIBUTION  OF  THE 
DIFFUSE  BACKGROUND  COMPONENT  OF  THE  PATH  RADIANCE* 
I.E.  RADIANCE  REFLECTED  FROM  THE  UNDERLYING  SURFACE 
AND  THEN  SCATTERED  BY  THE  ATMOSPHERE  IN  THE  DIREC¬ 
TION  OF  THE  SENSOR  OR  OBSERVER.  A  WEIGHTED  AREA 
AVERAGE  OF  THE  BACKGROUND  REFLECTIVITY  IS  REQUIRED 
GIVING  THE  M08T  WEIGHT  TO  THE  BACKGROUND  WITHIN 
0  TO  25  DEGREES  OF  THE  OBSERVATION  ANGLE.  ALB  WILL 
OFTEN  DIFFER  FROM  THE  LOCAL  OR  IMMEDIATE  BACKGROUND 
REFLECTIVITY  DBREF »  DESCRIBED  IN  ITEM  18  BELOW. 

BAT  SETS  THE  LOWER  ALTITUDE  LIMIT  OF  THE  PRIMARY 
OZONE  LAYER.  BAT  SHOULD  CORRESPOND  TO  TKE  BASE 
ALTITUDE  OF  THE  UPPER  LAYER  AS  ENTERED  IN  ITEM  11* 

AND  CORRESPOND  IN  GENERAL  TO  THE  TROPOPAUSE  HEIGHT. 
VALID  LIMIT8  FOR  THIS  ENTRY  ARE  10  TO  15  KM.  THE 
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APPENDIX  B 

OZONE  OPTICAL  DEPTH  IS  CALCULATED  AS  A  FUNCTION  OF 
WAVELENGTH*  ASSUHINO  AVERAGE  TOTAL  OZONE  A8  A 
FUNCTION  OF  LATITUDE  AND  SEASON,  OTHER  CONTRIBU¬ 
TIONS  TO  ABSORPTION  IN  THE  TOP  LAYER  SHOULD  BE 
INCLUDED  IN  THE  SPECIFIED  SINGLE  SCATTERING  ALBEDO 
FOR  THE  LAYER  (ITEM  11). 

THE  EXTRATERRESTRIAL  SOLAR  IRRADIANCE  MAY  BE 
ENTERED  IN  ANY  DESIRED  UNIT8  ♦  THE  RESULTANT  RADI¬ 
ANCE  WILL  THEN  BE  IN  THE  SAME  UNITS  /SR, 

A  NORMAL  SET  OF  UNITS  FOR  IRRADIANCE  WOULD  BE 
W*10E-6/NM*CM#*2,  IT  SHOULD  BE 
NOTED  THAT  APPARENT  CONTRAST  AND  CONTRAST  TRANS¬ 
MITTANCE  ARE  INDEPENDENT  OF  FAC, 

SOLAR  ZENITH  ANGLE  IS  AUTOMATICALLY  CALCULATED  FROM 
TIME.DATE  AND  LATITUDE  DATA, 

IF  THE  OPTION  TO  ENTER  A  SPECIFIC  SOLAR  ZENITH  ANGLE  IS 

CHOSEN  IN  5*  IT  IS  ENTERED  IN  ITEM. 

TH£  LATITUDE  AND  JULIAN  DATE  ARE  ALSO  REQUIRED  FOR 

DETERMINATION  OF  THE  AVERAGE  TOTAL  OZONE  AMOUNT 

< SEE  ITEM  3>, 

IF  LOCATION  IS  IN  SOUTHERN  HEMISPHERE*  KEEP  THE 
LATITUDE  ENTRY  POSITIVE*  BUT  ADD  183  TO  JULIAN  DATE. 

IT  SHOULD  BE  NOTED  THAT  ONLY  ROUGH  CORRECTIONS  FOR 
REFRACTION  AND  EARTH ' S  CURVATURE  ARE  MADErSO  THAT 
RADIANCE  CALCULATIONS  FOR  OBSERVATION  ANGLES  CLOSE  TO 
THE  HORIZON  ARF  SUBJECT  TO  ERRORS  ASSOCIATED  WITH  THOSE 
UNRESOLVED  EFFECT8. 
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ONE  HAY  INTRODUCE  A3  HANY  LAYERS  TO  DEFINE  THE 

OPTICAL  PROPPERTIES  OF  THE  ATMOSPHERE  AS  WARRANTED 
BY  THE  COMPLETENESS  AND  ACCURACY  OF  THE  AVAILABLE 

OBSERVATIONS  AND  FORECASTS*  THE  OPTICAL  PROPERTY 
PROFILES  FOR  THE  CLOUD  FREE  AREAS  < ABOVE*  BELOW* 

AND  BETWEEN  CLOUDS)  ARE  SET  UP  IN  ITEMS  10  THRU  14. 
CLOUD  LAYERS  ARE  SUPERIMPOSED  IN  ITEM8  15  AND  16 
AND  SENSOR  AND  TARGET  LEVELS  ARE  SPECIFIED  IN  ITEMS 
17-20.  EVEN  IN  THE  CASE  OF  AN  OVERCAST  SKY* 

SET-UP  A  CLOUD  FREE  ATMOSPHERE  INITIALLY  FOR  THE 
COMPLETE  VERTICAL  PROFILE.  FOR  OVERCAST  LAYER  (NO 
CLOUD  FREE  PATH)  THE  CLOUD  INFORMATION  IN  ITEM  16 
FOR  THAT  LAYER  WILL  SUPERSEDE  THE  CLEAR  ATMOSPHERE 
PROFILE  INITIALLY  POSTULATED  IN  ITEM  11  FOR  THAT 
ALTITUDE  INTERVAL.  A  CLEAR  ATM08PHERE  BOUNDARY  IN 
ITEM  11  8H0ULD  NOT  BE  CHOSEN  INTERIOR  TO  A  CLOUD 
LAYER  AS  SPECIFIED  IN  ITEM  16. 

ISPF(ILAY)  PRESCRIBES  THE  METHOD  USED  TO  INTRO¬ 
DUCE  THE  OPTICAL  PROPERTIES  FOR  EACH  LAYER  ( BASE 
ALTITUDE  ZL)  IN  A  PLANE  PARALLEL  HORIZONTALLY  HOMO¬ 
GENEOUS  ATMOSPHERE.  FIVE  METHODS  (OPTIONS)  ARE 
AVAILABLE  AS  DESCRIBED  BELOW.  OPTICAL  PARAMETERS  TO 
BE  SPECIFIED  AREt 

A  -  SINGLE  SCATTERING  ALBEDO. 

RATIO  OF  TOTAL  VOLUME  SCATTERING  COEFFICIENT 

(RAYLEIGH  AND  AEROSOL)  TO  TOTAL  EXTINCTION 

COEFFICIENT  ( RAYLEIOH  AND  AEROSOL  SCATTERING 

PLUS  AEROSOL  ABSORPTION). 
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Q  -  SCATTERING  RATIO. 

RATIO  OF  TOTAL  VOLUME  SCATTERINO  COEFFICIENT 
TO  RAYLEIGH  SCATTERING  COEFFICIENT. 

P(BETA)  -  AEROSOL  SINGLE  SCATTERING  PHA8E  FUNCTION. 

G  -  AEROSOL  PHASE  FUNCTION  ASYMMETRY  FACTOR. 


THE  OPTIONS  FOR  ENTRY  OF  THESE  VARIABLES  ARE  I 

1)  DATA  ENTRIES  FOR  EACH  LAYER  ARE  Q  AND  U  ONLY. 
P(BETA)  AND  0  ARE  DETERMINED  BY  3-TERM 
HENYEY-GREENSTEIN  FUNCTIONS  WHOSE  COEFFI¬ 
CIENTS  01 r  G2  AND  C  ARE  PRESCRIBED  IN  THE 
FA8CAT  CODE  A6  A  FUNCTION  OF  Q. 

2)  DATA  ENTRIES  ARE  0  AND  W  ONLY.  OPTION  MAY 
BE  USED  A8  AN  ALTERNATE  GENERAL  FORMAT  FOR 
FOG/CLOUDS.  P  AND  G  FOR  AVERAGE  DROP  SIZE 
DISTRIBUTION  ARE  RETURNED  FROM  MEMORY. 

3)  DATA  ENTRIE8  ARE  Q.  U.  G  AND  P  (SEE  ITEMS 
13  AND  14).  OPERATOR  INPUTS  ALL  VARIABLES 
INCLUDING  LOOK-UP  TABLE  FOR  SINGLE  SCATTER¬ 
ING  PHASE  FUNCTION  (SEE  ITEM  13). 
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4)  DATA  ENTRIES  ARE  Q*  U*  AND  THE  HENYEY- 
QREENSTE1N  PHASE  FUNCTION  COEFFICIENTS »  01* 

02*  AND  C. 

5)  DATA  ENTRIES  ARE  HADE  FOR  EACH  LAYER 
THROUGH  REFERENCE  TO  LOWTRAN  AEROSOL  MODELS. 

USINO  SUBROUTINES  BY  J.  HILL/AFL.  OPERATOR 
ENTERS)  (A)  EXTINCTION  COEFFICIENT  (SSO  MM) 

(RAYLEIGH  AND  AEROSOL  SCATTERING  PLUS 
AEROSOL  ABSORPTION  BUT  NOT  GASEUS  ABSORP¬ 
TION)  I  (B)  RELATIVE  HUMIDITY)  AND  (C)  LOWTRAN 

AER080L  MODEL  (RURAL*  MARITIME*  URBAN* 
TR0P08PHERE*  ADVECTION  FOG*  OR  RADIATION  FOO). 
REQUISITE  OPTICAL  PROPERTIES  ARE  RETURNED  BY 
MODEL  ALGORITHMS  AND  TABLE  INTERPOLATION. 

13114  IF  ISPF(ILAY)  OPTION  IS  31*  32*  OR  33  (OPTION  3 
ABOVE)*  BETA/P  TABLE  ARE  ENTERED  HERE. 

13  A  MAXIMUM  OF  2  CLOUD  LAYERS  MAY  BE  ENTERED.  IF  MORE 
LAYERS  ARE  OBSERVED  OR  REPORTED*  APPROXIMATE  THE 
SKY  CONDITION  WITH  A  2-LAYER  REPRESENTATION  OF  CLOUD 
COVER. 

16  FOR  EACH  CLOUD  LAYER*  ENTER  THE  BASE  AND  TOP  ALTI¬ 
TUDE  (KM  MSL )  AND  THE  FRACTION  OF  SKY  COVERED  BY  THE 
CLOUD  LAYER  IN  TENTHS.  THE  CLOUD  TYPE  (ICLOUD)  IS 
IDENTIFIED  BY  ONE  OF  THE  STANDARD  CLOUD  FORMS 
( CIRRUS/C IRRQ8TRATUB*  AL TOST RAT US/ ALTOCUMULUS  * 
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C  CUMULUS  *  STRATUS/$TRATQCUMULUS  »  AND  NI MB0STRATU8/ 

C  PRECIPITATION.  THE  RELATIVE  OPTICAL  THICKNESS 

C  ( ICLFAC)  HAS  3  OPTIONS!  AMERAQE  *  THICK  ►  AND  THIN. 

C  MODEL  ALGORITHMS  RETURN  THE  APPROPRIATE  OPTICAL 

C  DEPTH  FROM  ICLOHD  AND  ICLFAC  WITHOUT  DIRECT  REFER- 

C  ENCE  TO  THE  GEOMETRIC  THICKNE8S  OF  THE  CLOUD. 

C 

C  17118  INTRINSIC  PROPERTIES  OF  LOCAL  SLOPING  BACKGROUNDS 
C  AND  3~DIMENSX0NAL  TARGETS  FOR  SELECTED  DOWNWARD 

C  LOOKING  PATHS  OF  SIGHT  ARE  ENTERED  IN  ITEM  18. 

C  THE  TOTAL  NUMBER  OF  COMBINATIONS  OF  SENSOR-TARGET 

C  ALTITUDES  AND  TARGET-BACKGROUND  CONDITIONS  TO  BE 

C  DEALT  WITH  ARE  ENTERED  IN  ITEM  17, 

C 

C 

C  **##  IMPORTANT  NOTE  *.*«* 

C 

C  FOR  EXTRATERRESTRIAL  SENSOR  POSITION*  DESIGNATE 

C  SENSOR  ALTITUDE  (DSENS)  AS  100  KM. 

C 

C 

c  ****  IMPORTANT  NOTE  **** 

C 

C 

C  THE  ZENITiH  ANGLE  OF  THE  TARGET  NORMAL  AND  THE  ZENITH 

C  ANGLE  OF  THE  LOCAL  BACKGROUND  NORMAL  ARE  ENTERED 

C  DIRECTLY  WITHOUT  REGARD  TO  THE  SOLAR  ZENITH  ANGLE . 

C  HOWEVER* THE  AZIMUTHAL  ANGLES  OF  THE  TARGET  NORMAL  AND 
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THE  LOCAL  BACKGROUND  NORMAL  ARE  ENTERED  A3  DEPARTURES 
REPEAT  DEPARTURES  FROM  THE  SOLAR  AZIMUTHAL  ANGLE » 
WITHOUT  REGARD  TO  SIGN. 


C 
C 
C 
C 
C 

C  THE  INHERENT  LOCAL  BACKGROUND  AND  TARGET  RADIANCES 

C  ARE  DETERMINED  UNDER  THE  ASSUMPTION  OF  DIFFUSE  OR 

C  LAMBERTIAN  REFLECTANCE  OF  THE  IRRADIANCE  REACHING 

C  THE  OBJECT  AND  BACKGROUND  SURFACE  FROM  ALL  DIREC- 

C  TIONS.  ILLUMINATION  OF  TARGET  AND/OR  BACKGROUND 

C  CAN  BE  DESIGNATED  AS 

C 

C  1) 

C  2) 

C 

C  3) 

C 

C 

C  19120  THESE  ITEHS  DEAL  IN  A  SIMILAR  MANNER  WITH  INPUT 
C  DATA  FOR  UPWARD  LOOKING  PATHS  OF  SIGHT, 

C 

C  21122  ****  IMPORTANT  NOTE  **** 

C 

C  AS  OPPOSED  TO  ZENITH  OBSERVING  ANGLES i  THE  AZIMUTHAL. 

C  OBSERVING  ANGLES  ARE  ENTERED  AS  DEPARTURES  REPEAT 

C  DEPARTURES  FROM  THE  SOLAR  AZIMUTHAL  ANGLE. 

C 


IN  DIRECT  SUNLIGHT 

IN  CLOUD  SHADOW  WHICH  WILL  INCLUDE  DIRECT 
SOLAR  IRRADIANCE  PASSING  THROUGH  THIN  CLOUDS 

IN  LOCAL  SHADOW  WITH  NO  DIRECT  SUN  ILLUMIN¬ 
ATION. 
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C  OPTION  0  DEFAULTS  TO  AZIHUTHAL  ANGLES  OF  0  (UPSUN)r 

C  90  ( CROSSSUN)  AND  180  (DOWNSUN).  THIS  OPTION  MUST  BE 

C  USED  IF  ISEA* 1  < ITEM  2>* 

C 

C  OPTION  1  PERMITS  SELECTION  OF  ANY  3  AZIMUTHAL  VIEWING 

C  DIRECTIONS t  ENTERED  AS  THE  DEPARTURES  FROM  THE  SOLAR 

C  AZIMUTHAL  ANGLE  WITHOUT  REGARD  TO  SIGN. 

C 

C 

C*********************»*  OUTPUT  DATA  S*M*********t**«t******* 
C  TITLE 

C  AVERAGE  SURFACE  REFLECTANCE 

C  REPRESENTATIVE  WAVELENGTH 

C  BASE  ALTITUDE  OF  TOP  LAYER 

C  EXTRATERRESTRIAL  SOLAR  IRRADXANCE 

C  SOLAR  ZENITH  ANGLE 

C  JULIAN  DATE  (SOLAR  AZIMUTH  ANGLE  -  LOCAL  TIME) 

C  LATITUDE 

C  TOTAL  ATMOSPHERIC  LAYERS 

C  LAYER  NUMBER 

C  AEROSOL  TYPE 

C  ALTITUDE 

C  Qr ArGiGl *G2* CrVEXT » HUMIDITY  » MODEL  TYPE 

C  (ACCORDING  TO  AEROSOL  TYPE) 

C  TOTAL  CLOUD  LAYERS 

C  CLOUD  LAYER 

C  TYPE 

C  TOP  ALTITUDE  OF  CLOUD 

C  BASE  ALTITUDE  OF  CLOUD 
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AMOUNT  FRACTION 
RELATIVE  OPTICAL  THICKNESS 
DOWNWARD  PATH  OF  SIGHT 
SENSOR  ALTITUDE 
TARGET  ALTITUDE 
TARGET 

REFLECTIVITY 
NORMAL  ZENITH  DEGREES 
DIRECTION  AZIMUTH  DEGREES 
ILLUMINATION 
LOCAL  BACKGROUND 
REFLECTIVITY 
NORMAL  ZENITH  DEGREES 
DIRECTION  AZIMUTH  DEGREES 
ILLUMINATION 
UPWARD  PATH  OF  SIGHT 
SENSOR  ALTITUDE 
TARGET  ALTITUDE 
TARGET 

REFLECTIVITY 
NORMAL  ZENITH  DEGREES 
DIRECTION  AZIMUTH  DEGREES 
ILLUMINATION 

SUBLAYER 

Z  *  SCALE  HT  »DELTATRiDELTAT> 
MAJOR  LAYER  ASYMMETRY  PARAMETERS 

T » W  >  G » FL»FLP 


C 


SOLAR  ZENITH  ANGLE 
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C  LfZ.TAUSL.TAUPSL 

C  DELTA  EDDINGTON  RESULTS 

C  DOWNWARD  PATH  OF  SIGHT 

C  INHERENT  BACKGROUND  SURFACE  RADIANCE 

C  APPARENT  BACKGROUND  RADIANCE  TABLE 

C  INHERENT  TARGET  RADIANCE 

C  APPARENT  TARGET  RADIANCE  TABLE 

C  APPARENT  TARET  CONTRAST  TABLE 

C  UPWARD  PATH  OF  SIGHT 

C  APPARENT  BACKGROUND  RADIANCE  TABLE 

C  INHERENT  TARGET  RADIANCE 

C  APPARENT  TARGET  CONTRAST  TABLE 

C  APPARENT  TARGET  RADIANCE  TABLE 

C************t********t*tttt*tt***tt**tt*ttttt**t*t*t**t******** 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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COMMON 


/ATM0S/HS<40»3)  «  SOLAR  IRRADIANCE 

DI0(40>5)  «  DELTA  EDDINGTON  RADIANCE 
< AVERAGE  DIFFUSE) 

DI1(40»5)  -  DELTA  EDDINGTON  RADIANCE 
(VERTICAL  ASYMMETRY) 

HO < 40 )  -  DOWNWARD  SURFACE  IRRADIANCE 
P< 15»20)  -  SINGLE  SCATTERING  PHASE  FUNCTION 
RAD  =  RADIANS 
PI  -  3.1413? 

BETA ( 15 )  =  SCATTERING  ANGLES 
NB  -  TOTAL  NUMBER  BETAS 
THETA(20)  b  VIEW  ANGLES 

NT  -  TOTAL  NUMBER  OF  VIEW  ANGLES 


ABETA( 20)  *  COMPUTED  SCATTERING  ANGLES 

PHASE  ■  PHASE  FUNCTION 

UHS<40>*  SOLAR  1RRADIANCE  < BETA. LE . 25 > 

IDELTA  >  INDEX  ON  CLOUD  LAYER 
NPHI  *  0  DEFAULT  AZIMUTH  ANGLES  USED 
1  -  ENTER  3  AZIMUTH  ANGLES 
PHINT< 3)  -  IF  NPI.I.NE.O  THREE  AZIMUTH  ANGLES 
TPHIS  >  SOLAR  AZIMUTHAL  ANGLE 
/CDAT/CZL< 20)  «  LAYER  BASE  ALTITUDE  <UITH  CLOUDS) 

CO (20)  »  SCATTERING  RATIO  -  WITH  CLOUDS 

CG(20)  =>  ASYMMETRY  PARAMETER  -  WITH  CLOUDS 

CA<20)  »  SINGLE  SCATTERING  ALBEDO  -  WITH  CLOUDS 

CGI (20)  «  USER  INPUT  ASYMMETRY  FACTOR 

CG2  (20)  -  USER  INPUT  ASYMMETRY  FACTOR 

CLC ( 20 )  «  ASYMMETRY  PARAMETER 

ICSPF ( 20)  -  AEROSOL  TYPE 

/CLOUD/ICLOUDUO)  ■  CLOUD  TYPE 

ALTOP(IO)  »  TOP  ALTITUDE  OF  CLOUD  LAYER 

ALTBT < 10)  -  BASE  ALTITUDE  OF  CLOUD  LAYER 

AMTFR(IO)  »  AMOUNT  FRACTION  OF  CLOUD  LAYER 

NCLOUD-  TOTAL  NUMBER  OF  CLOUD  LAYERS 

ITOPUO)  -  CLOUD  TOP  ALTITUDE  INDEX 

IBOT < 10 )  »  CLOUD  BOTTOM  ALTITUDE  INDEX 

XCLFAC( 10) r  1 “AVERAGE  OPTICAL  THICKNESS  OF  CLOUD 

2- THICK 

3- THIN 

DTCLD<5)  -  CLOUD  TYPE  DESIGNATION 

CLDFAC  <  3 )  -  CLOUD  FACTORS  FOR  OPTICAL  THICKNESS 

PCF < 2 )  -  PARTLY  CLOUDY  FRACTION 
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/DAT/Q(20>  -  SCATTERING  RATIO 

ZL<20>  "  LAYER  BASE  ALTITUDE 
S< 40)  -  SCALE  HEIGHT 

BRO  -  SURFACE  RAYLEIGH  SCATTERING  COEFFICIENT 
DTR(40)  «  RAYLEIGH  OPTICAL  THICKNESS 
Z(  40)  -  ALTITUDE 
DT  <  40 )  *  DELTA  TAU 

FL ( 20 > "FORWARD  SCATTERING  FRACTION  MONO- 
DIRECTIONAL. 

FLB( 20) "FORWARD  SCATTERING  FRACTION* ISOTROPIC 

TAUSL<40)  ■  OPTICAL  DEPTH 

TAUPSL<40)  -  DELTA  EDDINGTON  OPTICAL  DEPTH 

THETAS  -  SOLAR  ZENITH  ANGLE 

CTHES  ■  COSINE  THETAS 

NC  »  TOTAL  THETAS  «1 

FLP<20>-  FRACTIONAL  SCATTERING  INTO  FORWARD  PEAK 
FAC  -  EXTRATERRESTRIAL  IRRADI ANCE 

/DATE/ JULIAN  ■  JULIAN  DATE 

XLAT  *  LATITUDE  IN  DEGREES 

HOUR  -  LOCAL  TIME  (TENTHS  OF  HOURS) 

/DDAT/ZZL(20)  *  ALTITUDE 

QQ(20)  «  SCATTERING  RATIO 
OG  <  20  )  =  ASYMMETRY  PARAMETER 
AA ( 20 )  "  SINGLE  SCATTERING  ALBEDO 
GG1 (20)  »  USER  INPUT  H.G .ASSYMMETRY  FACTOR 
GG2( 20 )  »  USER  INPUT  H.G.  ASSYMTRY  FACTOR 
CCC ( 20  >  "  ASYMMETRY  PARAMETER 
IISPF (20)  -  AEROSOL  TYPE 


USED  FOR  ARRAY  SWAP 
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/LUIO/ 

LUW  -  LOGICAL  UNIT  MEMBER  FOR  OUTPUT 

LUR  -  LOGICAL  UNIT  KIMBER  FOR  INPUT 

LUST  -  LOGICAL  UNIT  NUHBER  FOR  STORAGE  OF  OUTPUT 

/OPTIK/ 

HSPF (20)  -  SPF  FLAGS 

1  -  P-F<Q> 

2  «  FOG  AND  CLOUD 

3  -  USER  INPUT  OF  P/B  TABLE 

4  *  USER  INPUT  OF  Gl'G2rC 

G1 (20)  -  USER  INPUT  H.G.  ASYMMETRY  FACTOR 
G2(20)  «  USER  INPUT  H,G.  ASYMMETRY  FACTOR 
CC(20)  -  USER  INPUT  H*G»  PARTITIONING  FACTOR 
T  <  20 )  «  TAU 

BAT  ■  BASE  ALTITUDE  OF  TOP  LAYER 

/OUT/TP < 80 » 20)  ■  TRANSMITTANCE 

SN( 3 » 80 > 20 )  ■  NSTAR. » . < PHI »LAYER r THC . A) > 

SNINF( 3 » 20 » 3 )  -  SKY  RADIANCE  UP 
NSTRAD ( 3r 20 1 3 )  -  SKY  RADIANCE  DOWN 
HERIT  «  INHERENT  RADIANCE 


C 

C  /PARAM/SUB<40) » DI A<40) t SUP < 40) rCON< 40) t C< 40) »K<20> r 

C  P<20) t ALFA <20 )» BETA (20) »EKT(20) »  EKT 1 ( 20) »ET(20) 


/REFLEC/ 
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ISEAT(20)*INDEX»  VIEW  ANGLES  FOR  SEA  REFLECTANCE 
N7  =  TOTAL  SEA  REFLECTANCE  VIEW  ANGLES  CHOSEN 

/SAVE/ 

GSAVE (20)  »  ASYMMETRY  PARAMETERS 

SREF  *  SEA  SURFACE  SOLAR  REFLECTANCE 
SHS  »  SOLAR  XRRADI ANCE  AT  SURFACE 

SPHASE-PHASE  FUNCTION  FOR  REFLECTED  SOLAR  ENERGY 
SUHS  -  XRRADX ANCE  AT  SURFACE  (BETA*LE»25> 

/SPFT/TBETAH0»10)  -  CONSTANTS  FOR  BETA  TABLE 
TTP( 4 0 v 1 0 >  «  TRANSMITTANCE 
IENDT(IO)  -  INDEX  ON  B»P 

TTB( 25)  «  ALPHA  ARRAY  FOR  SPF  TABLE  TITLE 

/SSALB/U<20)  «  SINGLE  SCATTERING  ALBEDO  (SAVED) 
/TABP/ITBPRaFLAGS  » PR I NT  SELECTED  TABLE  DATA(SKYRAD) 

TABALT (2)  -  ALTITUDES  FOR  EACH  TABLE 
ALBPR  =  OBJECT/BACKGROUND  ALBEDO 

ITFLE  »  BASE  FILE  NUMBERS 
INF0(27>  »  ALPHA  TITLES  FOR  TABLES 
INFIN  2)  ■  FLAG  FOR  INFINITY  LEVEL 
NDLOOK  *  TOTAL  DOWN  PATH 

NULOOK  -  TOTAL  UP  PATH 
DTARG(IO)  «  TARGET  ALT.  DOWN  PATH 
DSENS(IO)  »  SENSOR  ALT.  DOWN  PATH 
DTAREFUO)  -  TARGET  REFLECTIVITY  DOWN 
DZNORMUO)  -  TARGET  NORMAL  ZENITH  ANGLE  DOWN 

DANORM  (10)  «  TARGET  NORMAL  AZIMUTH  ANGLE  DOWN 
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DBREF < 10)  «  LOCAL  BACKGROUND  REFLECTIVITY  DOUN 
DZSLOP(IO)  =  BACKGROUND  NORMAL  ZENITH  ANGLE  DOWN 
DASLOP( 10) "BACKGROUND  NORMAL  AZIMUTH  ANGLE  DOWN 
IDLUM<10>  *  BACKGROUND  ILLUM. ,  1 "SUNLIGHT 

2-CLOUD  SHADOW 
3  -  LOCAL  SHADOW 

IDTRG(IO) 

1- SUNLIGHT 

2- CLOUD  SHADOW 

3- LOCAL  SHADOW 

UTARGUO)  -  TARGET  ALT.  UP  PATH 
USENS(IO)  «  SENSOR  ALT.  UP  PATH 
UTAREF < 10)  -  TARGET  REFLECTIVITY 
UZNORM(IO)  ■  TARGET  NORMAL  ZENITH  ANGLE  UP 
UANORM(IO)  »  TARET  NORMAL  AZIMUTH  ANGLE  UP 
IUTRG  < 10 ) 

1- SUNLIGHT 

2- CLOUD  SHADOW 

3- LOCAL  SHADOW 

/TDIT/TD( 80, 20) -DIFFUSE  TRANSMITTANCE 

/XRAD/TS<  20)  -  OPTICAL  DEPTH  TO  BOTTOM  OF  LAYERS 

A( 20)  -  SINGLE  SCATTERING  ALBEDO  FOR  LAYERS 
G<  20 )  »  ASYMMETRY  PARAMETER  FOR  LAYERS 
ALB  -  AVERAGE  SURFACE  REFLECTANCE  (ASSUMED) 

UIN  -  COSINE  OF  SOLAR  ZENITH  ANGLE 
F-SOLAR  FLUX  CN  PLANE  NORMAL  TO  INCIDENCE 


10  -  AVERAGE  DIFFUSE  RADIANCE 

11  -  VERTICAL  ANOMALLY  DIFFUSE  RADIANCE 
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C  FDD  -  DIFFUSE  FLUX  DOWN 

C  FUP  -  DIFFUSE  FLUX  UP 

C  FDTOT  *  TOTAL  DOUNUELLXNG  IRRADIANCE 

C  FNET  =  NET  IRRADIANCE 

C  GRAD  -  IRRADIANCE  GRADIENT 

C  ENERGY  «  SPHERICAL  IRRADIANCE 

C 
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